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INTERFEROMETRIC MEASUREMENTS OF WAVE-LENGTHS 
IV. THE ACCURACY OF MEASURED WAVE-LENGTHS 
M. G. Adam and S. Nichols 
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(Received 1957 December 5) 


Summary 


Solar absorption lines are very different in character from the corre- 
sponding laboratory lines, and the wave-length accuracy obtainable for them 
is considerably less. It is moreover possible that some part of the red shift 
and change in red shift with line strength originates in the micrometer 
measurement of the broad absorption and narrow emission lines. An investi- 
gation is made of the accuracy and consistency obtainable in measuring arc 
and solar wave-lengths by the method of circular channels. It appears from 
this investigation that red shifts can be obtained which are significant to better 
than o-oor A. 





Introduction.—T he recent interest in the interpretation of solar wave-lengths 
(9, 10, 12) makes it more than ever important to consider the reliability of the 
observations against which the various theories must be tested. For a number 
of years we have been developing and using a new method of wave-length 
measurement and have been especially concerned with the reduction of systematic 
errors such as revealed themselves in earlier accepted wave-lengths. We now 
test our method, both for its potential accuracy with narrow-line spectroscopic 
sources, and also for the actual dependence we may place on our measurements 
for vacuum-arc and solar lines. 

A wave-length accuracy of one or two ten-thousandths of an angstrém can 
be achieved, in the visible region, for narrow-line laboratory sources (5), but 
in the measurement of solar wave-lengths there are several factors which set 
a limit of uncertainty considerably greater than this. The limitations are 
inherent in the nature of the solar spectrum, not only because we are concerned 
with continuous spectrum interferometry, but also because of the considerable 
intrinsic breadth of the absorption lines. 

In continuous spectrum measurements we are restricted to the use of only 
very short spacers in the interferometer (14), perhaps less than one-tenth of those 
used for the emission line sources. ‘This reduces the resolution obtainable; it 
means that the chromatic variation of phase change upon reflection at the silvered 
surfaces of the plates is important, and it magnifies the effect of any errors due to 
irregularities in the interferometer plates. The accuracy obtainable, by the 
method of circular channels, as far as this restriction and general instrumental 
limitations are concerned, we have investigated by measurements on the very 
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narrow lines emitted by a mercury isotope lamp. These measurements extend 
over the whole visible spectrum and moreover provide us with the phase change 
corrections and chromatic variation in camera focal length, which are required in 
increasing the wave-length range of our solar and arc determinations (Section 1), 

The Fraunhofer lines are so different in character from the laboratory arc 
lines with which they are compared that the question arises to what extent the 
observed solar wave-lengths, and so the red shifts and systematic variations in 
red shifts, are due to the observer’s personal judgment in the micrometer 
measurement of emission and absorption lines of very different widths. Williams 
and Middleton (15) consider that it is not possible to make micrometer 
measurements to an accuracy greater than 1/40th of the half-width of the line. 
But a solar absorption line of moderate strength (say Rowland number 8) has 
a half-width of some o-200A and on this criterion a wave-length uncertainty 
of o-oo5 A. In the green this is already half the Relativity shift. We have tried 
to answer this question in the first place by independent measurements on the 
same arc and solar plates by two observers, in fact using different machines, and 
secondly by measuring the same arc and solar lines with three different spacers 
in the interferometer. For the arc line changing the spacers chiefly changes the 
resolution, but for the solar spectrum it also alters the appearance of the plates 
quite considerably, both in the breadth of the channels in which the lines are 
measured, and in changing the amount of inter-order scattered light. ‘These 
measurements are described in Section 2 and used to form an estimate of the 
reliability of the measured red shifts. 


1. Measurements of mercury lines from an isotope lamp 


1.1. Wave-lengths obtained with short spacers and the correction for the dispersion 
of phase change.—The source in this instance was a G.E.C. water-cooled cold- 
cathode lamp containing approximately o-1 mg of Hg 198 with pure argon at 1omr 
pressure. ‘The flow of water was maintained continuously in order to avoid 
self-reversal and line broadening, and the currents were kept below 25 milliamps. 
The wave-lengths from this type of lamp are slightly longer than those from the 
electrodeless lamps, but values for both types of lamp have been very. carefully 
determined by Barrell (5). 

The procedure in measuring wave-lengths by the method of circular channels 
has already been fully described in Papers I and III (1, 3). The Oxford spectro- 
scope was used undiaphragmed in its one prism form (11) and carefully focused 
for each of the strongest mercury lines, from 5790A to 4046A, in turn. ‘The 
luminous column of the mercury lamp was imaged to fill the spectroscope 
collimator while giving full illumination over the height of the slit, as for the arc 
spectra (1). A 33 cm field lens in the spectroscope focal plane together with the 
interferometer collimating lens imaged the spectroscope collimator as a circle 
3°8 mm diameter between the interferometer plates (3). Similarly in the cadmium 
side arm a circular diaphragm before the lamp was imaged to give an identical 
and coincident light patch in the same plane. ‘The cadmium pattern was 
photographed on the plate in the image of the usual cross-shaped diaphragm, 
with two side arms as described in Paper III (3), to take care of the different 
spreading of the beams from lines at different positions in the spectroscope focal 
plane. Following the investigations described in Paper III a new interferometer 
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stand has been constructed with a different method of supporting and adjusting 
the plates to reduce as far as possible any distortion of the plates in setting up 
the interferometer. 

Fused quartz plates with a wedge angle of 2}’ were used in the etalon. ‘The 
smaller angle (cf. 15’ in the earlier work) reduces the shift of the pattern centre 
with wave-length, as does also the smaller dispersion of fused quartz compared 
with crystalline quartz. With these plates the pattern centre shift at 4046A 
leads to a wave-length change of only o-o009 A for the outermost rings of an 
outlying line, the effect on the final measured wave-length is practically negligible, 
and definitely so if a line is measured first on one side and then on the other of 
the cadmium pattern centre. ‘Three spacers were used in the interferometer, 
10 mm, 2°5 mm and 1-5 mm fused quartz rings, which covers the range of spacers 
which would normally be used for solar wave-lengths. ‘The five strongest mercury 
lines AA §790°7, 5769°6, 5460-8, 4358-3, 4046-6 A were measured with each spacer. 
Except for the pair of yellow lines which were photographed together always 
with one on either side of the pattern centre, the minimum deviation of the 
spectroscope was changed slightly between exposures so that each line was 
measured in four positions close in and far out on the red side of the pattern 
centre and close in and far out on the violet side. 

Four plates were measured for each line for the 2-5 mm and the 10 mm spacers 
and eight plates for the 1-5 mm spacer. ‘This was done in the usual way (1, 2) 
on a rotating stage micrometer. From the measured diameters fractional orders 
of interference e« were deduced for the line by the method of least squares; this 
gives incidentally a least squares value AD*? = 4f?A/t, as in Paper I (1), where 
f is the camera focal length, ¢ is the plate separation and A the wave-length. In 
the cadmium pattern the cross was used to give a least squares value of AD*? 
from which fractional orders were deduced for the three pseudo-cadmium 
‘‘lines’’, i.e. the vertical arm of the cross and the two side arms. In deducing 
wave-lengths from the line-e we assume that ¢ varies linearly with line position 
on the plate from the centre cadmium ‘“‘line’’ to the appropriate side “‘line’’: 
all measured lines lie between the centre cadmium line and one side arm. Because 
our exposures of standard (cadmium) and unknown (mercury) are simultaneous 
we need only consider the chromatic variation of the refractive index of air with 
atmospheric conditions. Pressures and temperatures were recorded for each 
exposure, and we can estimate the effect on our wave-lengths from Babcock’s 
Tables (4). The correction begins to be important at 4046 A, but under our 
conditions even here the effect is only of the order of 0-0003 A, and this has been 
ignored. 

We have calculated wave-lengths, eliminating phase change correction, from 
our results for three spacers in the following way: suppose that for any one 
spacer the order of interference for the cadmium red line (wave-length A, 
is p, (integral) and e, (fractional) and correspondingly p, and ¢, for the mercury 
line (wave-length A). If A is the phase change correction to «, corresponding 
to zero correction for the cadmium red line at 6438-4696 A then 

(po+€2—A)A= (Pi +e cas 
Le. (Potes)A=(pyteq)Agy +A.A. 
If we write k =. A, kis so small compared with the first term on the right-hand side 
(2t), that the uncertainty of A in k is quite negligible. Our unknown wave-length A 
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is then given by 
= (Pite)Aca 4 k 
Pot €e Pot €2 


The first term on the right-hand side we will call the apparent A (i.e. neglecting the 
phase change correction). Now for the plates taken with any one spacer (p. + €,) 
changes very little and k/(p.+¢€,) is constant to the accuracy with which it is 
required in this equation. ‘Thus from any number of plates we can take, for the 
spacer, a mean value of the apparent wave-length (A app.) and a mean equation 





k 
Dos = = [A app.] mean value- 


Our three spacers then provide three such mean equations at each wave-length, 
in the two unknowns A and k, for which we may make a least squares solution. 

Table I shows in column 2 the fractional part of the absolute wave-length 
so obtained and in column 3 the value of A given by this solution. Column 1 
shows Barrell’s wave-length (5) for the similar lamp. When we remember that 
Barrell’s wave-lengths have been obtained with six spacers from 2:5 cm—12°5 cm 
compared with 1-5 mm-—ro-r mm in the present instance, these results may be 
regarded as very satisfactory in that they differ from Barrell’s by well under 
o-oo1 A over the whole visible range. 


TABLE I 
Measured wavelengths and phase change corrections 


Barrell Observed 


A(A) A A C) 
4046°5716 5722 0°0496 0°0485 
4358°3377 °3379 "0309 ‘0317 
5460-7536 "7542 “0076 ‘0069 
5769°5992 "5989 ‘0051 "0054 
5790°6635 *6640 00058 0°0040 


The phase change correction A we determine in this way is not the best we can 
get from our material. ‘The most reliable values of phase change come from the 
fractional orders measured with short spacers, since de/dA= —2t/A*, any wave- 
length uncertainty has a larger effect on the « values derived from longer spacers. 
We have therefore taken Barrell’s accurate wave-lengths and used them with the ¢’s 
from the 1-5 mm spacer to establish the best value of the phase change corrections ; 
this procedure gives the corrections 6 in column 4 of Table I. These values of 5 
are higher than Barrell and ‘Teasdale-Buckell’s phase-change corrections (6), 
or the often-quoted values given by Childs (8), but not so high as the values 
deduced from Schultz’s experimental results (13). The 5, A curve is, however, 
undoubtedly peculiar to any particular pair of silver coatings and we could 
hardly expect any better agreement. Our plates were kindly coated for us (by 
evaporation) in Dr Kuhn’s laboratory. ‘The mean reflectivity at 15500 A was 
gI per cent, mean transmissivity 63 per cent. The smooth curve drawn through 
the plot of 6 and A as given above in Table I is then our definitive curve for work 
with this pair of silver films; the one in fact used for the arc and solar lines in 
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Section 2 was deduced, slightly differently, from the first four of the 1-5 mm 
plates and is negligibly different from the 6, A curve above. 

1.2. Chromatic variation in camera focal length and its application to wave-length 
determinations.—The narrow lines of the mercury 198 lamp permit the deter- 
mination of the quantity AD *? from the few measured arcs of the line interference 
pattern, but for ordinary vacuum-arc lines and especially for solar lines the values 
of AD*? so obtained would not be sufficiently accurate. We do better to obtain 
AD*? for the line from the cadmium red line pattern, knowing the ratio of the 
camera focal length at wave-length A to that at 6438-4696 A. Our measures in 
Section 1.1 enable us to determine once and for all the chromatic variation of 


(=) ; (*2-) 
A Hg A Cd red 


which is in fact the ratio (f,/fg433)”._ All the spacers should give the same value 
for this quantity; the values used were obtained from the plates taken with 
the 2-5 mm and the 10mm spacers, which we should expect to give more 
reliable values. The two sets agreed at any one wave-length to two parts in 1o* 
(see requirement in Paper I), and the mean values at each wave-length, together 
with the value previously obtained at 5085-8A (Paper 1), were used to draw a 
smooth curve showing the chromatic variation of the square of the camera focal 
length. ‘This curve has the same form as the curve supplied for the axial 
chromatic aberration of the camera lens, both giving minimum values at about 
4700 A, where (f;,/fease)” is 0°99863. 

We can now test this curve and the whole procedure we propose to use for 
our arc and solar lines by re-reducing our measured diameters in 1.1 using the 
ratio of camera focal lengths to determine AD*? for any line from the cadmium 
red value on the plates and applying the phase change correction from the 6, A curve 
in Table I to the resulting fractional order. ‘The wave-lengths then given by the 
three spacers are shown in Table II, in the form observed wave-length—Barrell’s 
wave-length, Barrell’s wave-length being given in column 1. We see that none 
of the residuals exceeds 0-001 A and the mean for any one spacer is less than half 
this value. We may then conclude that our method of wave-length measurement 
even with spacers as narrow as 1-5 mm is capable of giving absolute wave-lengths 
with an accuracy better than o-oo1 A from as few as eight (or even fewer) plates, 








TasB_e II 


Wave-length accuracy using phase change corrections 5 and adopted variations of camera 
focal length. 


(Observed wave-length—Barrell’s wave-length) 


Barrell 1°5 mm 2°5 mm 1omm 

A(A) spacer spacer spacer 
4046°5716 +0:0008 +0000! +0:0008 
4358°3377 + ‘ooo! — ‘0005 + 0008 
5460°7536 — *0003 — ‘0004 + *0004 
5769°5992 — -0008 — +0006 + ‘0002 
5790°6635 +0:0000 —0-0006 — 0-000! 

Mean +0-0000 — 0'0004 +0°0004 
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2. Measurements on vacuum-arc and solar lines 

The measurements in Section 1 were all made on narrow lines and our results 
in Table II show no systematic changes with spacer, i.e. with resolution. The 
2°5mm and 1omm plates were in fact measured by M.G.A. and the 1-5 mm 
plates by S.N. on the same machine. If there were any differences due to 
measuring habit we should expect them to be greatest for the 1-5 mm plates 
and four of these were measured by M.G.A. as well as by S.N., two at 
5460°8 A and two at 4046:6A. ‘The mean difference in the wave-lengths obtained 
from these plates was o-0005 A (algebraic) and 0-0006 A (arithmetic). We can 
therefore conclude that these narrow lines show no systematic effects due to 
personal equation and changes in resolution. But the broader vacuum-arc and 
much broader solar lines are more susceptible to such effects and require separate 
investigation. Arc and solar plates for this purpose were obtained in the green. 
The experimental arrangements were essentially the same as those described 
in Section 1, except that the spectroscope was used in its five prism, diaphragmed, 
form, two 60° prisms, one 30° prism, in a Littrow arrangement. The collimator 
image between the interferometer plates was now rectangular and this image 
was made coincident and central with the circular patch from the cadmium lamp. 

2.1. Comparison of measurements by different observers.—A test for the effect 
of measuring habit on the wave-lengths derived for arc and solar lines was made 
on the three iron lines 5266-6, 5269°5 and 5270:°4A. ‘These are moderately 
strong iron lines in both Sun and arc, but 5270-4 A is blended in the Sun and 
has been measured now in the arc only. Plates for the Sun and arc were obtained 
with the 2-5 mm spacer and five plates for each were measured and reduced by 
both M.G.A. and S.N., S.N. measuring on the original rotating table micro- 
meter and M.G.A. on a new and improved form of this instrument. 


Taste III 
Comparison of measurements for arc lines 


(Differences Ay—A,4) 
Mean A(A) — 526675534 = 5269°5346 = 55270°3538 


Arc Int: 30 60 30 
Plate <A 882 — 00009 — 00006 +0:Coo! 
A 883 + -0000 + ‘0003 + ‘Ooo! 
A 884 + +0003 — ‘0002 + ‘0002 
' A 885 + -o000 + ‘0007 + ‘0006 
A 888 — 00005 — 0°0002 +0000! 
Mean — 0°0002 +0-0000 +0°0002 


The agreement between the two independent measures of the arc plates 
is shown in Table III where the difference between the wave-lengths as measured 
by S.N. (Ay) and the wave-length as measured by M.G.A. (A,) is shown for 
each line and plate. The wave-length at the head of each column is the mean 
value for all the measures corrected to zero pressure in the vacuum-arc; the 
intensities are those listed by Burns and Walters (7). ‘Table IV shows the 
corresponding results for the solar lines, the mean values are corrected for 
Farth’s motion (for intensities we give the Rowland numbers of the lines). We 
see that the agreement for the measurements of the arc plates is very good, the 
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mean values for the five plates differing by only one or two ten-thousandths of 
an angstrom. For the solar plates the differences are greater, but the mean 
measured wave-lengths agree to less than o-oo1 A which is a good deal better 
than might have been expected from Williams’ criterion. 


TABLE IV 


Comparison of measurements for solar lines 


(Differences Ay—A,) 
Mean A (A) 5266°5605 5269°5454 
6 


Rowland No. 8d? 

Plate A 878 +0000! —0°'0021 
A 879 + +0031 — +0006 
A 881 — +0008 + ‘Ooo1 
A 889 + -0008 — 0018 
A 890 +0°0002 +0:0005 

Mean +0:°0007 —0-0008 


2.2. Wave-length determinations using different spacers.—In investigating 
the effect of changing spacers in the interferometer we must be careful to choose 
solar lines as far as possible free from blends. For this reason we have made 
the test primarily on the iron line 5324:193 A Rowland number 7 which shows 
an apparently undistorted profile in the Utrecht Atlas, and Rowlands’ Revision 
shows no line of number greater than —3 likely to affect it. Along with this 
line we have measured the iron lines at 5328-053 A Rowland number 8d? and 
5328-544. Rowland number 4; these are not as free from possible blending 
effects, but they were considered measurable at least with the 1-5 mm and 
2°5mm spacers. With a 5 mm spacer the inter-order scatter for the solar 
spectrum is considerable and must in particular make the measurements of a 
close group of three lines unreliable (there is a chromium line of Rowland number 
2 between the two iron lines at 5328-1 and 5328-5 A). 

Wave-lengths have been obtained for these lines using the 1-5 mm, 2°5 mm 
and 5 mm spacers in the interferometer. The results for the arc lines are shown 
in Table V; columns 3, 4, 5 show for each spacer the difference between the 
mean value for all the spacers in column 2, and the measured wave-length, with 
the standard deviation of the measured wave-length. The measures refer to 
seven or eight plates for each spacer, except for the weakest line 5328-5 A, for 
which the number of observations is six in column 3, three in column 4 and five 
in column 5. The mean values in column 2 are corrected to zero arc pressure. 


TABLE V 
Arc wave-lengths measured with different spacers 


(Mean A— measured A) 
Arc Mean A I°5 mm 2°5 mm 


5mm 

Int. (A) spacer spacer spacer 

30 5324°1777 —0°0010 + 0'0003 +0-°001!I +0°0004 +0°'0000 + 0‘0001 

5° 5328°0372 —0'0018 +0'0003 +0:'0020 + 0'0004 —0'0002 + 0°0001 

15 5328°5313 — 0°0003 + 0°'0003 +0°0007 + 0°0004 — 0'0003 + 0°0005 
Mean —0*0010 +0°'0013 —0°0002 
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Although the differences between the wave-lengths measured with different 
spacers are larger than for the mercury lines, the change with spacer length is 
not of the systematic sort we might fear, for example from a wrong phase change 
correction, and the average differences from the mean values are only of the order 
o-oo1A. ‘The lines were moved relative to the pattern centre, but there was 
not as much variation in line position as for the mercury line plates, and probably 
not enough to eliminate the accidental errors of the experimental set-up. 
The corresponding results for solar lines are shown in Table VI, in this case 
eight plates were measured for each spacer. The mean values in column 2 are 
the means (corrected for Earth’s motion) from the 1-5 mm and 2-5 mm spacers 
only, since, as we have pointed out above, the 5 mm spacer values are likely to 
be uncertain (although we see that the isolated line 5324-2 is not apparently 
very much affected). The average differences of the 1-5 mm and 2:5 mm 
spacer wave-lengths from the mean value are again only of the order of o-oo1 A. 


TaBLe VI 
Solar wave-lengths measured with different spacers 


(Mean A— measured 4) 


Row. Mean A* I°5 mm 2°5 mm 5mm 
No. (A) spacer spacer spacer 
5324°1868 —0:'0008 + 0:0006 +0:'0007 + 0°0004 — 0°0013 +0°0003 
8d? 5328-0463 —0'0016 + 0°0005 +0°0016 +0:0004 — 0'0024 + 0°0005 
4 5328°5376 — 0°0006 + 0:0004 +0:0006 + 00005 — 00022 + 0°0005 
Mean — 00010 +0°0010 —0'0020 


* From the 1-5 mm and 2:5 mm spacers. 


The fact that they are in the same sense and of the same size as those for the arc 
lines supports the suggestion that in both cases they arise from the accidental 
errors of the particular plate orientation and line positions used for these particular 
arc and solar measures. ‘This means that the red shift as determined by different 
spacers is even less subject to error than the absolute wave-lengths. ‘Table VII 
shows the red shifts for the three lines as given by the 1-5 mm and 2-5 mm 
spacers; the 5 mm spacer red shifts are not included because this spacer appears 


‘ 


TaBLe VII 
The red shift given by different spacers 


Row. Mean A (A) 1°5 mm 2°5 mm 
EI. No. Solar spacer spacer 
Fe 7 5324°1868 8-9 mA 9°5 mA 
Fe 8d? 5328-0463 8-9 9°5 
Fe 4 5328-5376 6-6 6°4 


to be already too long for the spectroscopic resolution (14). We see that for 
any one line the difference between the red shifts given by the two spacers is 
well under o-oo1 A. 
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Conclusion 


On the basis of our tests described in Sections 1 and 2 we can conclude first 
of all, that as far as experimental technique is concerned, our procedure by the 
method of circular channels eliminates, or reduces below o-oo1 A, the ‘‘ scale 
error’’ of any region in the visible spectrum. ‘That is the error arising because 
the wave-lengths in the region as a whole are not properly related to the primary 
standard. Moreover a fairly small number of plates suffices to give a wave-length 
accuracy for individual lines of o-oo1 A. In the second place we find that the 
‘personal equation”’ of different observers in deriving wave-lengths from the 
peculiar “‘ circular channel ’’ patterns is very small, both for the narrow vacuum-arc 
lines and also for the broad solar lines. For three arc lines and two solar lines 
any systematic difference in the wave-lengths determined by different observers 
is only 0-0002 A for the arc lines, and even for the solar lines the mean values 
for individual lines do not differ by as much as o-oor A. 

Finally we have some evidence that the micrometer measurement of broad 
lines is more meaningful than has been suggested (12). ‘Thus different observers 
measuring the same plates arrive at very closely the same wave-lengths, and 
also changing the interferometer spacer gives differences, (small and non- 
systematic), which are very closely the same for the broad solar lines as for the 
narrow arc lines, provided the spacer is not too long for the spectroscopic 
resolution. ‘The significance of the micrometer measurements is the most 
difficult and the most fundamental of our problems. There is no doubt that the 
full story of the solar red shifts is to be found in detailed contour measurements 
of high resolution spectra on an absolute wave-length scale. Such a programme 
is, however, impracticable for any but a very few lines, and in any statistical 
approach we are bound to rely on micrometer measures. Our investigations 
here described show that, provided the interferometer spacer is properly related 
to the spectroscopic resolution, red shifts consistent and reliable to under 0-001 A 
can be obtained for unblended solar lines. This is an accuracy probably more 


than sufficient to test the various theories with which we try to account for the 
shifts. 


University Observatory, 
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INTERFEROMETRIC MEASUREMENTS OF WAVE-LENGTHS 
V. THE RADIAL CURRENT INTERPRETATION OF SOLAR RED SHIFTS 
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Summary 


The relativity-radial current theory of solar red shifts was thought 
earlier to imply currents with impossibly high velocities, but with a new 
model atmosphere, based on 0:4 granules of high brightness contrast, 
Schréter finds that the shifts can be explained with streaming velocities of 
only 2-3 km sec~!. There is quite a good agreement between his theoretical 
predictions and centre to limb observations of the shift. Our wave-length 
measures at the solar centre and in the vacuum-arc can also be used to test 
his theory. ‘They enable us to determine the resultant streaming velocity as 
a function of line strength; but these velocities do not agree with Schréter’s 
predictions. We then examine the wave-lengths dependence of the shifts 
using the Allegheny Observatory—Bureau of Standards measures; this test 
also fails to support the radial current interpretation. 





Introduction.—Doppler shifts, due to radial currents in the solar atmosphere, 
have most generally been considered to account for the failure to find the 
exact Relativity shift in the Fraunhofer lines. St John originated this explanation 
in 1928 (13) and since that time the theory has been discussed and modified 
by several writers. McCrea and Mitra (9, 10) showed that the change in 
streaming velocity with depth in the photosphere must be taken into account, 
and it was pointed out in an earlier paper that the velocities deduced must be 
mean values for upgoing and downgoing streams (1). But when both these 
factors were taken into account and when the microscopic circulation was 
associated with granulation, we found that the individual streaming velocities were 
impossibly high. In a recent important paper Schréter (12) develops a more 
elaborate two-stream model of the photosphere using the latest granulation data, 
and he overcomes this objection to the radial current theory. 

Schréter assumes that the observed mean velocities come from very much 
smaller granules than was previously supposed; he considers granules 0”*4 
diameter with high brightness contrast. The contrast is assumed to vary across 
the disk from a brightness ratio, granular to intergranular, of 1-45 at the centre 
rising to a maximum of 1°81 95 per cent of the way across (cos = 0-3), and then 
falling rapidly. In the two-stream model the high contrast is in accordance 
with high temperature differences between the streams; in Schréter’s model 
the differences vary with depth reaching a maximum of goo” abs. at an optical 
depth in the continuous spectrum of t=0°5. The streaming velocities are then 
calculated for any depth by turbulence theory, from the temperature deviations 
from the mean value. The velocities in the effective layers are now only of the 
order 2-3 kmsec~! up and down; not at all improbable values. By detailed 
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calculations of line contours as formed in each stream Schréter shows that the 
composite line becomes asymmetrical, and finds the associated apparent shift 
for lines of different strengths and the change in shift from centre to limb. In 
particular he reproduces fairly closely the centre shifts and centre-to-limb variation 
in the shifts which were discussed in 1948 (1). 

Observations of red shifts across the solar disk give the most obvious test of 
the radial current hypothesis, but our measures of absolute wave-lengths at the 
centre (coupled with the vacuum-arc measurements) give alternative and no 
less significant tests of the theory. With the new measures presented in this 
and the preceding paper we have now determined absolute wave-lengths over 
the spectral range 4200-6600 A, and within that range the line strengths vary 
from practically the weakest to the strongest measurable lines of neutral elements. 
Some description of the new measurements for strong lines is given in Section 1 
of the paper, and these red shifts, together with the earlier results, form the 
observational material upon which Schréter’s hypothesis is examined in 
Section 2.1. In Section 2.2 the radial current interpretation is considered for 
the very extensive measurements from the Allegeheny Observatory—Bureau of 
Standards work. The paper concludes with a very brief further discussion of 
the main features of our measured red shifts, since neither the Oxford nor the 
A.O.—B.S. measures can yet be reconciled with the radial current interpretation. 


1. The observational material 


We have previously measured absolute wave-lengths in both Sun and arc 
for 20 lines, ranging from Rowland numbers 1-8, in spectral regions from 
5080-6600 A. ‘To these lines we have now added those shown in Table I. 
The new lines were chosen to extend our measurements to the strongest 
Fraunhofer lines and, excepting the Balmer lines, we have now included the 
strongest lines of neutral elements in the visible spectrum. On the logarithmic 
scale of intensities which we use, the Balmer lines extend the range only very 
little, and their characteristic rounded contours would make it difficult to 
measure the wave-lengths by micrometer methods. ‘The next strongest, Mgb 
lines and the calcium resonance line, have sharp cores, and with heavy photographic 
densities it is not unduly difficult to measure them. ‘The additional lines in 
Table I are numbered 21-34 (following the first 20) for convenience in referring 
to them in the later parts of the paper. ‘The equivalent widths in the last column 
are in milliangstréms as given by Allen (4,5). ‘There are several sources of 
equivalent widths for the stronger lines, but for the sake of consistency we have 
used Allen’s measures for all 34 lines we shall discuss. 

‘The wave-lengths in both Sun and arc were measured by the method of 
circular channels with details of the arrangements as in Paper IV. For the 
green region the spectroscope was used in its 5-prism form; at 4227 A where the 
spectroscopic dispersion is higher it was used in its 3-prism form. ‘The wave- 
lengths below for lines 21, 22, 23, 30, 31 all refer to a 2:5 mm spacer in the 
interferometer. In the Mgb region both solar and arc observations include also 
some plates taken with a 1-5 mm spacer. The plates were measured and reduced 
as explained in the earlier papers, the ratio of camera focal lengths and the phase 
change correction for the different regions were obtained from the work described 
in Paper IV. Table II shows the absolute wave-lengths now obtained for the 
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TAaBLe I 


Additional lines measured 


Rev. Rowland Table Eq. width 
Line (Allen) 
No. El. Int. Wave-length in A mA 
21 Ca 20d 4226°742 1230 
22 Fe 4 4227°442 170 
23 Fe 6 4233 °613 220 
24 Mg 15 5167-330 752 
25 Fe 5 5167°510 152 
26 Fe 6 5171612 151 
9 Mg 20 5172°700 1200 
28 Ti 2 5173°751 75 
29 Mg = 30 5183-621 1600 
30 Fe 6 5266°565 220 
31 Fe 8d? 5269°552 388 
32 Fe 7 5324°193 254 
33 Fe 8d? 5328°053 341 
34 Fe 4 5328°544 196 


solar and vacuum-arc lines 21-31 with the standard deviation of the mean value 
and the number of observations. For lines 24-29 the values from the two spacers 
used have been given equal weight in forming the mean wave-length. The 
wave-length measurements for lines 30-34 have already been discussed in 
Paper IV. The results given in Table II for lines 30, 31 are those obtained by 
M.G.A. alone; the wave-lengths for lines 32-34, measured by S. Nichols 
alone, have already been given in Paper IV, Tables V and VI. 


Tas_e II 


Measured wave-lengths: solar and vacuum-are 


Sun Arc 

(1) (2) = (3) (4) (5) (6) (7) (8) (9) (10) (11) 

21 Ca 4226-7 "7385 +0°0006 5 -7268 +0°:0003 8 00000 7268 
22 Fe 4227°4 "4369 «= 00006 5 "4270 ©0°0002 5 0°0004 4266 
23 Fe 4233°6 6123 070007. 5 6031 ©0003 § ~=—- 00003 6028 
24 Mg 5167-3 *3308 «=6o°0OII 67 *3236 + 0:0008 6 = 0:0006 *3230 
25 Fe _ 5167°5 "5012 o:0008 7 ‘4879 0°0007 7 0'0003 +4876 
26 Fe 5171°6 ‘6062 00007 12 "5958 070005 7 0:'0003 "5955 
27 Mg 5172°7 6966 00006 12 6859 00006 )=666—Ss( 00006 6853 
28 Ti 5173°8 ‘7490 070008 12 ‘7417 070004 8 o-0001 “7416 
29 Mg 5183°6 6155 o-0cc6 8 ‘6063. 010009 9g 0:'0003 “6060 
30 Fe 5266-6 “5601 ©0009) #5 5541 0°0003 ««5§~—s«0 0006 *5535 
31 Fe 5269-6 *5458 0°0007 5 °5348  0°0003 5 0:'0002 5346 


(1) Line number; (2) El.; (3) Wave-length in angstréms; (4) Mean obs. solar wave- 
length; (5) St. dev. of mean; (6) Number of obs.; (7) Mean arc wave-length; (8) St. dev. 
of mean; (9) Number of obs.; (10) Press: correction in angstréms; (11) Vac.-arc wave- 
length. 
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The solar wave-lengths in Table II, column 4, are a mean value for a slit 
height of 18 mm at the centre of the 180 mm image; they are corrected for the 
Earth’s rotation and orbital motion. 

With the exceptions noted below, the mean vacuum-are wave-lengths in 
column 7 were obtained under our usual arc conditions, pressures not greater 
than 10 cm of mercury, currents less than 8 amps and pole separations not less 
than 10 mm. Self-reversal in the calcium resonance line was most easily avoided 

Tasce III 

Solar red shifts 
(1) (2) (3) (4) (s) (6) (7) (8) (9) 
I Fe 5068°8 383 2°93 z'P°-e’D 9°6 —I'l 25°6 
2 Fe 5074°7 1094 4°20 y°F°-e®G 12°8 +2°0 20°9 
4 Fe 5078907 1092 4°28 y°F°-f°G 8-2 —2°6 17°3 
4 Fe 5079°2 66 2°19 a5P—y5Pp° 9°6 —1'2 20°9 
g Fe 5°79°7 16 0°99 a®F-z5F° 9°8 —1°'0 20°5 
6 Fe 5083°3 16 0°95 a®F-z5F° 10°4 —0'4 21°8 
Fe 6008 -6 982 3°87 z®D°-e® F 8-4 —4°3 14°3 
8 Mn 6013°5 ay 3°06 z§P°-e8S 7°6 —5s'I 12°6 
9 Mn 6016-6 27 3°06 Z®P°—e8S 73 —5'5 15°3 
10 Mn _ 6021°8 27 3°09 Z8P°—-e®S 6°1 —6°7 15°3 
II Fe 6024°1 1178 4°53 y°F°-f°G 9°9 —2°9 17°6 
12 Fe 6027°0 1018 4°06 c®F-v3G° 6-7 —6°1 IO'l 
13 Ca 6493°8 18 2°51 3°D-3d4p*F° II°3 —2°5 22°8 
14 Fe 64.9498 168 2°39 a®H-z>G° 12°9 —o'9 27°1 
15 Ca 6499°6 18 2°SI 3°D-3d4p°F° 10°6 —3°2 14°2 
16 Fe 6569°2 1253 4°71 y’D°-¢g°F 7°0 —6°9 I1'9 
17 Ca 6572°8 ~ 3 0°00 4'S-4°P° 4°6 —9°3 5°0 
18 Fe 6575°0 206 2°58 b*F-z’G° 7°0 —6°9 gl 
19 Fe 6592°9 268 2°92 a®G-y*F° 12°9 —I'l 17°0 
20 Fe 6593°9 168 2°42 a®H-z'G? 15°7 +1°7 I1‘2 
21 Ca 4226°7 2 0°00 4'S-41P° ry +2°7 291°0 
22 Fe 4227°4 689 3°40 Z>Fo—e’F 10°3 +13 40°2 
23 Fe 4233°6 152 2°47 z’D°-e’7D 9°5 +0°5 52°0 
24 Mg 516773 2 2°70 3°P°-4°S 7% 6-32 «6a 
25 Fe 5167°5 CF 1°48 a®F-z'D° 13°6 +26 29°4 
26 Fe 5171°6 36 1°48 a®F-zF° 10°7 -—0%3 29°2 
27 Mg 5172°7 2 2°70 3°P°-4°S 113 +03 232°0 
28 "ER 5173°8 4 0°00 a®F-z3F° 74 —3°6 14°5 
29 Mg 5183°6 2 2°70 3°P°-4°S 9°5 —1°5 308°7 
30 Fe 5266°6 383 2°99 z’P°-e’D 6°6 —4°6 41°8 
31 Fe 5269°6 15 0°86 a®F-—z'D?® II‘2 +0°0 73°6 
32 Fe 5324°2 553 3°20 z2°D)°-e®D 9'2 —2°1 47°7 
33 Fe 5328°1 15 o'91 a5F-z5D° 9°2 —2:'1 64°0 
j 34 Fe 5328°5 a7 1°55 a®F-z2°D° 6°5 —4°8 36°8 


(1) Line number; (2) El.; (3) Wave-length in angstréms; (4) Mult. number; (5) Ex. 
Pot. in volts (low); (6) Transition; (7) Red shift A mA; (8) Resid. shift § mA; (9) 10°W/). 
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by using lower currents and pressures and close pole-pieces, and two of the 
calcium observations and one of the five observations for lines 22 and 23 were 
obtained under these conditions. The close pole-pieces caused no significant 
change in the wave-lengths obtained. ‘Three of the eight calcium observations 
were for the unreversed line and these wave-lengths agreed well with those 
obtained for the sharp absorption core of the self-reversed line. The corrections 
to zero arc pressures were derived, for all the elements, from Babcock’s expressions 
for the depressions of energy levels in the iron atom (6); the corrections are 
all small, as can be seen from column 10. 

Finally in Table II] we have collected together all the red shifts which we 
have so far obtained. ‘The red shifts for lines 1-20 are from Table V of Paper II 
(2); for lines 21-31 they are the differences of columns 4 and 11 in Table II 
above; for lines 32-34 they are the mean red shifts for the 2-5 mm and the 
1*5 mm spacer in ‘l'able VII of Paper IV (3). Iam very grateful to Mr Nichols 
for permission to include his results for these last three lines. ‘Table III shows 
also the physical properties of the lines which we shall need for our subsequent 
discussion. Columns 4, 5, 6 are from Miss Moore’s Revised Multiplet Table (11) 
A in column 7 is the solar red shift in milliangstréms, and 6 in column 8 (the 
residual shift) is the red shift A minus the predicted relativity shift (2-12A x 107°). 
6 is thus the quantity for which we have as yet no agreed physical theory and for 
which we must account. Column g gives a measure of line strength in dimension- 
less units ; it is Allen’s equivalent width for the line W, divided by the wave-length 
A, both measured in the same units; W/A is the unit in which Schroéter expresses 
line strengths. 


2. The radial current hypothesis 


2.1. Observed velocities—On the radial current theory the residual shift 6 
in Table III is a velocity shift and for lines formed at the same geometrical height 
in the solar atmosphere the velocity has a constant value. In general the higher 
the level of formation the smaller the outward velocity we should expect, and 
indeed for very high lines the shift can become red according to Schréter’s 
investigation of resultant line contours. Over medium line strengths it seems 
reasonable to expect the velocity shift to decrease steadily as the effective level 
of formation rises, i.e. with increasing line strength. This is only a generalization 
and ignores questions of abundance, excitation potential and wave-length changes 
in the continuous absorption coefficient, but these should be of secondary 
importance. Fig. 1 shows the observed 6/A as a function of line strength in the 
form log1o®W/A. For convenience in discussion we will indicate by (W/A) 
the quantity 10°W/A. Let us assume that there is a linear relation between 
6/A and log (W/A) which we express by 


5/A=c(log (W/A) — log (W/A) ) (1) 


in which c is constant and (W/A), is the line strength at which the resultant velocity 
is zero. Our quantity 6 can then be expressed in the form 


§=ad+cAlog (W/A) (2) 


where a is —clog(W/A))._ In making a fit to this expression, it seems fairer to 
the hypothesis to omit from the solution the strongest lines 21, 24, 27, 29 (open 
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circles in Fig. 1). ‘They must be formed in the highest layers of the atmosphere 
and possibly do not continue the steady progression of velocity shift and line 
strength which we may expect for the less strong lines. 

From ‘Table III we have then 30 observational equations of the form (2) 
and a least-squares solution of these gives 


5 = — (1764 + 0-392) x 10-®A + (0-gg9g + 0°299) x 10~®A log (W/A) (3) 
in which 6 and A are measured in the same units. Or in form (1) we have 


5/A=0-999 x 107° (log (W/A) — 1-765) (4) 
which is the straight line shown on Fig. 1. Equation (4) gives zero velocity for 
lines where log (W/A) = 1-765, i.e. lines of Rowland number about 8. Our faintest 
line is Rowland number 1, log (W/A)=0-699; for such a line equation (4) gives 
6/A= — 1-07 x 10-® corresponding to an upward velocity of 0-32 kmsec™!, and 
all stronger lines must have smaller velocities than this. 





6/A x 10° T | i ! 
0-8 a 
0-4 
0-0 
i§-4 


-0:8 











| | | | I | | A, i 
0-6 0-8 1-0 1-2 14 1-6 1-8 2:0 22 2:4 26 2-8 
Log W/A x 10° 


Fic. 1.—Residual velocity shift (8/X) as a function of line strength. 





Schroter tested his theory of red shifts against two series of solar observations, 
the Oxford measures of 14 lines at 6100 A of mean equivalent width 0-088 A (1), 
and the Potsdam observations (7) of g lines at 4400 A of mean equivalent width 
0-113 A. The Oxford measures give absolute wave-lengths from centre to limb, 
while the Potsdam observations give only the change in wave-length across the 
disk, and Schréter assumes for them an absolute value of the centre red shift. 
He then shows that both series of observations can be fitted on the radial current 
theory with a resultant outward streaming velocity at the centre of the disk (v) 
of o-42kmsec~!. Fig. 2 shows the Oxford observed A’s, expressed in km sec~, 
as a function of cos 0, 6 being the angle between the line of sight and the normal 
to the solar surface; the range in cos@ and the estimated uncertainty of the 
velocity is shown with each observed point. The theoretical relation is given 
by Schréter’s equation (37) 
Ag (in km sec~') = 0°636—v. B,.cosé 
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in which 0-636 is the Relativity shift in kmsec~! and By, which depends mainly 
on the granulation contrast, is given in his Fig. 8. Curve (a) on Fig. 2 is 
Schroter’s predicted curve with v=o-42 kmsec~!. But our equation (4) above 
gives for the Oxford lines (log (W/A)=1-16) a resultant velocity at the centre 
of the disk of o-184kmsec~, and for the Potsdam lines (log (W/A)=1-41) a 
velocity of o-10o7kmsec~. If we substitute these velocities in Schriter’s 
theoretical expression we obtain the upper curves (b) and (c) respectively in 
Fig. 2, which no longer agree with the observations. The curve (b) would be 
still higher if we had included the four strongest lines in our solution for 
equation (4). 


Red Shift 
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Fic. 2.—Red shift (km sec~1) as a function of disk position. 





Our velocities now deduced are based on absolute measures of both solar 
and arc wave-lengths, which are used to give the general form of the velocity—line 
strength relation. In the 1948 results the wave-lengths were measured by the 
older method of parabolic channels using atmospheric oxygen lines as standards, 
and the vacuum-arc wave-lengths were from a variety of sources in the literature. 
The first five of the 1948 lines have re-appeared in our later measures by the 
method of circular channels (lines 8-12 in Table III, from Paper II), and the 
new centre red shifts are corrected in the direction which the collected observations 
now indicate. ‘The uncertainty in the absolute values of the red shifts in 1948 
means that the observed curve in Fig. 2 may be shifted bodily parallel to the 
velocity axis, but the shape of the curve should not be affected. Such a 
displacement would give limb red shifts well in excess of the Relativity value which 
is shown by the dotted line on Fig. 2. We must note, however, as was pointed 
out in 1948, that the observations in Fig. 2 correspond to low resolution on the 
solar disk and Miss Hart’s results (8) indicate that the rapid increase in red 
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shift must be closer to the limb. The exact form of the centre-limb change in 
red shift is clearly a matter of first importance and we hope soon to have a better 
determination of it, from observations made with the large image given by the 
new solar telescope at Oxford; but these plates have not yet been measured. 

Schréter points out that, while line strength is the most important factor, the 
red shift depends also on excitation potential. It is possible that some of the very 
considerable scatter of the observed points in Fig. 1 may be due to this effect. 
We have therefore examined the residuals of our 30 observational equations for 5 
as a function of excitation potential. ‘The correlation coefficient between the 
residuals and the lower state excitation potentials is, however, quite meaningless, 
being only +0-0080 for the 30 observations. 

2.2. Wave-length dependence of red shifts.—In the last section we assumed 
that our residual shifts were velocity shifts, but we find that the velocities they 
indicate do not account for the observed centre-limb behaviour of the red shifts. 
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Fic. 3.—Line strength-wave-length distribution of the A.O.-B.S. red shifts. 


We will now go further and enquire whether the solar red shifts show the essential 
feature of velocity shifts in their dependence upon wave-length. For a given 
level of formation in the solar atmosphere, say a given line strength, the resultant 
streaming velocity should be constant and the residual shifts 5 should be 
proportional to the line wave-length. Our results in Table III, which were 
designed to cover a large range in line strength, do not give enough observations 
at any one line strength to make this test in its simplest form, but we can use 
again the very large collection of red shifts in the Allegeheny Observatory—Bureau 
of Standards investigations. Fig. 3 shows the line strength-wave-length 
distribution of the 158 observations already discussed in Paper II*. We see that 
if we restrict ourselves to line strengths between 80mA and 125 mA (between 

* The A.O.-B.S. wave-lengths are for the integrated solar disk, and these must be slightly longer 


than centre wave-lengths. Estimating the difference between the two from the observed centre- 
limb change in Fig. 2 and limb darkening observations we find that it is less than o-oo1 A. 
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the dotted lines in Fig. 3) we can include a large number of red shifts with a 
wave-length range from 4100-6800 A, and for this analysis we will take log (W/A) 
as constant; the maximum range in log (W/A) is 1-09-1-47. 

For the 76 lines included within these limits Fig. 4 shows the residual shift 8, 
(observed red shift—Relativity shift), as a function of wave-length. We see that 
the scatter of the points is considerable and the general distribution gives little 
indication of the strictly linear dependence of 5 on A demanded by the radial 
current hypothesis. On this hypothesis the relation between the two is of the 
form 5=mA where m is constant. The mean value of 5/A for the 76 points is 
—o55 x 10~°, and the line (a) on the figure shows the relation between 6 and A 
we should expect from this value. The line is obviously a very inadequate 
representation of the distribution; the scatter of the observed points about it is 
clearly not random, but varies systematically with wave-length. Over the first 
half of the wave-length range there is a marked increase of 5 with A instead of a 
decrease, so much so, that at 5000A the residual shifts are positive, indicating 
an ingoing resultant velocity instead of an outgoing one. In the red the residual 
shifts have systematically larger negative values than the relation indicates. 
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Fic. 4.—Residual shifts as a function of wave-length. (A.O.-B.S. observations). 


The observed 4, A relation obviously requires representation by some form 
more complicated than a linear one. But since a linear relation is predicted by 
the radial current theory we have made a least squares solution for the line of 
closest fit. ‘This is the line (b) on Fig. 4. The equation is 


5 =8-o1 + 2°73 — (2°11 +.0°51) x 10-®A (5) 


in which both 6 and Q are in milliangstréms. If in this expression the second 
term is regarded as corresponding to the resultant velocity demanded by the 
radial current hypothesis, how are we to account for the large and quite well 
determined constant term? On the radial current theory 6 is directly propor- 
tional to A, and the first term in (5) should not be significantly different from 
zero. We must conclude then that radial currents and the Relativity shift are 


not, by themselves, sufficient to account for the observed red shifts and their 
variation with wave-length; some additional hypothesis must be included. 
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Conclusion 


The radial current hypothesis is the only physical theory put forward to 
account for solar red shifts, and Schréter’s very careful investigation on the basis 
of very small, high contrast, granules is made in accordance with all the latest 
theories of the solar atmosphere. Moreover his model atmosphere is able to 
account for other observational phenomena. But when we examine the theory 
for red shifts over a wider range of observational material than Schroter 
considered, both that presented here and the very numerous measures given by 
the Allegheny Observatory—-Bureau of Standards, we find that we are not yet 
justified in accepting his explanation. We lack, at present, the observational 
evidence that the shifts are, in fact, velocity shifts, but if nevertheless we interpret 
them as such, we find that the velocities deduced from the centre of disk values 
do not accord with the change in red shift to the limb. 














Log equivalent “ ergs ”’. 


Fic. 5.—A/A* as a function of log equivalent “‘ ergs”’. 


Observationally our only certain information about the centre shifts is that 
they increase with line strength. In Paper II we found a high statistical correlation 
between the shift expressed in cm~! and the logarithm of the line strength expressed 
in equivalent ‘‘ergs’’. Our range in log equivalent ‘‘ergs’’ is now more than 
doubled and Fig. 5 shows the plot of these two quantities for the 34 points of 
Table III; the factors for obtaining equivalent ‘‘ergs’’ are 1-101 for the new 
lines 21-23 and 0-964—-0-929 for the new lines 24-34. Although the scatter of 
the points in Fig. 5 is considerable, the correlation coefficient 1s + 0-587 and for 
the 34 observations this indicates a significance level better than o-oo1: the 
computed regression lines are shown on the figure. 

We may use our observations, as set out in Fig. 5, to point out some features 
of the shifts. The outstanding and somewhat disappointing feature is the very 
large scatter of the points. ‘The preceding Paper IV shows that we may hope to 
get red shifts accurate to better than o-oo1 A. As far as accidental errors now are 
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concerned the standard deviations in ‘Table II show that the maximum uncertainty 
in A for our strongest and most difficult lines is only 0-o014 A, the vertical lines 
through some of the points in Fig. 5 show the corresponding estimated uncertainty 
in A/A? for these lines, and similarly for some of the weakest lines. It appears then 
that the scatter of the points in Fig. 5 is outside the range of observational 
uncertainty. If it were produced by the peculiar perturbations of atomic energy 
levels we should expect to find members of the same multiplet showing the same 
deviations. We have joined together such lines in Fig. 5 and we see, that among 
the medium to strong lines, there is no evidence for such an effect. Our measures 
include now three resonance lines, 17, 21, 28 in Table III, but these show no 
special peculiarities in their red shifts. It is of course possible that the scatter 
is due to unsuspected blends in the solar spectrum, but it is greater, rather than 
smaller, for the stronger lines and these lines should be much less affected by 
any faint blends. We cannot but conclude that there remains some substantial 
physical effect in the solar atmosphere, influencing the red shifts, which we have 
not yet taken into consideration. So far as the observations go, we have seen 
that it is certainly of such a magnitude as to obscure the presence of shifts due to 
radial currents. 


I gratefully acknowledge the constant help and advice I have received from 
Professor Plaskett; I should also like to thank Miss J. A. Bailey for her valuable 
computational assistance. 


University Observatory, 
Oxford: 


1957 November 21. 


References 


(1) M. G. Adam, M.N., 108, 446, 1948; Comm. Univ. Obs. Oxford, No. 26. 
(2) M. G. Adam, M.N., 115, 367, 1955; Comm. Univ. Obs. Oxford, No. 49. Paper II. 
(3) M. G. Adam and S. Nichols, M.N., 118, 97, 1958. Comm. Univ. Obs. Oxford, 
No. 65. Paper IV. 
(4) C. W. Allen, Mem. Comm. Solar Obs., 5, 1934. 
(5) C. W. Allen, Mem. Comm. Solar Obs., 6, 1938. 
(6) H. D. Babcock, Ap. F., 67, 240, 1928. 
(7) E. F. Freundlich, A. v. Brunn and H. Briick, Z.f. Ap., 1, 43, 1930. 
(8) A. B. Hart, M.N., 114, 17, 1954; Comm. Univ. Obs. Oxford, No. 43. 
(9) W. H. McCrea and K. K. Mitra, Obs., 57, 379, 1934. 
(10) K. K. Mitra, Obs.) 59, 160, 1936. 
(11) C. E. Moore, A Multiplet Table of Astrophysical Interest, rev. ed., Princeton, 1945. 
(12) E. H. Schréter, Z.f. Ap., 41, 141, 1957. 
(13) C. E. St. John, Ap. ¥., 57, 195, 1928. 











all 


ee or le ee ee ee, a, 





XUM 
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Summary 


82 photographs of the globular cluster NGC 121 on the outskirts of the 
Small Magellanic Cloud, covering a total range of 9 years but mostly con- 
centrated in three consecutive seasons, are used to derive periods and light- 
curves of 5 variables within 1’ of the cluster. Of these, three are short- 
period variables (periods between 0-50 and 0-64 days) of Bailey a and 6 type 
with the usual relations between period, amplitude and shape of light-curve. 
The other two, among the reddest and brightest of the cluster-stars, have 
periods of 140-2 and 112:4: days. The light-variations are presented in the 
form of step-estimates on a “ basic ’’’ sequence. 

Four other short-period variables have been found on the plates which 
are too far from the cluster to be members. ‘They are sufficiently faint to be 
regarded as possibly the first instance of RR Lyr variables in the general 
field of the Small Magellanic Cloud. 





NGC 121 (0" 24:4, —71° 48’; 1950) was discovered by John Herschel but 
was first resolved as a cluster on Harvard photographs (1) and recognized as 
a globular cluster by Mohr (2). The Harvard workers accepted the cluster as 
a member of the Small Magellanic Cloud. ‘This view was also adopted by 
Gascoigne and Kron (3) in their observations of colours and magnitudes of 
Cloud clusters; for NGC121 they found P=11-44, P-—V=+0-70. The 
cluster lies 35’ nf the globular cluster 47 Tuc. 

When the faint (19th mag) variables were discovered in the cluster at the 
Radcliffe Observatory (4) it was at first suspected that the cluster might lie 
beyond the Small Cloud. However, Baade’s revision of the distance scale, and 
the discovery of other clusters in the Large Cloud with equally faint variable stars, 
led Thackeray and Wesselink (5) to deduce increased distances of both Clouds 
in good agreement with Baade’s revised scale. At the time of writing, no radial 
velocity of NGC 121 is available, but there seems to be no reason to doubt that 
it is in fact a member of the Small Cloud, especially since the 21 cm radio 
observations have shown that the boundary of the Cloud extends beyond the 
position of the cluster. 

The discovery of these short-period variables was also important in that it 
indicated the presence of Population II objects in both Clouds. In this paper 
periods and light-curves of these variable stars are presented. Special interest 
attaches to such data since the discovery that the RR Lyr variables in the galactic 
nucleus (6) differ from those in many globular clusters as studied by Bailey. 
It is therefore considered worthwhile to publish results for NGC 121 based on 
step-estimates of magnitude, prior to the final establishment of a sufficiently 
faint magnitude sequence. 
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Observations.—A total of 82 blue photographs of NGC 121, all taken with 
the Radcliffe 74-inch reflector, have been used for estimating the variables, 
Of these 29 were taken by Dr A. J. Wesselink, 9 by Dr M. W. Feast and the 
remainder by the writer. Exposures were on Kodak 10320 emulsion, mostly 
of 30 min, but a few were of 10 or 15 min. A few red exposures are also available, 
one being shown in Plate I (6) for comparison with a blue exposure. Owing to 
the great need for series of plates on closely adjacent nights and at large hour 
angles, many plates were taken in conditions far from perfect. ‘The faintness 
of the variables implied that good estimates could only be made on plates of first 
quality. Consequently, some of the original plates could not be used at all for 
estimates and in many cases all that could be said of some variables was that they 
were fainter than some limit brighter than minimum. 

Blink discovery of the variables.—Pairs of plates were examined in the blink 
comparator belonging to the Royal Observatory, Cape, which has very kindly 
been made available on loan to the Radcliffe Observatory by H.M. Astronomer. 
in all a total of nine variables has been discovered (five in the cluster). A few 
others, recorded as possibly variable on blink comparisons, were rejected on 
further examination. Variables 1 to 4 were discovered by the writer at a very 
early stage while concentrating on the cluster itself. V5 was discovered by 
Dr Wesselink while examining the plates for another purpose, and V6 to Vg 
were discovered by the writer while carrying out further blinks of the whole 
area (30’ x 40’) covered by the plates. 


TABLE I 
Positions of variables and comparison stars in field of NGC 121 (seconds of arc) 
Variable AX AY r 
I + 38 es 38 
a “+ §7 + 8 58 
3 _ 2 + 40 40 
4 + 51 — 6 51 
5 — 129: — 288: 316 
6 — 68 — 147: 162 
7 — 510: + 581: 773 
8 ot 2 — 40 40 
9 — 19: — 421: 421 
Comparison 

A — 22 + 108 

B + 138 — 37 

‘a + 43 one 

b + 44 + 3 

c -—. 2 + 53 

d + 67 + 22 

e + 59 — 2 


The blink-search for variables has been carried out in such a way that no 
accurate estimate of the completeness of discovery can be made—as was possible 
in the discovery of the numerous variables in the Sculptor system (7). The 
chance of discovery of each variable in NGC 121 appears to be smaller than in 
the Sculptor system, for which a considerable number of very high quality plates 
were available. Nevertheless it is believed that during the course of blink- 
comparison of 15 pairs of plates the discovery is complete within 4’ of NGC 121 
in so far as the resolution of these plates allow. Further out, where seven pairs 
of plates have been searched and where the chance of discovery is smaller, it 
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cannot be claimed that the four variables found are the only ones above the limit 
of the plates. 

Identification of variables.—The positions of the nine variables relative to the 
cluster were measured on a Toepfer micrometer (also kindly loaned by the 
Royal Observatory, Cape). ‘The measures are listed in Table I, together with 
the distances of each variable from the cluster centre. ‘This table also lists 
positions of seven comparison stars on the same system of coordinates. Variables 
1-4 and 8, together with the comparison stars, are marked on Plate I(a). The 
fields of variables 5, 6, 7 and g are illustrated in Fig. 1. 


' — 

















Fic. 1.—Fields of outlying variables near NGC 121. 
Scale: 1 mm=6"-7. South is at top, following to right. 


Owing to the limited field of the Toepfer micrometer only one CPD star 
could be included in the measures; consequently there is some doubt concerning 
the orientation and scale of the measures in Table I. Nevertheless, the combined 
use of measures and charts should amply suffice for identification of all the 
variables. 

Estimates.—The basic sequence of comparison stars A, B, a to e, was used for 
all the variables nearest to the cluster, viz. Nos. 1-4, 6 and 8. Each variable was 
estimated as so many steps brighter or fainter than a given comparison star, the 
number of estimated steps between each pair of comparison stars being left 
free to vary from plate to plate. At the conclusion of the estimates it was found 





sequence) 


Estimcte (“basic” 





120 


that the average step differences 


de 2, and hence we have the following sequence :— 


Star A B a b c d e 
Step fe) 3 8 12 14 18 20 


If on a particular plate a variable were estimated as c1v2d (instead of c2v2d) 


the resultant estimate was 15 + (i 


Estimates of eight out of the nine variables have been set out in a table deposited 


with the Society. 


Derivation of periods.—¥or a variety of reasons it proved to be far more 
difficult to secure satisfactory long series of plates of NGC 121 on the same or 
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were as follows: AB 3, Ba 5, ab 4, be 2, cd 4, 


nstead of 16). 
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Fic. 2.—Light-curves of variables in NGC 121. 
Ordinates : individual estimates on the “ basic’? sequence. 
Abscissae : phase, with arbitrary zero. Observations covering 0:2 P are repeated right and left 
(a) V1; P-'=0'007132 , 
(b) V8; P-1=0-008g0. ® } 


(c) V2; P-!=1-8783. 
(d) V2; P-!=2-8782 (incorrect) 
(e) V3; P-1=1°56284. 
P-!=1-96651. 


(f) V4; 
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consecutive nights than of the Sculptor system. One reason was that NGC 121 
never attains a zenith distance less than 46° at Pretoria. Consequently the 
derivation of periods for these variables presented considerable difficulty in some 
cases. ‘lhe period of V4, which has a very sharp rise to maximum, could be found 
easily, especially as the series of plates included a fortuitously large number of 
maxima. 

‘The individual observations of the five variables considered as cluster-members 
are plotted against phase in Fig. 2, according to the periods found. V2 is plotted 
twice, for P~! = 1-8783 considered to be the true value, and also for P~! = 2-8782, 
which is considered incorrect. For the two long-period variables, closely adjacent 
observations are grouped together. 


"TABLE IT 


Mean light-estimates of NGC 121 variables. 
(Estimates on basic sequence) 


Vi V2 V3 V4 V8 
Ph Est n Ph Est n Ph Est fv Ph Est n Ph Est n 
0°65 14°2 4 o'80 20'2 9g o-77. 17°F «10 oO°1t 20°77 o-os G60 3 
0°75 15°2 4 098 19°83 10 O'O!I 19°5 13 0°29 20°2 7 o'16 8:5 3 
0°83 9°5 2 O°13 16°5 2 O'I5 15°5 6 0°44 16°90 5 O°92. L2°§t 1 
0°95 5°5 4 o°T7 t0°s: 2 0°22 12°5 6 0°50 10°2 12 0'26 8:0 3 
0'06 o'O 4 0°26 13°72 9 0°34 15°5 9 0°56 11°96 o-s7 G7 4 
org —ro 1 0°34 16:0 4 0°44 16°2 5§ 0°66 15°77 5 O°s4 42°55 7 
0°25 oOo 2 O45 17°5 «4 oss 177 6 o-97G 1774 «4 o-7T 46 5 
0°35 23. 3 0°53 181r 5 0°65 188 5 0°90 I19°0: 3 O'Or 5°! § 
0°44 5°8 2 0°66 19°5 7 
0°55 9°0 3 
Vi V2 V3 V4 V8 
P 140°2 0°53240 0'63986 0°50852 112°4: 
JD 24300004 4219. 277°36 4277°59 4331°28 3619* 
Mmax —1°O 10°3 12°5: 10°2 4°5 
Mmnin +1§°2 20°2 19°5 20°7 115: 
m 7% ¥5°2 16:0: 15°4 8:0: 
(m—M)/P 0:28 o'10 0°23 o'10 0°4: 


* Epoch of minimum for V8. 


Mean estimates of these five variables at different phases are tabulated in 
‘Table II, together with the mean estimate at maximum and minimum. It follows 
that the median estimate of the three short-period variables is 


m (est.) = 15°44 
in the “basic” sequence, i.e. almost precisely midway between comparison 
stars c and d. 
For the two long-period variables we have, with double weight to V1, 


m (est.) = 7°4. 


Notes on individual variables: cluster members. 

Vi. Period 1404-2, light-curve Fig. 2(a). ‘This star appears to be a regular 
long-period variable and is considerably redder than any other cluster member 
with the exception of V8. Its minimum photographic light is almost precisely 
equal to the median magnitude of the short-period variables. In the diagram 
of the light-curve, different cycles are denoted by different symbols and the 
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small scatter may be real. The rise to maximum is steeper than the decline and 
contains a suggestion of a still-stand. In considering the shape of this and the 
other light-curves it must be borne in mind that if the step-differences vary in 
magnitude along the sequence, as may well be the case, then distortions of the 
curves will occur. 

V2. Period o4-53240, Fig. 2(c). ‘The derivation of this period proved to 
be more difficult than for any of the other cluster members. However, it is believed 
that the reciprocal period P-!=1-8783 as used in Fig. 2(c) is correct. The 
observations are replotted in Fig. 2 (d) for a reciprocal period P~! = 2-8782 and the 
scatter is appreciably greater (the value 2:8783 gives an equally bad plot). The 
dashed curve shows the decline from maximum as recorded on a long series on 
one night (J D 2434277). 

V3. Period 04-63986, Fig. 2(e). This variable, with period longer than V2 
and V4, has aslower rise to maximum, and apart from two well-established bright 
maxima has a considerably smaller amplitude. Thus the star is of Bailey type 4, 
as compared with V2 and V4 of type a. 

V4. Period o4-50852, Fig. 2(f). A typical a variable with steep rise to 
maximum and relatively large amplitude. 

V8. Period 1129-4?, Fig. 2(b).. A very red star with an apparently regular 
period, but there is more evidence of fluctuations about the mean curve than 
for V1. The light-curve is rather peculiar in that more time is spent near 
maximum than minimum. ‘There is some resemblance to the rounded light-curve 
of TT Cen(P=g0"%-7). Itis considered more likely to belong to the red long-period 
variable class than to be an eclipsing variable, although Dr Wesselink has 
reminded me of the similar light-curve of the eclipsing variable W Cru (P= 1984:5). 

A possibility remains that there is a double cycle with period 225 days, with 
uneven maxima and minima. Unfortunately one set of minima is very poorly 
represented on this period. If there is a double cycle in a period of 225 days, 
the star bears some resemblance to the RV ‘Tau class, but no galactic RV Tau 
star is known with a period greater than go days. 


Non-members of NGC 121 


Vs. Although no series of observations on the same night covers the full 
rise from minimum, the initial rise is shown on JD 2433977 and a rapid decline 
is shown on several occasions. ‘The observations are well represented by 
P-1!=1-89037; this indicates a Bailey a-type variable with m— M=o:12P. 
Although a “‘local”’ sequence was used for the estimates, the median magnitude 
is found to be about the same as that for V2 to V4. However, the star is too far 
from NGC 121 to be a member. 

V6. ‘The nearest to NGC 121 of the variables considered as a non-member. 
At an angular distance of 162” from the centre of the cluster, a minimum linear 
distance of 40 pc is implied, and is it very doubtful that a globular cluster would 
extend this far. According to Kron (8), NGC 121, which appears as No. 2 on 
his list, has a radius of 10 pe. 

Only three maxima were observed, JD 2433654 (two plates), 2434217 and 
2434304. On the latter occasion a rapid decline occurs, while three days later 
the star appears to have been caught on the rise. No period has been found to 
fit all the observations satisfactorily, but a reciprocal period P 1=1-660 is 
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considered possible if one allows considerable scatter in the brightness at 
maximum. 


V7. ‘The distance from the centre of the plates is so great that the images 
are seriously affected by coma and the estimates are extremely difficult. The 
behaviour on JD 2434277 marks it as a short-period variable but no derivation 
of period has been attempted. 

Vg. This variable is subject to rapid fluctuations but no certain period has 
been found. As for V7, the individual estimates were rather difficult. ‘The 
amplitude is probably smaller than for the other variables and may vary with 
time. ‘There is no indication of a steep rise to maximum. The behaviour suggests 
that of a c-type variable and a reciprocal period P-!=2-6209 appears to be the 
most probable value. However, P~!= 1-623 also fits the observations moderately 
well. ‘The mean magnitude is about the same size as that of the other short-period 
variables or if anything slightly brighter. 

This star is the nearest of the nine variables to the globular cluster 47 Tuc 
whose centre lies 26’ away. It is not considered likely that this or any of the other 
variables is connected with 47 Tuc. 

Conclusions.—The five variables 1, 2, 3, 4 and 8 are all within 1’ of the centre 
of NGC 121 and are to be considered as certain cluster-members. ‘The three short- 
period variables (2, 3 and 4) follow the familiar relationship between period, 
amplitude and shape of light-curve as found in globular clusters in the Galaxy, 
notably as found by Martin (g) in w Cen. Their properties are not those of the 
short-period variables found by Gaposchkin (6) in the galactic nucleus. 

As compared with a density of one variable per 0-63 square minutes in the 
area centred on NGC 121, we have found only one variable per 300 square minutes 
in the surrounding neighbourhood. Allowing for incompleteness of discovery, 
this latter density is certainly under-estimated, possibly by a factor two. V6 is 
the nearest of these outliers to NGC 121 but it is too far to be considered a member. 
We therefore consider V5, 6, 7 and g to be all field variables. Since all four of 
these variables have short periods and have median magnitudes approximately 
equal to those within NGC 121, the interesting possibility emerges that they 
represent the first instances of RR Lyr-type variables in the general field of the 
Small Magellanic Cloud. The only way to test this hypothesis is to examine 
other fields in the region of the Cloud. 

By contrast with the above minimum figure of one variable per 300 square 
minutes, Shapley (10) in his study of RR Lyr variables within the Galaxy found 
an average density of one variable per 6 square degrees (with a magnitude limit 
of about 16) the greatest concentration being in MWF 209 where he found 
one per 2°5 square degrees. Baade (11) in a search for faint RR Lyr variables 
in high galactic latitude found five of about 17 mag. in 8-5 square degrees. ‘These 
concentrations are so much less than that found in the neighbourhood of NGC 121 
that even allowing for the fact that the magnitude limit of the present search is 
some two magnitudes fainter than Baade’s, we find considerable support to the 
hypothesis that variables 5, 6, 7 and g of this paper are actually outlying members 
of the Small Cloud. 

The similarity of NGC121 to globular clusters within the Galaxy in the 
behaviour of its variable stars proves the existence of a Population II component 
to the Small Cloud, and the existence also of a general field of type II now appears 
to be indicated. 
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I am greatly indebted to Drs Wesselink and Feast for their valuable assistance 
in securing additional photographs of NGC 121, so essential to the derivation of 
periods, and further for reading and commenting on this paper. 


Radcliffe Observatory, 
Pretoria. 
1957 August 31. 
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Summary 


Radial velocities of 314 early B type stars with calculated distances 
(based on photoelectric colours and MK classifications) greater than 1 kpc 
are analysed for differential galactic rotation. The incorporation of Radcliffe 
velocities of southern stars renders the distribution in longitude much more 
uniform than in earlier analyses. At distances greater than 2 kpc there is even 
a marked deficiency of northern material. 

First-order analyses of stellar and interstellar velocities are carried out 
for various distance groups (Section 5, Table III). A negative K-term 
(Section 6) is found in the region 1 to 2 kpc, but more distance groups do not 
suggest that this is due to a general contraction of the galaxy. ‘The possible 
connection with the deviation of the vertex is briefly discussed. ‘The K-term 
in the region 1 to 2 kpc is diminished considerably when a systematic correc- 
tion, apparently applicable to the Lick velocities of faint B stars (Section 3), 
is taken into account. The effect of this correction is also to increase the 
derived values of the Oort A constant. The group of stars at distances 
greater than 3-0 kpc yields an abnormally low value of A, but this is regarded 
as chiefly a result of selection leading to a large proportion of overestimated 
distances in this group. 

The velocity dispersion found for the inner part of the galaxy is not 
found to have an abnormal value (Section 7). 

The Camm linear approximation is used (Section 8) to derive a value of the 
Oort A constant. The final value A =17-5 +1°5 km/sec/kpc is adopted. 

Two partly independent methods are used to derive a value of R, the 
distance to the galactic centre; the most probable value is 8-9 kpc (Section 10, 
Table IX). A value of 0-5 for the ratio of the axes of the velocity ellipsoid 
appears to be appropriate to distant B stars. 

The calibration of the MK luminosity scale is tested (Section 11) for 
stars in the range Bs to Og5. 

A curve of differential rotation w(R) is derived (Section 12). The results 
show no significant deviation from the radio results and suggest that w” is posi- 
tive in the range of R covered by the observations. 





1. Introduction.—In determining the constants of galactic rotation according 
to the Oort-Lindblad theory the radial velocities of the O and B type stars have 
always played an important role. A most complete and authoritative discussion 
of the early observations and analyses was given by Plaskett and Pearce (1) in 
1936. Since that time several other analyses have been undertaken, some 
depending partly on new observations, others modifying the method of analysis 
in various ways. The principal deficiency in all this work has been the lack of 
velocities for faint B type stars in the southern Milky Way. ‘This deficiency has 
now been remedied by observations at the Radcliffe Observatory (2, 3) and 
in fact the position has been reached that for the more distant stars there is a 
greater concentration available for analysis in the southern than in the northern 
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Milky Way. In view of this substantial increase in observational material it 
is obviously desirable to carry out a new analysis. This is particularly appropriate 
at the present moment since some recent investigators have deduced rather 
widely discrepant values of Oort’s rotational constant A. 

Besides the new southern radial velocities we have also been able to take 
advantage in the present work of two other significant improvements in the 
observational data referring to B type stars. These are the recent determination 
of photoelectric magnitudes and colours, and of spectral types on the MK system 
for large numbers of B type stars in both hemispheres. 

2. Data selected for analysis.—In the present work we have attempted to 
restrict the observational material to those stars with reliable magnitudes, colours, 
types and radial velocities. ‘To this end a list was first compiled of all stars in 
the General Catalogue of Radial Velocities (GCRV) (4) with spectral types Bs 
or earlier and from this list stars were selected which satisfied each of the following 
three criteria: 

(1) The quality of the velocity is a, b or c. 

(2) The absolute magnitude of the star is known on the basis of the 
calibration of the MK system (e.g. Of stars were omitted as their 
absolute magnitudes are rather uncertain). 

(3) The photoelectric magnitude and colour of the star are known on the 
B,V system or can be converted to that system (U magnitudes are 
also known for many of the stars discussed but we have made no use of 
them in the present work). 

Criteria (2) and (3) effectively reduce the list to stars occurring in the catalogues 

of Morgan, Whitford and Code (5) and Hiltner (6). 

To this list were added stars from the recent Radcliffe Catalogues (2, 3) 
which satisfy approximately the same conditions. The following points con- 
cerning these stars are to be noted: 

(a) If the velocity was considered constant it is included; otherwise the star 
was rejected. It is possible that this is a somewhat more stringent condition 
than (1) above. 

(6) For stars taken from Table II of (2) and (3) the MK type given there is 
adopted. There is some overlap between the Radcliffe MK types and those of 
other workers (principally (5) and (6)). Since the overlap is only fragmentary 
it has seemed best to retain the Radcliffe types throughout for these stars. 

(c) The magnitudes and colours adopted for stars from (3) (Table II) are 
mainly due to Oosttérhoff (7). A few values were, however, taken from other 
sources (e.g. (5§)). The magnitudes and colours adopted for stars in (2) 
(Table II) are mainly unpublished measurements at the Cape Observatory, 
kindly supplied by Dr R. H. Stoy. Since photoelectric magnitudes and colours 
were not available for about 40 stars in Table II (2) these stars had to be rejected. 

(d) For some of the stars in the first list (that selected from the GCRV) new 
velocities have been published in the Tables I of (2) and (3). Generally we have 
adopted the mean of the Radcliffe and GCRV velocities unless one or the other 
was known for some reason to be uncertain. It should be noted that for these 
stars the Hiltner (6) or Morgan, Whitford and Code (5) spectral types were 
adopted for the sake of consistency in preference to the Radcliffe types. 

For the stars in the list prepared as above, it is possible to determine for each 
star its absolute magnitude, the correction for interstellar absorption, and the 
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| it distance. In the case of those stars which are in Hiltner’s catalogue (6), Hiltner 
ate has himself given his deduced distance moduli and these have been adopted in 
her our work. ‘The intrinsic colours used by Hiltner are given in his Table 6 and 
are essentially the same as those given by Morgan, Harris and Johnson (8). 
ake The absolute magnitude calibration of the MK system used by Hiltner is given 
the by Johnson and Hiltner (g). For stars not in Hiltner’s catalogue the magnitudes 
ion and colours of Morgan, Whitford and Code (5), Oosterhoff (7) and Stoy 
em (unpublished) were converted to the B,V system by the following (usually 
accepted) formulae: 
an (a) Oosterhoff’s C,; B-—V=1-107C,+1°371, 
rs, (6) Morgan, Whitford and Code’s C,; B-—V=2-07C,+0°30, 
r™ (c) Stoy’s SCI and SPg B—V=0-85 SCI + 0-26, V = Spg —0-g2 SCI — 0°06. 
Bs The same set of intrinsic colours used by Hiltner was then used to deduce the 
ing interstellar reddening and absorption. A value of 3-0 for the ratio of total to 
selective absorption was used. In deducing the absolute magnitudes of the stars 
from the MK types the calibration of Morgan and Keenan (10) was adopted. 
the This calibration is slightly different from the one adopted by Hiltner. The 
vale differences have been discussed by Johnson and Hiltner (g) and from this 
discussion it might perhaps be concluded that the Morgan—Keenan calibration 
the is the more suitable of the two for the group of stars studied here and for this 
“— reason we have adopted it for stars not found in Hiltner’s catalogue. In any 
a case the differences between the two calibrations are small, probably considerably 
smaller than the uncertainties of the quantities involved; the fact that a slightly 
— different calibration has been used for these stars as compared with those taken 
from Hiltner’s Catalogue will not affect the analysis in any significant way. It 
3) may be noted here that for stars classified as luminosity class I at Pretoria we 
ag have used a mean of the Ia and Ib luminosities. 

Finally from the catalogue prepared in the above manner, those stars with 
on heliocentric distances (r) equal to or greater than 1 kpc were selected for 
-_ analysis. From the point of view of studies of galactic rotation the motions of 

the more distant stars are of considerably greater interest, at the present time, 
wis than ones near the Sun. Furthermore, whilst a great deal of the previous work 
al has been restricted to regions nearer the Sun than 1 kpc we now have sufficient 
ary data to investigate more distant regions. It may also be noted that in restricting 

ourselves to stars beyond 1 kpc we avoid the difficulties associated with the 
one Sco—Cen cluster. 
= With all these restrictions we are left with a catalogue of 319 stars. Of these, 
(2) four stars were rejected because of abnormally high galactic latitude (b> 25°) 
ry, and one (HD 151397) was rejected because of its very high velocity (from four 
uss plates) in the opposite sense to galactic rotation theory. The interstellar velocity 
ill from this star was however retained for analysis. We are thus left with 314 stars 
_— of which 268 have measured radial velocities of the interstellar Cat lines. Of 
ave | these 314 stars 122 depend on Radclifie velocities alone. The distribution of the 
eer | stars is shown in Fig. 1, in which the stars are shown projected on the galactic 
on plane. It should be pointed out that although the original selection was for 
ere spectral types earlier than B5 in the GCRV (and a similar restriction on HD types 

for the southern stars) a few stars in the final list have MK types later than Bs. 
- No attempt has been made to use the present material to deduce the local 
he } solar motion. We have adopted the value of +20 km/sec towards R.A. 18", 
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Dec. +34°, and have corrected all the velocities for this motion using the charts 
prepared by Smart (11). 


(( 


Oc 





. 


Fic. 1.—The distribution in the galactic plane of the B type stars discussed in the present paper. 
The position of the Sun is indicated by the large cross at the centre of the figure. The complete 
circles are drawn at 1 kpc intervals. The large filled circle indicates the position of the galactic 
centre if this is at Ry=8-2 kpc 1)>=328°. The region used in the Camm analysis (Section 8), 
| R—R,| < 1kpce, is enclosed between the two circular arcs. The 9 stars marked by small crosses 
are those used to determine Ry by the second method of Section 10. Some points represent more 
than one star. : 

3. Comparison of Radcliffe velocities with other Observatories.—The inter- 
comparison of velocities determined at various Observatories which appears 
in Table I of the introduction to the GCRV suggests that no serious systematic 
errors exist for the B stars, and accordingly no systematic corrections were there 
applied. Examination of the Catalogue shows that, while there is a greater degree 
of observational overlap for faint B stars than for faint stars of other types, the 
intercomparisons of velocities rest mainly on stars brighter than 7™-o. Yet it 
is the material on stars fainter than this limit which is of special interest in modern 
analyses of galactic rotation; and it is for these fainter stars that the use of 
different spectrographs and cameras (as mentioned in the GCRV) may invalidate 
the general derivation of systematic differences between various Observatories. 
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It is of considerable importance to the present analysis to investigate whether 
appreciable systematic errors exist for the fainter B stars. 

Observations of standard stars with the Radcliffe Cassegrain spectrograph 
have shown a satisfactory agreement with the GCRV system (2, 12, 13), 
although there is some indication that with the c camera (49 A/mm) the Radcliffe 
velocities are I to 2 km/sec too positive compared with the GCRV. ‘The B stars 
present a special problem in that Petrie (14) has shown that at 49 A/mm the 
important lines He 4471 and 4026 require large corrections to the wave-lengths 
usually employed, which vary with luminosity or line-width. The necessity 
for such corrections has been confirmed at the Radcliffe Observatory when 
coarse-grained emulsions are used. 

Concurrently with the Radcliffe observations of southern O and B stars 
which are used in this paper, we have observed 97 B stars previously observed 
elsewhere. These include 23 stars suggested by Petrie for a triple overlap 
between Victoria, Lick (Neubauer) (15) and Radcliffe. 

The accompanying Table I shows the comparison of Radcliffe velocities of 
B stars with the other main contributors—Victoria, Mt Wilson and Lick. 
Victoria A refers to published velocities from the GCRV, Victoria B to a selection 
of stars for which revised velocities (based on the new Victoria system of wave- 
lengths) have been kindly supplied by Dr Petrie. Lick A refers to Lick velocities 
from the GCRV for stars brighter than 6™-5; they are mainly the results of the 
Lick Chile expedition. Lick B refers to stars fainter than 6™-5 ; these are mainly 
determined by Neubauer ((15) or only published in GCRV) and two stars 
brighter than 6™-5 taken from (15) are included in this group rather than Lick A. 


TABLE [| 


Comparison of Radcliffe B star velocities with those of other observatories 


No. of Residual No. of Residual Mean of 

Observatory Stars - (stars) and Stars with Ca 11 and Stars and 
s.e. Ca II s.e. Ca 11, and 

s.e. 

Victoria A 23 +2°0 +0°9 10 +r +1°7 +1°8 +0°8 
Victoria B 12 +o73 +14 +073 +174 
Mt Wilson 22 +3°7 +12 17 +06 +2°0 +24 +1°0 
Lick A $2 —0o4 +13 —o4 +13 
Lick B 33 +114 +1°9 24 +6°3 +2°8 +9°2 +1°6 


Residuals are in the sense Radcliffe minus other observatory. 


In all cases except Lick B the agreement between Radcliffe velocities and 
other determinations is satisfactory. This is particularly so for the new Victoria 
results. The difference between Mt Wilson and Radcliffe could be fully accounted 
for if the Victoria wave-lengths for He 4471 and 4026 were applicable to the 
Mt Wilson material, but whether this is likely to be so cannot be stated from 
available information. The systematic difference from the fainter Lick velocities 
is large and is moreover important for the discussions of the present paper. It is 
not possible at the present time to explain these differences satisfactorily, but in 
view of the good agreement of the Radcliffe results with the other Observatories 
and also with the Lick observations of brighter stars, we feel justified in ascribing 


9 
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the result to a systematic error affecting the Lick B velocities. The final solutions 
in the present paper are based on corrections of +10 km/sec applied to the 
Lick B velocities; but many solutions are also carried out for comparison without 
applying any correction. It is clearly desirable that further observations of the 
faint Lick stars should be made, preferably by a number of different Observatories. 
4. The methods of analysis for galactic rotation.—It will be convenient to collect 
together in this section some of the equations (mostly well known) on galactic 
rotation which will be discussed in later sections. ‘The following notations are 
used : 
V = radial velocity corrected for solar motion. 
l, b= galactic longitude and latitude (Ohlsson’s system). 
1, =longitude of galactic centre. 
A=I1-|,. 
r = distance of star from Sun. 
R=distance of star from galactic centre. 
R, = distance of Sun from galactic centre. 


p=r/Rp. 
w, Wy=angular velocity about galactic centre for any star and at R,, 
respectively. 


eR, €yR, = velocities of expansion of galaxy at distance R and R, respectively. 

k= 4/a/b=ratio of axes of the velocity ellipsoid. 

Since the system of B stars is highly concentrated to the plane of the galaxy 
we shall consider only motions in that plane and also neglect variations of the 
laws of force with distance from the plane. Then in general the state of motion 
in the galaxy may be represented by a rotation and an expansion or contraction 
of the system. It may be shown (e.g. Trumpler and Weaver (16, p. 581)) that, 
neglecting a star’s peculiar motion, its. radial velocity (V) corrected for solar 
motion may be represented by the equation 

V = R,(w — wo) sinAcos b— Ry(e — €9) cosA cos b + er cos? b. (1) 

In practice it is not possible to employ (1) with the available material to 
solve for (w—wy) and (e€—€,) as functions of R and in all work of this kind some 
approximation to (1) has to be employed. It may be mentioned that the equation 
deduced by Edmondson (17) for spiral motion in the galaxy is a special case 
of equation (1). 

If we assume that there is no general expansion or contraction of the galaxy 
then (1) reduces to 

' V = Ry(w— wy) sinAcosb. (2) 
This is Bottlinger’s well-known equation which has been used extensively in 
recent years, especially in analysing the results of the 21 cm radio observations. 
Its application to the present data will be discussed in a later section. The 
assumption that there is no general expansion or contraction of the galaxy (zero 
** K-term’’) is a basic one in most of the 21 cm work and it is very desirable to 
check it as far as possible with the optical data. For this reason it is important 
to discuss approximations of (1) which do not assume zero K-term. 


If we restrict ourselves to regions where |R— R,| is small then we may write 
approximately 


w= w+ (R—Ry)a’, (3) 
e=eo+(R—Ry)ey. (4) 
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Dashes indicating differentials with respect to R. 
If we further follow the normal practice and put 
A= — 4a,’ Ry) =Oort’s constant (5) 
and 
C= — }€, Ry (6) 
we find 
V =e9r cos? b —2A(R— Ry) sinAcosb—2C(R— Ry)[(r/Ry) cos?b—cosAcosb]. (7) 
This is Camm’s linear approximation (18, 19) and may be employed in practice 
to yield values of A, €¢) and C. Now, and R are related by the exact equation 
R?= R,? +17? cos? b—2Ror cosAcos b. (8) 
If, in addition to the approximations (3) and (4) above, we assume that r/R, 
is small, then we may write approximately 


R—R,)=-—rcosAcosb. (9) 
The general form (1) then further reduces to 
V = Apr cos* b + Ayrsin 2(A—h) cos*b (10) 
where 
Ay = €o + €0'Ro, (11) 
A,= b Rov (wo)? + (eo (12) 
and tan h=€,'/w9’. (13) 


Equation (10) may perhaps be called the generalized Oort equation. The well- 
known Oort equation assuming zero K-term follows immediately from (10). 
It is 
V = Arsin 2A cos? b. (14) 
It has frequently been assumed that the K-term arises in some other way 
than from general expansion or contraction of the galaxy and it is then introduced 
as a constant term in equations (1), (2), (7), (10) and (14). For instance (14) 
is then written 
V=K+ Arsin 2A cos*b. (15) 


In an attempt to derive a value of R, from the present material (Section 10) 
we have found it convenient to use a formula due to Titus (20) : 


AR I 
V(A)= 4 — 
) ee t | sina (16) 
Equation (16) is based on the assumption of circular orbits in the galactic plane 


and general dynamic theory (21) which gives a motion at distance R from the 
galactic centre of the form 





V(R)=pR/ (a0 + b2oR*). (17) 
This approach to the problem may be criticized on the grounds that 
equation (17) postulates a certain force law. Another possible approach has 
been suggested by Weaver (22) in which the determination of R, is likely to 
rest on a highly restricted number of radial velocities. In this method, the 
Camm linear approximation to the Bottlinger formula is used without the terms 
involving non-circular motion. Then, with |R—R,| small, we may write 
} V=R,(R—R,)wo’ sinAcosb= —2A(R—R,)sinAcosb. (18) 
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We shall now proceed to apply some of these equations to the observational 
material. 

5. Generalized Oort analysis.—The Oort approximation, that is a formula of 
the form of (10), (14) or (15), depending on how the K-term is treated, has 
frequently been employed in investigations of galactic rotation. It is strictly 
applicable only to small values of r. Weaver (23, 24, 25) has claimed that very 
different results may be obtained from this approximation than from equation (2) 
and he has attributed this to ‘‘mathematical bias’’ in the use of the Oort 
approximation. It is therefore of considerable interest to analyse the present 
data according to (10) for comparison with the other methods of analysis dealt 
with later. Whilst the data will be analysed in this section according to (10) 
it will become clear that the analysis also applies as it stands if equation (15) should 
appear more appropriate. 

The stars were divided into three groups (A, B and C) according to distance 
as is shown in Table II, where the distribution amongst the various MK luminosity 


TaB_e II 


Statistics of stars used in Oort-type solutions 


Group A B Cc AB 
Distance range (kpc) ro<r<2'0 2°0<r<3:0 r>3°0 1°§ <7 <2°5 
Mean distance (kpc) 1°44 2°34 4°31 1°94 
Adopted 7 (kpc) rs 2°43 4:0 2°0 





Number of 


Stars 172 95 47 197 
Ca 11 velocities 141 86 4I — 
Stars luminosity class I KY 43 31 49 
Stars luminosity class II & III 48 26 15 35 
Stars luminosity class IV & V 67 #7 ° 39 
No luminosity class (early O type) 20 9 I 14 


classes is also shown. ‘The group AB for which a solution has also been carried 
out is of course not independent of groups A and B. 

Normal points were formed in each of the seven groups (four groups of stars 
and three of interstellar gas) for each 10° of longitude. Unit weight was given 
to all GCRV velocities and to all Radcliffe velocities except those catalogued 
with a colon, to which half weight was given. In solving equation (10) for a 
particular group it is very desirable to avoid distortion due to stars in one longitude 
being at a different mean distance from stars in another longitude. The normal 
points were therefore corrected linearly (i.e. assume (10) to apply) to distances 
I°5, 2°3, 4:0 and 2-0 kpc for the groups A, B, C and AB respectively. The Call 
velocities were also corrected in the same ratio. In groups A, Band AB the changes 
resulting from this correction are negligibly small. In group C the changes 
involved are somewhat larger due to the greater range of distances and the smaller 
number of stars in this group. The normal points for the seven groups are 
plotted in Figs. 2 to 8. In view of the discussion given earlier of the Lick 
velocities of faint B stars it has been felt necessary to investigate the effect of the 
proposed corrections on the galactic rotation constants. Where two points 
occur, in Figs. 2 to 8, at the same longitude the more positive refers to the 
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al ' 
corrected results and the other to the uncorrected results. In view of this 
of duplication we have therefore 14 groups to investigate. Equation (10) was solved 
as by the normal method of least squares for each of the 14 groups. The results 
ly and their standard (r.m.s.) errors are given in Table III whilst the curves in 
Figs. 2 to 8 are drawn according to these solutions. It should be noted that in 
ry g g 
2) Taste III 
rt Least square solutions of Oort-type analysis 
nt 
alt Group r Ay Aor Igth A Ay 
: ) Corrected solutions 
A (stars) I's +27°'8+1°8 —4'2+13 325°+2° +18-6+1°2 —2°8+-0'9 
A (Ca 11) (0°75) +14'9+1°2 —2'2+0°9 328°+3° +18-6+1°6 —2°8+1°2 
ce AB (stars) 2°0 +35°'2+2°2 —1'4+1°6 328°+2° +17641'1 —o7+08 
ty B (stars) a3 +37°6+4:2 +66+2°8 320°+3° +16°341°8 +2912 
B (Ca 11) (1°15) +204+2°7 +3:0o+2'1 319°+4° +17°7+253 +2:6+1'8 
C (stars) 4°0 +35°747°0 —o6+5'8 326°+7° + 8-9+1°8 —o2+1°4 
C (Ca 11) (2°0) +-27°-742°2 +oo+1°8 330°°+3° +138+11 +o'4+0°9 
Uncorrected solutions 
A (stars) 1°5 +2qit21 —6:8+1°5 326°+3° +161+1°4 —4°5+1°0 
A (Ca 11) (0°75) +13-2+1'2 43407 332°+3° +16°5+1°6 —5'7£0°9 
+5 AB (stars) 20 +32°5+2°4 —48+1°6 331°+3° +16°241°2 —24+08 
B (stars) 2°3 +346+4-4 +3:2+3°0 322°+3° +15°0+1'°9 +1413 
B (Ca 11) (1°15) +17°6+2°5 +o2+t2'1 323°+4° +15°342°2 +o2i1°8 
C (stars) 4°0 +32°1+7°4 —38+5°9 327°+7° + 8-0+1°8 —Totis 
= C (Ca 11) (2-0) +24:2+2°1 —2:341°6 332°+3° +121+1'0 —12+08 
kpc km/sec km/sec km/sec/kpc km/sec/kpc 
Errors are standard errors. 
+40 7.5 6 2 i 6.6 2 356 2 3 19'S 
2 8 2 a 3 ' 3 i 2 es 1 18 
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ok Fic. 2.—Mean radial velocities, corrected for solar motion, for stars in Group A reduced 
1e i to a mean distance of 1-5 kpc and shown as a function of galactic longitude. The first order (Oort) 
ts \ least squares solution and the level of the K-term are shown. In longitudes where Lick B velocities 
have been employed (see Section 3) open circles indicate the uncorrected means and the dotted 
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lines the corresponding least squares solution. 
shown at the top of the figure. 


The number of stars contributing to each point is 
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Fic. 5.—Notation as in Fig. 2. Stellar velocities in Group B reduced to 2°3 kpc. 
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Interstellar Ca 11 velocities in Group B reduced to 1-15 kpc. 
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7.—Notation as in Fig. 2. Stellar velocities in Group C reduced to 4:0 kpc. 
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Interstellar Ca 11 velocities in Group C reduced to 2°0 kpc. 
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computing the standard errors of A, and A, no account has been taken of possible 
errors in the mean values of r adopted. Thus the errors quoted are likely to 
be underestimated. 

Whilst a discussion of some of these results is best deferred to a later section, 
when they will be compared with the results derived by other methods, certain 
points may be discussed at this stage. 

‘The mean value of (/, + A) for stars and gas is 327° for the uncorrected solutions 
and 324° for the corrected solutions (group AB not being used in these means). 
However, both in view of the greater number of stars in group A and doubts as 
to the validity of the equation (10) for groups as distant as B and C, it may perhaps 
be best to confine attention to group A. In that case (/,+h) is 329° for the 
uncorrected solutions and 326° for the corrected solutions, again taking the mean 
of gas and stars. It will be clear from the standard errors quoted in Table III 
that the results do not give any evidence for a significant deviation from 327°:7, 
the value of the longitude of the centre recently given by van Woerden, Rougoor 
and Oort (26) from 21 cm studies. Within the accuracy of the results therefore 
we may put h=o. It then follows from equations (11), (12) and (13) that 
Ay=€) and A,=A (the Oort constant) whilst €,’ =o. 

It is of interest to note that in groups A and B the value of Ar is very closely 
half the stellar value for the Catt velocities and in fact this has been assumed to 
be the case in calculating A from the Catt results. The results in group C are 
of too low an accuracy to tell whether the same result holds good there. Whilst, 
of course, there is other evidence (optical and radio) for the concentration of the 
interstellar gas in rather thin (~4 kpc thickness) spiral arms it would appear 
that this has no noticeable effect on the present analysis, which would be satisfied 
by a uniform distribution of gas out to. at least 2-5 kpc. 

The values of A (=A,) given in Table III show a decrease of this constant 
with increasing distance. The low values given by group C are discussed below. 
Setting aside group C it does not appear that the differences between groups A 
and B are large enough to be regarded as significant. From groups A and B, 
gas and stars, we find mean values of A of +15-6km/sec/kpce for the uncorrected 
solutions and +17-8km/sec/kpe for the corrected solutions. In view of the 
discussion of the corrections given earlier and also of the results to be given later 
in this paper the higher value of A is to be preferred. The considerable difference 
between the two results may perhaps partly explain the anomalously low value 
of A (1 3°2) recently derived by Weaver (24) from B stars since the Lick velocities 
requiring correction probably play a considerable part in his analysis. 

It is noteworthy that in view of the inclusion of the new southern results the 
distribution of stars in longitude is much more even than in previous analyses 
of B star velocities and the results should therefore have greater reliability. 
The greatest weight of course attaches to group A, which has more stars and a 
more even longitude distribution than the other groups. We note that, particularly 
in groups B and C, there is an unfortunate deficiency of stars from northern parts 
of the Milky Way. Reference to Table II shows that groups A and B contain 
quite different percentages of dwarfs and supergiants and it is therefore pleasing 
to note that substantially the same value of A is obtained for the two groups. 

We must now discuss the anomalously low value of A found for group C. 
The result must be regarded as unacceptable unless it is confirmed by much 
more extensive data. 
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It can hardly be expected that the Oort approximation will hold for stars 
beyond 3 kpc. The approximation (g) will fail, and further, as will be shown 
in Section 12, approximation (3) also fails for many stars of group C. However, 
failure of these approximations can hardly be the root cause; in Section 10, 
where second-order effects are considered, a similarly low value of A is still 
found for group C. 

It is probable that the result is due to a statistical effect of selection. Many 
of the stars in the group were originally selected for observation on account of 
their large calculated distance moduli; and the group itself represents the 
tail-end of a frequency distribution of calculated distance moduli m—a— M, 
a being written here for the absorption correction. The group is thus selected 
for small a and high luminosity; it will contain a relatively high proportion of 

(i) stars for which the absorption is abnormally small or has been under- 

estimated, 

(ii) stars either truly or erroneously classified as supergiants. 

Errors in m can be neglected compared with errors ina or M. But, with a range 
of some four magnitudes between main sequence stars and supergiants, it seems 
highly probable that the chief contributory cause to uncertainties in the distances 
arises from errors in the luminosity classifications, and possibly to failures in the 
calibration of the MK system. However, it may be noted that the selection effect 
considered here is particularly sensitive to systematic errors in the correction for 
absorption (e.g. in the ratio total absorption to colour excess). If this ratio has 
been underestimated, a will be systematically underestimated, especially in 
group C, and the selection effect will operate to an even greater extent, reinforcing 
errors of over-estimated luminosities; if on the other hand the ratio has been 
taken too large, the errors in a will to some extent balance the errors in M. 

As previously noted, the standard errors of A in the sixth column of Table III, 
ignore the effects of errors in # and are therefore certainly underestimates, especially 
for group C, composed of the smallest sample of stars. 

As regards the Radcliffe material, revised luminosities are being derived by 
Petrie’s method (27) at the Cambridge Observatories through the cooperation 
of Professor Redman. Such measures will provide a useful test of the suggestion 
that the low value of A found for group C arises chiefly from errors in the 
distances. 

6. The K-term.—The K-term was originally introduced as a constant, as 
in equation (15). Clearly the values of K appropriate to each group are the 
same as the values of Apr in Table III. The work of Plaskett and Pearce (1) 
on the A-term may be summarized as follows: 

(a) For the nearby B stars K was positive. 

(6) K decreased with increasing r within the range studied, that is to a 
maximum value of r of about 1 kpc. 

(c) The K-term is probably in the main due to systematic motions of groups 
of stars (especially the Sco—Cen cluster). 

(d) The small residual left after correction for the effect of (c) is probably 
due to Einstein effect in the B stars. 

A recent review of the available data (28) by Weaver substantiates the main 
conclusions of Plaskett and Pearce. 

Confining attention for the moment to group A we find a significant and 


negative K-term. The actual value of Ayr depends of course on whether or not 
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the corrected or uncorrected solutions are adopted but in either case the term is 
significantly different from zero, a value for Ayr (=K) of —4:2+ 1*3 km/sec 
being most likely for these stars. It is therefore interesting to note that the other 
investigations of galactic rotation for distant objects also yield a negative K-term. 
The results of Seyfert and Popper (29) on faint B stars are: 


67 B stars velocities Ar= + 34°34 1°5 km/sec/kpc, 
K=—4:0+ 1-0 km/sec. 
68 Cait line velocities Ar= + 16-6 + 2-6 km/sec/kpc, 
= —3°6+1-8 km/sec. 
(Errors are probable errors.) 


It is worth pointing out that there are only six stars in common between 
Seyfert and Popper’s work and this paper so that the results are virtually 
independent. 

Trumpler (30) has discussed the velocities of 39 galactic clusters in the 
distance range 0-4 to 4°6 kpc. He finds 

A=+15'0+1°0 km/sec/kpe and Ay=— 4:3 km/sec/kpc. 
Trumpler noted that this apparently implied a contraction along the galactic 
plane. 

Joy (31) found a K-term of —3-8 km/sec from his cepheid velocities and 
a value of Ay of —2-4 km/sec/kpc was found by Stibbs (32) in his more recent 
work on cepheids (if /,=328° is assumed). 

Considering these results alone it might appear possible to interpret them 
as indicating a real contraction of material (at least of interstellar gas and young 
stars) in the galactic plane. The linear velocity of contraction at the Sun is then 
AyR, which for Ry=8-2 kpc is about 25 km/sec using the corrected solution. 
It is of some interest to enquire into the relationship of the K-term and the 
well-known deviation of the vertex in star streaming. This is especially so since 
Tannahill (33) has shown that the B stars show the same deviation (about + 15°) 
as stars of other spectral types. A formula relating the constants of galactic 
rotation and contraction (or expansion) to the deviation of the vertex has been 
derived by Shiveshwarkar (34) and when a correction is made for an error in 
sign (Lindblad (35)) it may be used with the present results. In this way, with 
either the corrected or uncorrected set of constants, a deviation of the vertex 
of about +15° is obtained in good agreement with observation. Even though 
Shiveshwarkar’s formula may only be valid under certain conditions of mass 
distribution in the galaxy (see Heckmann (36) and Lindblad (35)), nevertheless 
Lindblad’s qualitative statement that a deviation of the vertex in the direction 
of increasing longitude demands a slight contraction of the material in question, 
would appear to be still applicable. 

The failure to detect a large negative K-term in groups B and C must of course 
greatly weaken the view that the K-term is due to a general contraction of the 
galaxy. Clearly the solutions for groups B and C are of considerably less weight 
than solution A, nevertheless groups containing 95 and 47 stars respectively are 
not negligible. We note furthermore that the results of Seyfert and Popper 
referred to above are, from the values of Ar quoted, comparable in distance 
with group B and yet the K-terms obtained from the two investigations are of 
opposite sign. Moreover, reference to Table III shows that there is no smooth 
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run of Agr with # in our own solutions. Seyfert and Popper’s work is deficient 
in southern stars, while our solutions B and C are deficient in northern stars. 
A poor distribution of stars in longitude may have a very marked effect on the 
deduced value of the K-term as is demonstrated in an actual case by Plaskett 
and Pearce (1). 

If we seek to interpret the K-term by Plaskett and Pearce’s original idea of 
the effect of community motions on the galactic rotation curve it might appear 
that we are now dealing with a sufficiently large volume of space for such 
community motions to be attributable to larger groupings of stars than individual 
clusters. In particular we might bear in mind in this connection the possible 
effect of spiral structure and the possibility of the radial ejection of gas and stars 
from spiral arms. The recent radio discovery (26) of a spiral arm near the centre 
of the galaxy, expanding outwards with a velocity of 50-100 km/sec, clearly 
shows that at least in certain regions quite large radial motion in the galaxy may 
be expected. 

7. Velocity dispersion.—For stars in the region of 1-2 kpc (group A) and for 
those 10° longitude intervals containing more than four stars each, the mean 
dispersion in residual velocity has been calculated. The results which are given 
in Table IV are uncorrected for observational errors (including undetected 
binaries and errors in the adopted distances). Whilst there is a considerable 
scatter in the results it is not felt that this is significant. The particularly large 
values at 44° and 94° are affected by one or two discrepant velocities. 

The mean dispersion of 14°8+2:2 km/sec (using the corrected Lick B 
velocities) is not far from the value of 12 km/sec previously derived from rather 
different data (2). The figures in Table IV give no support to Weaver’s 


TaBLe IV 
Velocity dispersion for stars in Group A 
1 Weight 09(km/sec) 
44° 8 34°7 
53° 5 15°6 
86° 6 7°3 
94° 4 44°0 
103° 12 II*4 
155° 7 5°0 
174° 18 13°8 
193° 6 23°4 
246° 6°5 I1°4 
255° 9°5 17°0 
264° 6°5 19°7 
274° 9°5 17°3 
336° 19 Ye 
343° 18 9°5 
352° 5° 17°5 
mean 69 14°8+2°2 


‘ 


contention (37) that the velocity dispersion is high for stars ‘‘in the inner part 
of the galaxy’. Weaver finds for a small group of stars (14 in all) lying between 
271° and 338° with one star at 37°, a velocity dispersion of 36 km/sec ‘‘a value 
nearly 5 times as great as that observed in the outer and local arms’’. From 
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Weaver’s figure it would appear that most of his stars lie between 1 and 2 kpc 
from the Sun (near 1 kpc). But although our stars cover this same region we 
have only one star in common. 

In another paper (38) Weaver obtains a velocity dispersion of 35-9 km/sec 
from 17 supergiants in this same region (presumably including the 14 stars used 
in (37)). For 87 B stars between 330° and 350° he finds a dispersion of 
230 km/sec which again is taken as indicating an abnormally high velocity 
dispersion in this region. ‘The standard error of this value is not given but is 
probably at least as large as ours, which depends on nearly twice the number of 
stars (172). In view of this and other possible causes of uncertainty in the data 
available to Weaver, it would appear necessary to obtain a good deal more data 
before any variation of velocity dispersion with longitude can be firmly established. 
Our results are in good agreement with those of Nordstrém (39) who found for 
353 stars earlier than B5 and brighter than 6™-o a velocity dispersion of 
8-8 km/sec and for 379 stars earlier than Bs and fainter than 6"-o a dispersion of 
14°7 km/sec. His group of fainter stars is obviously more nearly comparable 
with ours and the dispersion in this group agrees very well with our value of 
14°8 km/sec. 

In view of the errors in radial velocities and distances associated with fainter 
stars, we do not consider that the difference between the velocity dispersion of 
B stars brighter and fainter than 6™ should be regarded as significant without 
a great deal of further investigation. 

8. The application of the Camm linear approximation.—We now proceed to 
analyse the data on the basis of the Camm linear approximation—equation (7) 
above. ‘The most recent use of this equation is the analysis by Stibbs (32) of 
cepheid radial velocities. He assumed R, = 8-2 kpc and solved equation (7) for 
two different values of /,: 325° and 328°. ‘The recent radio work (26) has shown 
328° to be probably the better of the two values so that we did not consider it 
necessary to carry out the solution for /4)=325° for the B stars. However, as 
the value of R, is still subject to considerable uncertainty, it was felt necessary 
to investigate the effect of a change of R, on the solution and for this reason the 
equations were solved for values of Ry of 8-2 kpc and 12-0 kpc. 

Instead of computing R for each star from equation (8) its value was measured 
on a diagram showing the stars projected on the galactic plane. In view of the 
various errors in the determination of 7, the values of R thus found are of more 
than sufficient accuracy. Since equation (7) is only valid for |R—R,| small 
we used in the analysis only those stars for which |R— Ry| < 1-0 kpe (cf. Fig. 1). 
This reduces the number of stars to 134 in the case of Ry = 8-2 kpe and 131 in 
the case of Ry=12-0 kpc. Of these stars 127 are common to both solutions. 

In the R=8-2 kpc solution there are only four stars which depend wholly 
or partly on the Lick B velocities discussed in Section 3, whilst in the Ry = 12-0 kpc 
solution there are five such stars. ‘The analyses were carried out without applying 
the corrections to these velocities which are now considered desirable. Since 
however there are so few stars requiring correction the final results will not have 
been appreciably affected. 

Equation (7) was solved by least squares to give the constants and their 
standard errors listed in ‘Table V which also contains the corresponding 
(J, = 328°) solution by Stibbs (32) for the cepheids. In applying the method 
of least squares the assumption was made both in our work and in that of Stibbs 
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” that the principal source of error is in the velocities (observational error plus 
a velocity dispersion) ; the errors being considered to be relatively less important 
iny and (R—R,). It is not known how far this assumption is justified. Whilst 
” it would, theoretically, be possible to compute regressions of both V on r and 
d ron V it would be difficult in practice to decide how to weight the individual 
of stars in each case and how to determine a weighted mean of the two regressions. 
ry Furthermore, in the case of r, systematic errors may well be important and thus 
i a real improvement on the present method of analysis does not appear practical. 
of However the limitation of the least square method used should be clearly recog- 
ta nized since, in any case, it may quite probably have led to an underestimation 
“i of the standard errors quoted in lable V. It appears from ‘Table V that the two 
j solutions do not differ significantly. The constants are thus rather insensitive 
* to the adopted value of Ry in this range. 
of 
of TABLE V 
le Results of least square analysis by Camm’s linear approximation 
of Solution 1 Solution 2 ome (32) 
(cepheids) 
sal R, (adopted) kpc 120 8-2 8-2 
of 1, (adopted) 328° 328 328° 
ut €9 (km/sec/kpc) — £2 269 — 20 t09 = 24. £¥3 
A (km/sec/kpc) +570 £15 Pry 3 EES +186 +2:°0 
~ C (km/sec/kpe) + 4°7 +2°4 + 3:4 +2°5 + oF £279 
7) Number of stars 
of From Group A 88 87 
or From Group B 31 35 
vn From Group C 12 12 
it 
Total é 131 134 96 
as 
ry Errors are standard errors. 
he 
We see that a A-term is present (non-zero values for €g and C). It will be 
-d remembered that a negative value of Ay (=€,) was found in the Oort analysis for 
he stars in the range 1 kpc<r<2kpc. Since the stars in the present analysis are 
re mainly from this group (cf. Fig. 1), but only of course a selection from this group, 
all the negative value of €, is not surprising. It will be recalled that the results of 
[). the more distant groups in the Oort analysis did not support the interpretation 
in of the K-effect as due to a dilation of the galaxy. Confirmation of this view may 
1S. be obtained from a consideration of the residuals in the two Camm type solutions. 
ly In both cases the residuals are not randomly distributed with respect to r but 
pe are predominantly positive for r < 1-4 kpc and predominantly negative for greater 
ng values of r. This would appear to confirm that the K-effect cannot be represented 
ce in the form used in equation (7) (general galactic expansion or contraction). 
ve | We may also note that application of formulae (6) and (13) leads to a value of h 
; of about +15” if the values of C in Table V are taken seriously. ‘This would not 
eir accord at all well with the Oort analysis given above which indicated h=o. ‘This 
ng may be taken as another indication that the K-effect is to be interpreted otherwise 
od | than by a general galactic dilation. 
bs | 
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9. The value of Oort’s constant A.—The Camm linear approximation, 
which we regard as the most reliable approach to the problem of the value 
of the Oort A constant, yields A=17-6+1°5 km/sec/kpc for R)=12-0 kpc 
and A=17-3+1°5 km/sec/kpe for Ry=8-2 kpc (Table V). With greater weight 
to Ry=8-2 we adopt A=17'4+1°5 km/sec/kpc. In an earlier section it was 
shown that the Oort approximation gave a value of 15-7 km/sec/kpc for Oort’s 
constant A if uncorrected values of the Lick B velocities discussed in Section 3 
were used. On the other hand the corrected velocities give a value of 
17°8 km/sec/kpc. As pointed out above, the Camm-type analysis happens to 
be essentially free of stars requiring this correction. Hence the fact that it gives 
a value for A close to the value for the “‘corrected’’ Oort solutions may be 
regarded as further evidence for the validity of the applied corrections. 

The present work giving greatest weight to the Camm analysis leads us to 
adopt 

A=+17°5+1°5 km/sec/kpc 
as the most probable value indicated by the analysis of the present data. We 
may now compare it with other recent determinations. 

Table VI lists several recent determinations of A. Stibbs’ value from the 
cepheids, though slightly higher than ours, cannot be considered as being 


TasLe VI 


Recent determinations of Oort’s constant A 


A (km/sec/kpc) 


and s.e. 
Present paper Camm analysis, 134 B-stars I7‘4+1°5 
Oort analysis 
267 B-stars 
227 Ca 11 velocities 17°8 + (1°5) 
Adopted mean I75H1°5 
Stibbs (32) Camm analysis, 96 cepheids 18-6+2°0 
Oosterhoff (40) Re-discussion of (32) 17*4+2°1 
Gascoigne and Eggen (41) Camm analysis, 37 cepheids 17*5+1°9 
Petrie et al. (42) 79 B-stars 
as 68 Ca 11 velocities 17°77 + 16% 
Weaver (23) 79 cepheids 10°8+1°8 
Weaver (24) 187 B-stars 13'2+1°8 
Weaver (24) Adopted mean IVs +1°5 
Bahng, Code and ' 
Whitford (43) B-stars and cepheids IIS +1°5 


* s.e. computed from published p.e. 


significantly different from it. Furthermore recent rediscussions of some of 
Stibbs’ velocities by Oosterhoff (40) and Gascoigne and Eggen (41) using more 
accurate distances for the stars, deduced from photoelectric magnitudes and 
colours, leads to a value of A practically identical with ours. It is of considerable 
interest to see that our value of A is in close agreement with that deduced by 
somewhat different methods by Petrie, Cuttle and Andrews (42) from quite 
independent data on B stars. These values all contrast markedly with the values 
obtained by Weaver (23, 24) and by Bahng, Code and Whitford (43) for B stars 
and cepheids. Both Petrie et al. (42) and Weaver (44) discuss the difference 
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in their values of A for the B stars and it has been suspected that this difference 
might partly be due to Weaver depending for absolute magnitudes on MK types 
and Petrie et al. on Petrie’s Hy method (27). However since we have based 
our analysis on MK classifications and yet are in good agreement with Petrie et al. 
it would appear rather unlikely that this is a major cause for Weaver’s discrepant 
value. 

We may conclude that, so far as present data on B stars and on cepheids goes, 
a value of A close to 17°5 km/sec/kpce is preferable both to recent anomalously 
low values and to the value of 19-5 km/sec/kpc adopted in reducing the 21 cm 
results (45, etc.). 

10. Distance to the galactic centre.—Referring to Section 4, equations (16) 
and (18) have been used in two rather different approaches to the problem of 
determining R, from the present material. 

(a) The symmetry of the function V(A) in equation (16) allows the grouped 
mean velocities as a function of A to be formed in the following way. Stars in 
a particular range A, to A, are grouped together with those between —A, and 
—A, (velocities of the latter being taken with negative sign). It is important to 
note that this treatment eliminates a general K-term so long as the weights of 
stars between A, and A, and between —A, and —A, are the same. 

The form of the function V(A) is not amenable to an analytical treatment by 
least squares. Instead the sums of the squares } p(O-C)? were directly computed 
for various values of the parameters A, Ry and x. The values of the three 
parameters giving minimum > may be regarded as the most probable, while the 
significance of the solution may be assessed by finding the range of other values 
which give 


ve I 
¥< Some Dd min 


where N is the number of observed mean velocities used, and m is the number 
of degrees of freedom (i.e. number of disposable parameters) (46). 

In the computations, /, was regarded as a known quantity, J, =327°-7. For 
the ratio of the axes of the velocity ellipsoid, Hins and Blaauw (47) have found 
the value 0-49 + 0-04 for faint low-latitude stars. Consequently, for most of the 
computations we have adopted 

K=0°5. 


However, as will appear later, adoption of other values of « in certain cases did 
indeed show that no significant diminution of } was achieved and therefore 
that «=0°5 is likely to be close to the true value. 

An interesting feature of the method is that the product AR, which appears 
as a scale factor in equation (16) can be compared directly with the 21 cm radio 
observations of the interstellar gas. Schmidt’s analysis (48) has given a radio 
value AR, = 156 km/sec. 

Preliminary applications. As a first approach to the problem, the identical 
material of groups A, B and C already analysed for first-order effects were 
subjected to the above treatment. It will suffice to describe the results in brief. 

Group A. min was found to diminish steadily from p=o0-18 to p=0-075. 
The significance test suggests that Ry > 13 kpc, and a value A = 18-2 is derived. 
The method fails, presumably due to: 

(i) the smallness of second-order effects, 





144 M. W. Feast and A. D. Thackeray Vol. 118 


(ii) the presence of a large negative K-term which is not completely balanced 
out in the process of taking means over a 180° range of A, 

(iii) the probable presence of systematic errors in observed velocities at 
certain A. 

‘The Lick velocities carry considerable weight in this group, and if they are 
corrected by + 10 km/s > is decreased by a factor 2. We regard this improvement 
of fit as a further argument in favour of applying this correction to the Lick 
velocities of faint B stars. The effect of the corrections is to increase A to 20°5, 
but > still diminishes (rather more slowly than before) from p=0-18 to 0:075. 

Finally this group was re-analysed using corrected Lick velocities and a 
K-term, K=— 4:0 km/sec, as suggested by the first-order analysis. In this 
instance we find, with «x =o0'5, 


Lmin Ro A AR 
3106 8-24 21°I 174 
2991 14°0 21°0 294 
3058 20°0 20°6 412. 


For the first time a minimum is found with R,<20, but the flatness of the 
minimum is such that no statistically significant distinction can be drawn between 
the three solutions. All three values of AR, are larger than that given by the 
radio observations. We conclude that this group is unsuitable for a determination 
of Ry, although it does tend to indicate a distance to the galactic centre greater 
than 8 kpc. 

Group B. With «=o'5 and uncorrected Lick velocities, Sin is found with 
p=o-20 and A=14:6. Thus, the mean distance of the group being 2-30 kpc, 
we find Ry= 11-5 and AR,=168. The significance test suggests that R, > 8-5 kpc. 

If the Lick velocities are corrected-by +10 km/sec > is diminished by about 
30 per cent, the most probable value of A is somewhat increased, and that of R, 
diminished. 

Group C. minis found with p=o0-40 and A=g-g9. ‘Thus Ry = 10°0 is indicated 
but the significance test shows that R, = 8-o is not excluded. ‘The low value of A 
found parallels that found in Section 5. 

All three above analyses suggest a value of Ry greater than 8 kpc. It is clear 
that the method has to deal with uncertainties due to different causes in the three 
groups. In group C, where the second-order effects are largest, one has to 
recognize the greatest errors in calculated distances. In group A the second-order 
effects are rather top small, despite a comparative wealth of material, to derive 
much idea about the true value of R,; besides, as we have seen, the region 
occupied by group A appears to be influenced by a marked general K-term of 
about — 4 km/sec. 

Final application: group BC’. For a final determination of what may be 
regarded as the best value of R, attainable from this material by the method 
described, we have formed a new group BC’ of 113 stars with the limits r= 2:0 
to 3°6kpc. The upper limit of 3-6kpc will exclude many over-estimated 
distances. The velocities were corrected by a factor to reduce them all to the 
mean distance #= 2°54 kpc. In accordance with previous experience of improved 
fit, the Lick velocities were corrected by +10km/sec. Moreover, Lick and 
McDonald observations were given half-weight. The resulting velocities were 
grouped into 15 means, tabulated in Table VII where \’=|/—327°-7|. In the 
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third column, p, represents the weight of the stars in the range (/—/,)=0 to 180°, 
and in the fourth p_ represents the weight in the range (/—/,) = 180 to 360°. The 
total weight appears in the fifth column, and the number of stars in the sixth. 
If a general K-term is to be applied to the material, then we must correct V by 


adding (p_—p,)K'/p. 


TasBie VII 
Mean velocities for Group BC’ 
A’ V P+ p- p n A 
8-4 +17°7 6 ° 6 10 +1°8 = 
16°7 +24°0 tj I 8 14 +3°2 7 
33°6 +51°5 2 4 6 7 +32 . 
47°8 + 30°9 ° 6 6 | +10 
55°3 +3671 ° 142 144 15 —y 
64°2 + 18-7 1} 4 83 10 —2°0 . 
726 +16°2 13 9 104 12 —3°0 . 
84°5 + 10°5 4 4 I 2 —34 sf 
95'8 —13°0 I I 2 2 — oo . 
104°9 —29'1 44 3 74 8 -— : 
117°6 —40°9 4 ° 4 4 =F * 
127°0 —28°7 6 ° 6 6 —o'2 
135°4 — 37°4 7 I 8 9 +0°5 
153°6 — 34 3 2 24 4 +0°9 
172°6 —II'5 ° 23 24 3 +03 


In Table VIII (a) appear computed means of }min with «=0'5 and a series 
of values of Ry, while in Table VIII (5) appear values when p=o0-30 is assumed 
together with various values of x. From the latter it is clear that «=1-o is 
excluded, although such a value has been suggested as appertaining to main 
sequence B stars (16, Tables 3, 5, p. 332). Nordstrém (39) however gives a 
value 0-59 for O to Bs stars fainter than 6"-o. The present material adds support 
to the value «—o-5 as appropriate to distant stars of B type, in agreement with 
Hins and Blaauw’s value previously quoted (47). 

Turning to Table VIII (a), with «=0-5 assumed, we find the best fit with 
p=0'29, A=15°4, Ry =8-76 and AR, = 135. 

According to the significance test, we find permissible values of p from 0-21 
to 0°365, or Ry= 12:1 to 7:0 kpc. 


TasB_e VIII 
Determination of minimum Sp(O-C)? for group BC’ 
(a) K=0°5 
p Ro A min AR» 
0°20 12°70 15‘! 3976 IgI 
0°25 10°16 15‘I 3682 153 
0°28 9°06 15°3 3607 139 
029 8-7 15°4 3604 135 
0°30 8°47 15°5 3606 131 
0°35 7°26 15°5 3790 122 
(6) p=0'3, Ro=8°47 
K A win AR> 
04 15°5 3639 131 
O'5 15°5 3606 131 
06 15°6 3849 132 
| exe) 14°8 10,722 125 


10 
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The observations of Table VII are plotted in Fig. 9, where the curve shows 
the predicted velocities based on equation (16) with the above parameters, 
Note that the node occurs at A=82°, while the higher amplitude of the curve 
in the left-hand part compared with the right-hand part forms the basis for the 
value «=0°5. Increasing the value of « towards unity has the effect of depressing 
the relative amplitude of the left-hand portion. 





+60 T T T r v T T T 


+40 


+20 


VELOCITY (km/sec) 
° 














4 4 4 4 ———EE 1 
° 


- EE 
100° 20° 140° — 160° 


1Al 
Fic. 9.—Velocities of 113 stars (r=2-0 to 3°6 kpc) fitted to equation (16) with Ryo=8'8 kpc 
and x=o'5. Filled circles are based on corrected Lick B velocities (Section 3). Open circles are 
based on uncorrected velocities. For numbers of stars contributing to each mean see Table VII. 
Note the greater amplitude of the curve for the region |A|< 80° than for |A|> 80°. 
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180 


The foregoing analysis depends on fitting all 15 means of Table VII to the 
theoretical curve, with Ry as a disposable parameter. As is well known, one of 
the chief effects of increasing p is to shift the node from A=g0° to smaller values. 
If the node A, can be determined sufficiently accurately, we obtain a direct measure 
of R, from the relationship 

Ry) = 47 secA, 
which is strictly accurate for circular orbits. Other things being equal, it is the 
ring of stars at distance R, from the galactic centre which can provide the best 
determination of R, from an analysis of radial velocities. It would therefore 
be legitimate in the present analysis to confine attention to those stars near go 
from the galactic centre. Sucha restriction of the material is adopted in method (4) 
below, but here as a compromise we have tried the following alternative. 

The last column of Table VII gives the difference in velocities computed for 
the curves of best fit based on p=o0-29 and p=o-20. Those entries where the 
difference exceeds 10 km/sec are marked with an asterisk; there are nine such 
entries which may be regarded as “‘ sensitive directions’’ for the determination 
of Ry. If we use only these nine entries we find the following values of }min a° 
a function of p: 


Pp Xin 
0°20 1600: 
0°28 1314 
0°29 1318 
0°30 1340 


"35 159° 
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The most probable value of Ry is 9:06 kpc. According to the significance test, 
values of }min up to 1314 7/6, i.e. 1533 are permissible, whence values of 
p=0'214 to 0°341 or Ry=11'g9 to 7-45 kpc are permissible. It is interesting to 
note that the range of values of R, in this solution, depending on a smaller weight 
of observations, has been narrowed by 0-6 kpc. 

(b) In using equation (18) we have employed a modification of a method 
due to Weaver (22). 

Weaver concentrated attention entirely on distant stars for which |R— Ro} 
is small, and the velocities of such stars are undoubtedly the most sensitive to R, 
independently of w’ and w”. If for any star or group of stars we know J, I, b, 
and r and adopt a value of /, (327°-7), we may calculate for a set of assumed values 
of Ry a corresponding set of values of w,’. If then we do this for several stars 
or groups of stars and plot the results in the Ry, wy’ plane we obtain a set of 
curves which should all theoretically intersect at one point which gives the 
required true values of Ry and w,’. Since it is obvious that only stars with large r 
can give a value of Ry of any accuracy, we restrict ourselves to r>3-0 kpc. Also 
since (18) uses the Camm linear approximation, we have the restriction 
|R—R,|<1t-okpe. Finally since for (R—R))+0, Vo, w,' is indeterminate 
(from (18)), we have the further restriction that |R—R,|>0°5 kpc. Satisfying 
these criteria there are only nine stars, all lying between /=256° and 263°, and 
with r ranging from 3°5 to 4:8 kpc (cf. Fig. 1). These nine stars define the mean 
curve shown in Fig. 10. The dotted lines indicate the uncertainty introduced 
by the standard error of V. The second curve in Fig. 10 is drawn according to 
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Fic. 10.—The determination of Roy by the second method of Section 10. The two solid lines 
show the Ro, wo’ relationships defined by 9 distant stars and by a value of Oort’s constant of 
17°5 km/sec/kpc. Ry is determined by their point of intersection. The dotted lines show the 
standard errors of the curves. 


the condition A = — }Rgw9'=17°5 km/sec/kpc. A is well determined, indepen- 
dently of Ry, from many relatively nearby stars (the Camm solution in Section 8 
depends essentially on stars with r between 1 and 2kpc). Again the dotted 
curves give the uncertainty indicated by the standard error of A. 
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We find that the most probable value of Ry is 8-7 kpc, while the permissible 
range of values is from 7°6 to 11-3 kpc. 

This application of the method depends on the comparison of the motions 
of many relatively close stars with the motions of nine distant stars in the most 
sensitive direction. 

It is of interest to note that of the nine stars used in method (6), in conjunction 
with the value A=17-°5 km/sec/kpc, only four are common to the group BC’ 
used in method (a). Hence the two methods may be regarded as largely 
independent. 

In ‘Table IX we compare the results for Ry found by the two methods. 


TaBLe IX 


Determination of Ro (kpc) 


Most Range of 
probable values 
value 
Method (a): 
BC’: all stars 8-8 7°O to 12°1 
BC’ restricted g'I 7°4 to II'9 
Method (b): 
A=17'5, 9 Stars 8-7 7°6 to 11°3 
Mean 8-9 


‘To summarize the results of this section: 

(i) the most probable value of Ry is 8-9 kpc, based on the MK scale of 
distances, 

(ii) method (a), dependent on an assumed force law for the galaxy, confirms 
the result found in the Oort analysis that the most distant stars give an abnormally 
low value of A; this is probably to be attributed to errors in the distances of many 
stars in group C (Section 5). 

(iii) a value of « =0°5, for the ratio of the axes of the velocity ellipsoid, appears 
to be appropriate to distant B stars. 

11. The use of B star radial velocities to determine absolute magnitudes.—If we 
assume that the radial velocities of the B type stars, corrected for solar motion, 
follow the simple Oort formula (14) then we may use the velocities to derive the 
relative absolute magnitudes of the various sub-groups of B stars. This follows 
since for each star we may determine (7A) from (14) and then reduce these values 
to some fixed apparent magnitude. If (rA), is the value of (rA) reduced in this 
way then clearly log (rA), is proportional to the absolute magnitude of the star. 
If a value of A is assumed, the actual absolute magnitude can be deduced. In 
carrying out this calculation we restrict ourselves to stars with |sin 2A|>0°5. 
There are 214 stars between Bs and Os which satisfy this condition. Adopting 
A=17'5 km/sec/kpc we obtain the results shown in Table X. The fourth column 
of this table gives the normal calibration of the MK system (10). 

Whilst the main purpose of the present paper is to use the B stars with an 
assumed absolute magnitude system for the determination of galactic motions, 
it is not without interest to reverse the arguments to see that the results obtained 
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are consistent with those normally adopted. Such an investigation is of especial 
interest in that Blaauw (49) has determined absolute magnitudes of main sequence 
stars of types B1 to B3 and our material may be used to extend the comparison 
with the MK calibration to still earlier types. As noted above, the relative values 
of M in the table are independent of any assumed value of A. 


TABLE X 


Calibration of B star absolute magnitudes using the Oort approximation and 
A=17'5 km/sec/kpc. 


Spectral Type No.of Stars Weight Normal MK Deduced 

Abs. Mag. Abs. Mag. 
(10) 
Main Sequence (Luminosity Classes IV and V) 

m m mh 
B2 to Bis 15 14 —3°1 —36 + 04 
Bi to Bors 30 23°5 —36 —4'2 + 03 
Bo to Og°5 8 6°5 —4°4 —46 + 06 
Early O types 29 26°5 —5'I —4°8 + o-2 

Giants (Luminosity Classes IT and III) 

mh mm mn 
Bs to B3 6 4°5 —4°0 —4°7 + 06 
B2 to Br-5 8 | —4°6 —3°1 + 06 
Br to Bors 27 22°5 —4°8 —4°3 £03 
Bo to O9°5 12 fe) (—5°4) —5'4 + 05 
O9 5 4°5 (—6°4) —43 07 

Supergiants (Luminosity Class 1) 

m m Mm 
Bs to B3 * #8 14 —6°3 —5°8 + 04 
B2 to Br'5 15 15 —6°4 —6:2 + 04 
Br to Bors 26 24'5 —6°4 —6°4 + o-2 
Bo to O95 15 15 —6°4 —5°7 + 0°5 
Og 3 4 —6°4 —6°3 + I°0 


Errors are standard errors. 


It will be noticed that we have not restricted ourselves to nearby stars as 
might perhaps have been considered necessary in view of the approximate nature 
of equation (14). However the results of the Oort type analysis given earlier 
shows that deviations from (14) for nearby stars, due to K-terms, are likely to 
be quite as important as deviations from it for distant stars due to ‘‘ second-order ”’ 
effects. ‘Thus stars are included irrespective of distance. The corrected Lick 
B velocities (Section 3) have been used. In general it can be said that the 
agreement between the last two columns of Table X is remarkably good. 

12. The differential rotation of the B stars.—If we neglect any effects due to 
a general expansion or contraction of the galaxy then the exact equation for the 
effect of differential rotation on radial velocities is the Bottlinger equation (2) 
above. The first attempt to use this formula to obtain (w — wa») as a function of R 
was apparently made by Camm (18, 19) using the data available to him on 
planetary nebulae and cepheids. Recently the method has been applied to 
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cepheids by Weaver (23) and Stibbs (32) and by Petrie, Cuttle and Andrews (42) 
to B stars. A determination of (w—w,)) as a function of R has also been made 


from 21 cm observations (45, 50). 


These latter observations are the basis of 


the theoretical model of mass distribution in the galactic system deduced by 
Schmidt (48). It is therefore of interest to analyse the present data in the same 
way for comparison with these previous investigations, notwithstanding the 


possibil 


+24 


+2\ 


+18 


+9 


+ 6 


(w— wo) (km/sec/kpc) 


ities of radial motion in the galaxy discussed earlier. 

















R (kpc) 





Fic. 11.—Relationship between R and (w—w ) for 218 stars. 
for A=17°5 km|sec/kpc (w’ =w' =constant). 
(48). 


The straight line is drawn 


The curve is drawn according to Schmidt’s model 
A value of Ry=8-:2 kpc is assumed. Some points represent more than one star. 
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In analysing the data according to equation (2) a longitude of 327°-7 was 
adopted for the galactic centre and a value of 8-2 kpc for Ry. This value of Ry 
was used so that a direct comparison with previous work could be made. All 
stars with |sinA cos b| < 0-23 were rejected from the analysis. ‘This left 219 stars 
of which 186 have measured velocities for interstellar Catt. Figs. 11 and 12 
show the relationship between (w—w,) and R in the case of the stars and the 
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Fic. 12.—Relationship between R and (w—wy) for 186 interstellar Ca 11 velocities. The 
straight line is drawn for A=17°5 km|sec/kpc. (w'=wo'=constant). The curve is drawn accord- 
ing to Schmidt’s model (48). A value of Ro=8-2 kpc is assumed. Some points represent 
more than one velocity. The six open circles represent 3 stars in each of which double interstellar 
lines have been measured. 
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interstellar gas. In the case of the gas the value of R was deduced by assuming 
r(Cat1)=r(star)/2. The small number of points near R=R, is due to the 
selection of the original data forr > 1-oo kpc. Ifthe same value of Oort’s constant, 
A, applied throughout the region investigated, the relationship between (w — wy) 
and R would be linear. The straight lines in Figs. 11 and 12 are drawn for 
A=17'5 km/sec/kpc. The curves in the two figures show the expected relation- 
ship on Schmidt’s model (48) of the galaxy and essentially represent the 21 cm 
observations. The figures give some indication that Schmidt’s model is a better 
fit to the points than a linear relationship. Whilst there is only a slight indication 
of this, it may at least be claimed that the observations discussed in this paper 
are not inconsistent with Schmidt’s model. This is a point of some interest since 
the curves obtained by Weaver (23) and Stibbs (32) for cepheids and by Bahng, 
Code and Whitford (43) for B stars agree neither with each other nor with Schmidt’s 
model. In the case of Weaver’s curve w” is negative. This point has been 
discussed by Schmidt (48), who concluded that such a state of affairs was very 
unlikely on theoretical grounds. In the case of Stibbs’ results the deviations 
may be connected with errors in the distances (r) used in the analysis. Lacking 
photoelectric colours and magnitudes for the majority of his stars, Stibbs had 
to depend on more approximate methods of determining their distances. This 
may well have led to fairly large systematic errors for the more distant stars 
which are the ones that show the deviations from Schmidt’s model. In fact a 
small increase in the average extinction per unit distance in the galaxy over that 
adopted by Stibbs would probably reduce considerably the discrepancy between 
his results and the theoretical curve. The results of Bahng, Code and Whitford 
have so far appeared only in abstract and it is impossible therefore to make a 
proper comparison between our results and theirs. However our results appear 
to be based on more extensive data. We cannot therefore regard any differences 
between the angular velocity curves derived from stars and 21 cm observations 
as definitely established. The results of Petrie, Cuttle and Andrews (42), 
which cover a somewhat smaller range of R than the present investigation, appear 
to be in quite satisfactory agreement with our results. 
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also to Mr L. Tilbury for calculations of the corrections for solar motion used 
for some stars. 


Radcliffe Observatory, 
Pretoria: 
1957 December. 


References 


(1) J. S. Plaskett and J. A. Pearce., D.A.O. Publ., 5, 1931-1936. 

(2) M. W. Feast, A. D. Thackeray and A. J. Wesselink, Mem.R.A.S., 67, 51, 1955. 

(3) M. W. Feast, A. D. Thackeray and A. J. Wesselink, Mem.R.A.S., 68, 1, 1957. 

(4) R. E. Wilson, General Catalogue of Stellar Radial Velocities (Carnegie Institution), 
1953. 

(5) W. W. Morgan, A. D. Code and A. E. Whitford, Ap.¥. Sup., 2, 41, 1955. 

(6) W. A. Hiltner, Ap.F. Sup., 2, 389, 1956. 

(7) P. Th. Oosterhoff, B.A.N., 11, No. 425, 1951. 

(8) W. W. Morgan, D. L. Harris and H. L. Johnson, Ap.7., 118, 104, 1953. 

(9) H. L. Johnson and W. A. Hiltner, Ap. 7., 123, 267, 1956. 











118 


ing 
the 
int, 
wo) 
for 


cm 
tter 
ion 
per 
nce 
ng, 
it’s 
een 
ery 
ons 
ing 
rad 
his 


ars 


hat 
een 
ord 
ea 
ear 
ices 
ons 
2), 


ear 
nan 


rk; 
sed 


on), 








No. 2, 1958 Analysis of radial velocities of distant B type stars 153 


(10) P. C. Keenan and W. W. Morgan (in) Astrophysics, ed. J. A. Hynek (McGraw Hill),. 
1951. 

(11) W. M. Smart, Charts issued by the R.A.S. 1923. 

(12) D. W. N. Stibbs, M.N., 115, 363, 1955. 

(13) D. S. Evans, A. Menzies and R. H. Stoy, M.N., 117, 534, 1957. 

(14) R. M. Petrie, D.A.O. Publ., 9, 297, 1953. 

(15) F. J. Neubauer, Ap. 7., 97, 300, 1943. 

(16) R. J. Trumpler and H. F. Weaver, Statistical Astronomy (University of California 
Press), 1953. 

(17) F. K. Edmondson, P.A.S.P., 67, 10, 1955. 

(18) G. L. Camm, M.N., 99, 71, 1938. 

(19) G. L. Camm, M.N., 104, 163, 1944. 

(20) J. Titus, A.7., 49, 1, 1940. 

(21) S. Chandrasekhar, Ap. F., 90, 131, 1939. 

(22) H. F. Weaver, A.f., 59, 375, 1954. 

(23) H. F. Weaver, A.7., 60, 201, 1955. 

(24) H. F. Weaver, A.7., 60, 208, 1955. 

(25) H. F. Weaver, A.7., 60, 211, 1955. 

(26) H. van Woerden, W. Rougoor and J. H. Oort, C.R. Acad. Sci. Paris, 244, 1691, 1957. 

(27) R. M. Petrie, D.A.O. Publ., 9, 251, 1952. 

(28) H. F. Weaver (in) Vistas in Astronomy, Vol. I, ed. A. Beer (Pergamon Press), 1955. 

(29) C. K. Seyfert and D. M. Popper, Ap. 7., 93, 461, 1941. 

(30) R. J. Trumpler, Ap. 7., 91, 186, 1940. 

(31) A. H. Joy, Ap.F., 89, 356, 1939. 

(32) D. W. N. Stibbs, M.N., 116, 453, 1956. 

(33) T. R. Tannahill, 17.N., 114, 460, 1954. 

(34) S. W. Shiveshwarkar, M.N., 95, 655, 1935. 

(35) B. Lindblad, M.N., 96, 69, 1936. 

(36) O. Heckmann, M.N., 96, 67, 1936. 

(37) H. F. Weaver, A.7., 58, 177, 1953. 

(38) H. F. Weaver, P.A.S.P., 65, 132, 1953. 

(39) H. Nordstrém, Med. Fr. Lund. Ast. Obs. Ser. II, No. 79, 1936. 

(40) P. Th. Oosterhoff, Reported at the Conference on the Coordination of Galactic Re- 
search, Stockholm, 1957. B.A.N., 14, No. 484, 1958. 

(41) S. C. B. Gascoigne and O. J. Eggen, M.N., 117, 430, 1957. 

(42) R. M. Petrie, P. M. Cuttle and D. H. Andrews, A.7., 61, 289, 1956. 

(43) J. D. Bahng, A. D. Code and A. E. Whitford, Sky and Telescope, 16, 529, 1957. 

(44) H. F. Weaver, Reported at the Symposium on the Galactic Distance Scale, Cambridge, 
Mass., 1957. 

(45) H. C. v.d. Huist, C. A. Muller and J. H. Oort, B.A.N., 12, No. 452, 1954. 

(46) A. Pannekoek and E. van Dien, B.A.N., 8, No. 297, 1937. 

(47) C. H. Hins and A. Blaauw, B.A.N., 10, No. 391, 1948. 

(48) M. Schmidt, B.A.N., 13, No. 468, 1956. 

(49) A. Blaauw, Ap.7., 123, 408, 1956. 

(50) K. K. Kwee, C. A. Muller & G. Westerhout, B.A.N., 12, No. 458, 1954. 





STELLAR GROUPS. II. THE ¢ HERCULIS, « INDI 
AND 61 CYGNI GROUPS OF HIGH-VELOCITY STARS 


Olin F. Eggen 
(Communicated by the Astronomer Royal) 


(Received 1958 January 1) 


Summary 


A catalogue of 700 high-velocity stars has been searched for objects that 
share the well-determined space motions of ¢ Herculis (22 stars), € Indi 
(14 stars) and 61 Cygni (15 stars). The space motions of these groups can 
be summarized as follows (the Hyades and Sirius Groups of Paper I have 
been added for comparison): 


Group A D V; U V W 

h m ° km/sec 
Sirius 20 44 —42°7 18°4 —14 ° —12 
Hyades 6 18 + 7°5 44°0 +40 -18 — 2 
¢ Her 5 48 — 16°5 74°5 +54 -45 —26 
e Ind 6 44 + 6°7 88-5 +79 —-39 + 6 
61 Cyg 6 28 + 03 108°0 +92 -—-53 — 6], 


where the values of (A, D) represent the convergent points of the apparent 
motions and V, is the total space velocity of each group with respect to the 
Sun (the vector space velocities are in the direction /= 148°, b=0°; 1=58°, 
b=0°; and b=+490° for U, V, and W respectively). A search through the 
same catalogue for stars that might be.‘‘ members ”’ of synthetic groups with 
convergent point at 5= —go0° and V,=74'5, 88-5, and 108 km/sec indicated 
that only 2 or 3 chance members are expected to be included in each group. 

The colour-luminosity array for the stars in all three high-velocity groups 
is very similar to that for the galactic cluster M 67. 





In 1911 Benjamin Boss (1912@) announced the discovery that a dozen stars 
or so of large proper motion formed a moving cluster with 61 Cygni. Further 
investigations by Boss (19125) and by Russell (1912) appeared to establish 
the reality of the cluster but several years later when Rasmuson (1921!) 
re-examined the evidence with the new material then available, he concluded 
that ‘‘there are great reasons to exclude from the list of moving clusters the 
61 Cygni-stream...’’. Still later, Chaudhuri (1940) came to thesame conclusion. 

The purpose of the present investigation is not to re-hash the available data 
for the stars originally assigned to the 61 Cygni cluster but rather to search for 
stars sharing the well-determined space motions of selected high-velocity stars, 
including 61 Cygni. The stars chosen, for the accuracy of their proper motions, 
radial velocities and parallaxes, are as follows: 


Star 7 Po pe U V W 
a (km/sec) a 
B Hyi 0'153 +22°9 2°255 +53 — 43 — 28 
¢ Her o'1o1 —69'9 0608 +54 —47 — 24 
61 Cyg 0296 — 64:0 5°205 +92 =S3 sl 
e Ind 0°285 —40°4 4°692 +79 =—39 + 6 
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where the vector space velocities (galactic), U, V and W, are directed toward 
1=148°, b=0°; 1=58°, b=0°; and b=+90°, respectively. The catalogue 
searched consisted of 700 high-velocity stars, mainly from the lists published 
by Miczaika (1940) and by Roman (1955) but with a few additional stars 
obtained from recently published lists of radial velocities. 

Since most of the trigonometric parallaxes for the high-velocity stars are 
much less accurately determined than those for the stars listed above, the method 
used for selecting group members was that described in a previous paper (Eggen 
1958; Paper 1). Stars were admitted to candidacy for membership in a group 
if the observed position angle of the proper motion, 9, agreed with that computed 
from the group motion to the extent that A@sinA<10°, where A is the angular 
distance of the star from the convergent point (A, D) of the group motion. The 
candidates were then admitted to group membership only if the difference, 
Po—pe, between the observed radial velocity, po, and that computed from the 
group motion, pe, was less than, or equal to, 5, 6, or 7 km/sec, respectively, for 
radial velocities of quality a, b or c, as defined by R. E. Wilson (1953). 

The first two stars listed above, 8 Hydi and ¢ Herculis, obviously have a 
common space motion and the convergent point of what will be called the 
¢ Herculis Group was obtained from the mutual convergent point of the proper 
motions of these two stars. The values of (A, D) for what will be called the 
61 Cygni and « Indi Groups are merely the apices of the motion of the defining 
stars, computed from the radial velocities, proper motions and parallaxes listed 
above. The three convergent points and total space velocities with respect 
to the Sun, V;, are as follows (those for the Hyades and Sirius Groups, discussed 
in Paper I, are added for comparison) : 


Group A D V, 
h m °) km/sec 
Sirius 20 44 —42°7 18°4 
Hyades 6 18 + 7°5 44°0 
¢ Her 5 48 — 16°5 74°5 
e Ind 6 44 + 6°7 88-5 
61 Cyg 6 28 + 03 108°0 


Eight per cent, or 54, of the 700 stars examined were found to be members 
of the three high-velocity groups as follows: ¢ Herculis, 23 members; « Indi, 
15 members; and 61 Cygni, 16 members. The members of the ¢ Herculis, 
e Indi and 61 Cygni Groups are listed in Tables I, II and III, respectively, 
together with the observed position angle of the proper motion, 0), and the 
normalized residual from the computed value, A@sinA; the observed radial 
velocity, po, and its residual, Ap, and quality Q; the total annual proper motion p; 
the group parallax zg, computed from 4-738u/V,sinA; the observed magnitude 
and colour V,, and (P—V),,; and the spectral type and visual absolute magnitude, 
M,. The proper motions are mainly uncorrected values from the General 
Catalogue (Boss 1937) except for a few stars from the Yale zones or unpublished 
Greenwich determinations. The radial velocities and their quality have been 
taken from Wilson’s (1953) catalogue. ‘The magnitudes and colours were 
observed on the (P,V), system or reduced to that system from observations 
by Roman (1955) or by Stoy and his collaborators (cf. Evans, Menzies and Stoy 
1957). ‘The spectral types are mainly by Roman (1955). 
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As a test of the possibility that the group members might be so classified by 
chance coincidences of 6, and po with 0. and pe, the same catalogue of 700 stars 
was examined for ‘‘members’’ of a synthetic group with convergent point 
5=-—go0°. The galactic coordinates of this convergent, 1=270° and b= — 28°, 
place it near the centre of that hemisphere of the sky toward which the velocity 
vectors of the high-velocity stars show a strong preference (Oort 1926). 
Eleven per cent, or 76 of the 700 stars examined for membership in the synthetic 
group gave values of A@sinA<10°. When the criterion that Ap< 5, 6, or 7km/sec 
for stars of quality a, b, and c, respectively, was applied, the following ‘‘ members ’”’ 
were obtained for three values of V,: 

V, (km/sec) = 74°5 88-5 108 

No. = 3 (1) 3 
The one ‘‘member’’ for V,= 88-5 km/sec, shown in parentheses, is a star with 
a small radial velocity of quality c that also satisfies the requirements for the 
group with V,=74:5 km/sec. Apparently, we might expect two or three 
accidental members to the groups in Tables I-III but it does not seem possible 
to dismiss all the members as spurious. 

Because of the proximity of the apparent convergent points of the three 
high-velocity groups, there is some ambiguity in assigning four stars to a 
particular group. These stars, together with the appropriate values of A@ sin), 
Ap (in km/sec) and zg, are as follows: 


‘ 


¢ Herculis 61 Cygni e Indi 
A@sinA Ap Ty A@sinA Ap TT, A@sinA Ap Ty 
n Ret —2 —3 0'016 —8 +5 0009 : 
HD 66171 ° +6 0°023 +3 +2 0°030 
HD102365 +7 © 0°070 ° +4 0°0855 
HD104988 +6 +2 0°023 —I —I 0°028 


The assignment of these stars to the groups in Tables I-III has been guided by 
the trigonometric values of the parallaxes. 

Also, there are two high-velocity white dwarfs that may belong to either the 
61 Cygni or « Indi Groups, but because of (1) the uncertainty in the observed 
radial velocities of these faint stars with poor spectral lines and (2) the possibility 
that the observed line-shifts do not represent the true radial velocity, it is not 
possible to assign them definitely to either group. These stars are as follows: 


61 Cygni e Indi 
7 (p.e.) A@ sin A Ap Ty A@ sin A Ap Tg 
— 32°5613 0-108 0°007 +8 (— 18) "104 +1 (—1) 0°122 
L1534-1 ©°059 0°005 +3 (— 6) 0068 +6 (+4) o*100 


The mean trigonometric parallaxes, 7, indicate that both stars may belong to 
the 61 Cygni Group. 

Since the observed parallaxes of the stars, except 61 Cygni, « Indi, and the 
four stars mentioned above, were not used in assigning the members to the 
groups, a comparison between these and the group parallaxes offers an additional 
test of the reality of the groups. The parallaxes are listed in Table IV, where the 
absolute trigonometric values, 7, were taken from the Yale Parallax Catalogue 
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by TaBLe I 
irs 
int Members of the € Herculis Group 
,, 
ity Star 6 A@sinA po Ap QQ up ig Ve (P-V)g Sp. M, 
A ) ° ° (km/ sec) a a“ m m m 
: TV Psc 80 +3 +6 +1 b o-115 0°0075 4°6v +1°6: M3III -—r1-o: 
tic B Hyi 82 —3 +23 © @ 2'255 O'I510 2°75 +0°52 GilIV + 3°64 
eC HDg166 1117 —2 —14 —-6 b o-1r5 0°0075 6°76 +1715 K3 III +1°14 
ag p Per 129 —6 +28 o b 0172 o-0120 3°4v +1°53 M4III —1:-2 
7 Ret 26 —2 +45 —3 a o189 o-0160 5:26 +0°85 Ko +1:28 
HD30455 153 +3 +55 —3 a 0°435 070445 6°97 +0°51 G2V + 5°22 
HD43899 = 5536 +5 +70 +1 d o-0go o-0155 5:6: — Ko +1°6: 
HD68788 212 —6 —-3 +1 b 0567 0:0360 8:38 +0°74 Ki V +6°16 
| HD71377 257 -I +65 +6 c 0°094 o-0095 5°54 +111 K2III +042 
th HD89668 254 +1 +35 +7 Cc 0°670 0°0465 9°41 +1°02 K3V +7°75 
he HDgo2z50 237 —-7 +11 +3 b 0129 0:0085 6:49 +1°01 Ki III +1°14 
ee HD 98824 259 +7 + 5 +3 b orr1g 0:0075 7:04 +0°96 Ki III +1-42 
le HD100733 262 ° +18 © a 0088 o-0055 5°4v +1°60 M3III —o-g: 
HD106364 262 +3 —10 +6 b o113 0°0075 6°87 +1°05 K2III +1:°25 
8:77 +049 FoV +3°15 
” HD119425 255 +2 —42 -6 b 0127 0:0085 6:26 +1°:22 K3III +090 
ad ¢’ Lup 224 +3 —29 +4 a 0°134 070095 3°52 +1°48 gMo — 1°60 
A, HD150275 240 ° —32 +3 b 0289 o-0210 6:34 +091 sg Ko + 2°95 
¢ Her 309 ° —-70 —I a 0608 o-0990 2°80 +0°54 GoIV  +2°78 
€ Oct 118 +2 +12 —4 a 0:'084 0°0055 4°8: +1°6: gM6 —1°5: 
HD209134 82 -9 —-18 +6 b 2-092 o-1410 5°56 +0°92 K3V +6°51 
HD219829* 103 -3 —14 —-2 b 0484 0:0310 8:00 +0°72 KoV + 5°46 
HD221354 84 —9 —25 —4 b 1°085 0:0720 6°75 +0°74 KoV + 6:04 
, *ADS 16645, Am=1™, P=108 yrs: not plotted in Fig. 1. ‘The dynamical parallax is 
55 0”°033. 
) 
by TABLE II 
h Members of the « Indi Group 
€ 
ed Star 6 AOsinA po Ap Q. sz Tg Ve (P-V)g Sp. M, 
ity F ° > (km/sec) a P m m m 
at + 66°34 97 ° +10 -—-6 c 1°76 0096 10°34 +1°30 dM3 + 10°25 
A Aur 141 —I +66 -—3 a 0845 0:0720 4°68 +o0:50 GoV- + 3:97 
S: HD39853 58 +9 +87 +6 b 0:'076 o-o10o 5°66 +1°46 Ks III + 0°66 
HD40409 14 +1 +25 -—4 a 0°557 0°0315 4°59 +0°98 sg K3 + 2:09 
HD72324 231 —4 +75 —-1 b 0-085 o-00go 6:36 +0°93 GoIII + 1°13 
HDo9196 = 263 -9 +38 +6 b 0107 o-0060 5°79 +1730 K4III — 0°32 
HD102365 274 ° +15 +4 a 1°585 0°0855 4°87 +0°54 G5 V+ 4°53 
2 HD104988 276 -I +14 -1 b o-521 0:0285 8:19 +066 G8V + 5:47 
0 HD112758 283 +7 — 4 +2 b 0°839 0:0450 7°54 +0°70 dKr + 5°81 
- HD115539 272 -7 — 8 -1 b 0:096 o-0050 7:21 +0°87 G8 III- + 0°70 
IV 
HD135101 296 +4 —39 +6 b 0°656 00410 6°68 +0°58 G5 V + 4°74 
he 7°53 +0°63 dG6 + 5°59 
he HD137704 294 —6 —48 -6 b o-118 0°0070 5°46 +1°33 K4III — 0°30 
al HD142574 295 —2 —61 -4 b orogr o-0065 5°44 +1°52 K4III — o-50 
‘ HD148704 231 —I —59 +4 b 0532 o-0405 7:28 +0°65 G5 + 5°32 
€ € Ind 123 ° — 40 © a 4°692 0:2850 4°72 +0°95 K5V_ + 6°96 
ue 











Star 


HD10145 
HD18702 
HD23183 
HD32023 
HD35783 
HD39091 
B Col 

HD 40460 
HDs3501 
HDs55526 
HD66171 
HD108076 
HD112943 
HD120467 
HD182572 
61 Cyg 


oP 
109 
106 
115 


76: 


165 
15 
7 
183 
351 
352 
204 
272 
a7Y 
254 


52 


Olin F. Eggen 


TaBLe III 


Members of the 61 Cygni Group 


Aé sin A 


Po Ap Q pz 7, 
. (km/sec) i . 
+3 + 16 +3 b 0°734 0°0325 
+9 + 66 +1 b 0°677 0°0375 
+3 + 78 +2 c 0°130 0:0080 
+6 +106 +6 c 0:068 0-0080 
o + 19 -—3 b 0278 0°0125 
+3 + 12 -—5 Db I°IIO 0°0495 
—-6 + 89 +3 a 0402 0°0290 
+9 + 96 +1 Cc 0°'078 0°0075 
o + 39:«6—-1 «Cb 0°239 O-o1l5 
+5 + 64 -—5 a 0°204 O'OII5 
o + 37 +6 b 0503 0°0230 
+2 —- 1 -—4 b 0'579 0:0260 
+1 — 12 +2 b 0768 0:0340 
—8 35 —1 b 1°814 0°0840 
© —100 —3 a 0960 0:0890 
° — 64 © a 5°205 0°2960 
TaBLe IV 


Vy 


m 
7°66 
8-12 
6°14 
g*I0 
7°69 
5°63 
3°06 
6°59 
48: 
5°15 
8-20 
8-02 
9°79 
8-3: 
5°14 
5°19 
6:02 
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(P-V) xz Sp. M, 
m m 
+0°58 G5V +5:22 
+076 KoV +6:00 
+o-92 KolIII +0°66 
+047 F8V + 3°62 
+035 F6V 3°57 
+0°46 Go +410 
+110 gKr +0°37 
+o:93 Ki III +0°97 
+1°5: Mo +o°1: 
+1714 K4 +0°45 
+o'51 G2V +5°01 
+046 GoV +5°10 
+1:06 dMo +7°45 
— dK6 +7°9: 
+068 G8IV  +4:°89 
+1:08 K5V +7°55 
+1':24 K7V +8-38 


Comparison of the trigonometric and group parallaxes for the stars in Tables I-III 
(Unit = 0”-oo1) 


Star 

TV Psc 

B Hyi 

+ 66° 34 
409166 
HDro0145 
HD18702 
p Per 

n Ret 
HD30455 
A Aur 
HD35783 
HD 39091 
B Col 
HD40409 
HDs53501 
HDs55526 
HD66171 
HD68788 
HD72324 
HD89668 


ag 


8 


151 


96 

8 
32 
38 
12 
16 
44 
a2 
12 
50 
29 
32 
12 
12 
23 
36 


9 
46 


%¢ 
13A 
143Y, 166C 
102M, 80Yh, 94G 
-—3G 
32M, 45G 
38M, 23Y, —10S, 41D 
2A, 30M 
7Y 
47M, 15V 
65A, 65M, 68S 
12G 
38C ' 
sY, 45C 
60Y, 8C 
15Y, 46C 
—1Y, 52C 
18M, 13G 
10M, 14G, 22D 
23W 
40M, 42Y,—10V, —4C 


Star 


HD102365 
HD104988 
HD108076 
HD112758 
HD112943 
HD119425 
HD120467 
HD129245 
HD135101 
¢’ Lup 
HD142574 
HD148704 
HD150275 
¢ Her 
HD182572 
61 Cyg 

e Ind 
HD219134 
HD219829 
HD221354 


89 
296 
285 
141 

31 

72 


94Y, 89C 

—10M, 15C 

17M 

79M, 35C 

26M, soY, 41V, 71C 

8A, —5sY, —12W 

83Y, 56C 

12G 

28A, 7M 

6Y 

27M, 38W 

54Y, 75C 

19A, 9G 

95A, 108M, 104Yk, 1248 
51A, 75M, 101 Yk 

296 mean 

285 mean 

140A, 159M, 181S 

6M, 44Y, —17k, 8S, 52D 
22M 


(Jenkins 1952) and the group parallaxes, 7g, from Tables I-III; the abbreviations 
used to designate the source of the individual trigonometric values are those 
given by Jenkins. Because of the large scatter in the trigonometric values, there 
is probably no star in Table IV that could be eliminated as a group member on 
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52 
17 
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37 
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89 
55 
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the basis of this comparison. It is noteworthy that of the 600 stars within 
20 parsecs of the Sun for which radial velocities have been determined (Gliese 
1957) eight per cent belong to the Hyades, ¢ Herculis, e Indi or 61 Cygni Groups. 

The colour-luminosity array for the stars in Tables I-III is shown in Fig. 1, 
where the ¢ Herculis, « Indi and 61 Cygni Group members are indicated by 
crosses, open circles and filled circles, respectively. ‘The standard main sequence 
(Eggen 1955) is shown with a half-width of o™-2. This colour—luminosity array 
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° ON)e 
Fic. 1.—Colour-luminosity array for stars in the € Herculis (crosses), 
e Indi (open circles), and 61 Cygni (filled circles) Groups. 
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is almost identical with that obtained from the members of the cluster M 67 
(Johnson and Sandage 1955) and, indeed, nearly every group star in Fig. 1 
can be matched in colour and in luminosity by a member of the M 67 cluster 
except for the five bright stars near M,~ —1™ and (P—V),~ +1™5 to +1™-6, 
which occur in the groups‘but not in the cluster. These M-type giants, members 
of the £ Herculis Group, are mainly known (irregular) or suspected light variables. 

In the notation of a previous paper (Woolley and Eggen 1958) the galactic 
orbits of the members of these high-velocity groups are of class D, that is, the 
pericentric distances of the stars are between two-thirds and one-half of their 
present distances from the galactic centre. 

I am greatly indebted to Mr J. Alexander, Miss L. Mather and Miss A. F.. 
Heather for much computational assistance. 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex: 
1957 December 30. 
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THE VELOCITY ELLIPSOID OF Ao STARS 
Jj. B. Alexander 
(Communicated by the Astronomer Royal) 
n of (Received 1957 December 28) 


Summary 


A kinematic average absolute magnitude of +0™-14 is derived for stars of 
spectral type Ao by comparing the line-of-sight velocities near the solar apex 
with the transverse velocities in a region of the sky 90° away. Assuming that 
every Ao star has the same absolute magnitude, the approximate distances 
can be found from the observed apparent magnitudes, and this method applied 
to 475 stars of this type brighter than 6™-5 is used to obtain the equatorial 
linear components of velocity from the proper motions and the radial velo- 
cities. The solar motion and the constants of the velocity ellipsoid are 
| evaluated, and the results are compared with those of other investigations. 
The galactic longitude of the major axis is found to be 354°, and the root- 
mean-squares of the peculiar velocities in the directions of the principal axes 
are 18-9 km/sec, 11-1 km/sec and 8-4 km/sec respectively. 





1. Introduction.—In any study of stellar motions it is always desirable to 
know the actual linear velocities of the stars concerned, and this requires a 
knowledge of the distances involved. This is only possible, in general, for 
nearby stars with well-determined trigonometric parallaxes. Reasonably accurate 
distances can be found from the apparent magnitudes if we are dealing with a 
sample of stars which has a small dispersion in absolute magnitude about a known 
mean value. The stars of type A would seem to satisfy this condition. Most 
previous investigations on the motions of this spectral class have been made 
using either proper motions alone or radial velocities alone since there are too 
few A stars with accurate trigonometric parallaxes to form a significant sample. 

In making any statistical study of the movements of a group of stars, it is 
advantageous to use as homogeneous a set of observational material as possible. 
In Wilson’s General Catalogue of Stellar Radial Velocities (1), from which the 
stars in this treatment are selected, we now have available a comprehensive 
collection of data reduced to a common system which is probably free from serious 
systematic error. Only Ao stars brighter than 6™-5 are included in the present 
investigation. ‘Types cAo and Aop are excluded as there is evidence to show 
that these have absolute magnitudes systematically brighter than those of the . 
normal Ao class. Known visual binaries are also not considered. In the interests 
of homogeneity, the radial velocities and spectral classification are taken exclusively 
from Wilson. 

The proper motions and apparent visual magnitudes used are those given 
| in the Boss G.C. ‘There may exist systematic errors in these proper motions 
which have a marked effect on the derived distributions of the linear components 
of velocity, but there is no simple way of allowing for these. On the other hand, 
corrections for accidental errors in the proper motions can easily be made. 


1} 
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‘There are some Ao stars brighter than 6™-5 in the Southern Hemisphere 
which do not appear in Wilson, but this bias in the sampling probably has no 
serious effect on the final results. 

2. Previous investigations.—In the various sections, reference is made to the 
papers listed below. In making any comparisons, due attention must be paid 
to differences in the spectral ranges of the samples. 


Author Material used 
Nordstrém (2) Radial velocities of 697 A-type stars brighter 
than 6™-o, 
Smart (3) Radial speeds of 643 A stars in Schlesinger’s 
Catalogue of Bright Stars. 
MacRae and Nevin (4) Radial velocities of 457 B8-A3 stars of average 


photovisual magnitude 6™-9. 
Wilson and Raymond (5) Proper motions of 4457 G.C. stars of spectral 
type B7—A4 and brighter than 7™-05. 
Delhaye (6) Proper motions of 5089 stars of spectral type 
B8-As5 and magnitudes between 6™-0 and 7-5, 


3. Absolute magnitude.—With the usual notation, the equatorial linear 
components of velocity of a star with coordinates («,5) are given by the 
relations below 

u= —(xp,/p)sinacosd— (ku,/p)cosasind+Rcosacosd | 
v= (kp,/p)cosacosd—(xpu,/p)sinasind+Rsinacosd +}. (1) 
w= (xp,/p)cosd + Rsind 
If p is the parallax in seconds, R the radial velocity in km/sec and pw, and p, are 
measured in seconds of arc per annum, then the value of « is about 4°74. 

The parallaxes can be found from the apparent magnitudes only if the absolute 
magnitudes are known. It is assumed that there exists some average absolute 
magnitude M, which when assigned to every Ao star yields reasonable parallaxes 
for statistical purposes. 

To determine Mp, a trial value of +o0™-5 is taken. With this figure, the 
values of u, v and w can be calculated from equations (1) for any star. Two 
regions are considered. 

For 54 stars in an area surrounding the solar apex (the exact direction of the 
solar motion is not critical for this purpose ), the arithmetic means of the equatorial 
linear components are as follows 

i= — 2:6 km/sec; d= +157 km/sec; w= — 8-2 km/sec. 
These systematic velocities caused by the solar motion are in this case dependent 
mostly on the radial velocity data with little weight attached to the proper motions. 
Hence an exact choice of M, is not important in obtaining @, @ and @ for this 


group of stars. 
A similar procedure is performed for 54 stars lying in a belt go° away from 


the apex. The corresponding values are 

ti=0°0 km/sec; d= +15'1 km/sec; w= —5°3 km/sec. 
Here, the means of the velocities depend critically on the choice of M,. By 
adjusting M, such that a? + ¢*+ is the same for both regions, the average 
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absolute magnitude of the Ao stars in Wilson is obtained. After allowance is 
made for the small effects of galactic rotation on the systematic motions of stars 
in a given region, it is found that 

My= +0™14+0™25. 

The above value of M, is used throughout this paper. ‘The following estimates 
of absolute magnitude as a function of spectral type and luminosity class are 
taken from a table obtained by Keenan and Morgan (7). ‘The results are derived 
from a study of cluster motions and trigonometric parallaxes. Considering that 


TABLE | 

Type Main Seq. Subgiants Giants Supergiants 

V IV III II Ib Ia 
B8 —0'5 “59 “fo ~“@% “og 7s 
Bo 0'o =O —2°0 —3°3 —5°5 —7°O 
Ao +0°3 —O°4 —1I —3°0 —4'8 —7°0 
A2 +1°2 +0:2 =e'7 =—27 =£9 —7°O 
A3 +1°8 +1°0 —03 —2°5 —4°6 —7°'0 


the majority of Ao stars are placed in luminosity class V, there is no contradiction 
with the value of M, found kinematically. 

4. Solar motion.—With M, taken to be +0™-14, the values of u, v and w are 
evaluated from equations (1) for each of the 475 stars. The components U,, V, 
and W, of the solar motion with respect to a group of N stars are determined 
from the relations 

NU, = — [u]; NVy= — [2]; NW, = — [x]. 
(Throughout this paper, square brackets denote summation with respect to the 
appropriate group of stars.) ‘The values of U,, V, and W, for the 475 stars are 
given in row A of Tablé II below. If 117 stars with velocities greater than 
35 km/sec with respect to the standard of rest defined by A are excluded, then 
we obtain the results B for the remaining 358 stars. Finally, if instead of using 
arithmetic means to find U,, V, and W, we take median counts of the 475 stars, 
then we derive the values C. The solar motion is taken to be that defined by B. 


TABLE II 
U, (km/sec) V,, (km/sec) W, (km/sec) 
A +0°8 — 16°2 +9°2 
B +1°1 —16°5 +83 
Cc +0°6 —17°5 +81 


These components correspond to a velocity of 18-5 + 1-6 km/sec towards the 
following apex (1950): 
Gig = 273 °°8 + 2°°2; 59 = +26°°6 4+ 2°0. 
With respect to galactic coordinates the corresponding values are 
[=21°-4; b= +17°:5 (Lund pole). 
The galactic linear components of the solar motion are 
X= + 16-4 km/sec ; Y,= +6°4 km/sec; Zy= +5°6 km/sec. 


tt” 
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The zero of galactic longitude is defined in terms of essentially local characteristics. 
For general purposes, it is more convenient to take axes determined by the 
position of the centre of the galaxy in galactic longitude 328°. Consider two 
axes € and 7 in the directions (/= 148°, b=0°) and (/=58°, b=0°) respectively. 
‘Then the components of the solar motion in these directions are 


&)=—10'5km/sec; jp = + 14:2 km/sec. 


For a group of stars randomly distributed over the celestial sphere, the first- 
order terms of galactic rotation have no effect on the derived constants of the solar 
motion. ‘The sample of Ao stars, besides exhibiting a marked concentration 
towards the plane of the galaxy, tends to show a deficiency of stars for an area 
in the Southern Hemisphere compared with the diametrically opposite northern 
region. This latter irregularity has some effect on the systematic means of the 
velocities, but as a rough estimate puts the resultant errors in Xy and Y, at less 
than 0-2 km/sec, no attempt is made to allow for this in obtaining the components 
of the solar motion. As with previous authors, galactic rotation is again ignored 
when the constants of the velocity ellipsoid are obtained. 

From the radial velocities of 874 stars of which about 54 per cent are of type A, 
MacRae and Nevin obtain an apex at (270°, + 36°). With this direction assumed, 
the magnitude of the solar motion with respect to the A stars alone is found to 
be 18:37 +1-11 km/sec. Nordstrom obtains a value of 16-3 + 0-9 km/sec towards 
(263°°8, +22°:2) for his A stars, although the corresponding velocity for the 
B8—Bg group is as high as 21-2 km/sec. 

From proper motions alone, Wilson and Raymond place the apex at 
(270°:2+1°°3, +26°-0+1°2). 

5. Axes of velocity ellipsoid—The existence of such groups as the ‘Taurus 
and Ursa Major clusters shows that the distribution of the peculiar velocities 
of the Ao stars in the neighbourhood of the Sun cannot be represented by ary 
simple mathematical expression. Suppose, however, adopting Schwarzschild’s 
ellipsoidal hypothesis, that the velocities of our sample can be represented 
reasonably well by a frequency function of the form 


exp{— (aU? +bV?+cW?+2fVW+2gUW+2hUV)}. 


This is later seen to be a plausible assumption when it is found in Section 7 that 
the histograms of the velocities in the three directions do not depart appreciably 
from the corresponding Gaussian forms. With an ellipsoidal distribution, 
steady-state theory applied to the rotating galaxy demands that the axis of greatest 
mobility for components in the galactic plane is directed towards the galactic 
centre. Previous work indicates a significant deviation of the observed axes of 
the velocity ellipsoid of A stars from their ‘‘equilibrium’’ position. 

In this section, let (U, V, W) be the peculiar components of equatorial linear 
velocity of a star found by adding the values (Uy, Vo, Wo) derived in Section 4 
to the corresponding components (u, v, w) relative to the Sun evaluated from 
equations (1). Consider the direction whose direction-cosines with respect 
to equatorial axes are (/,m,n). The the sum S of the squares of the components 
of peculiar velocity in this direction is [(/U +mV +nW)?], where as before the 
square brackets denote summation. Using Lagrange’s Method of Undetermined 
Multipliers, we obtain the following equations for the directions in which S has 
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stationary values. 
({U?]—A)l+ [UV]m+ [UW]n=0 | 
[UV]l+([V?]—A)m+ [VW]n=0 
[UW]l+ [VW]m+ ([W?]—A)n=0 
P+m+n=1 | 
(A is the undetermined multiplier). 
Denote by A the determinant 
([U?]—A)_ [UV] [UW] 
[UV] ((V*]-A) [VW] 
[UW] [VW] ({W}-a) 
Then each of the three roots of the cubic in A, A=o, when in turn substituted 
in equations (2), gives a set of values (/, m,n) corresponding to a direction in 
which S has a stationary value. The situation is analogous to that of the inertia 
tensor in rigid-body dynamics. Corresponding to the inertia quadric we have 
in this case the quadric 
[U2 ]x? + [V2] y? + [W?]s? + 2[{VW] yz + 2[UW]xz + 2[UV ]xy =1, 
and by the above method we are in fact locating the principal axes of this. 
The process therefore yields three mutually perpendicular directions whatever 
the distribution of the peculiar velocities may in fact be. In particular, if we are 
considering the ellipsoidal form 
exp {—(aU?+bV*+cW?+2fVW+ 2gUW+2hUV); 
then the stationary values of S occur in the directions of the principal axes of 
the velocity ellipsoid. 

The equatorial (1950) and galactic coordinates of the directions found by the 
above method applied to the 475 stars are given in the first three rows of Table III 
below. If 131 stars with peculiar speeds greater than 30 km/sec are ignored, 
the axis of greatest mobility is moved about 12° further away from the direction 
of the galactic centre. The results for the remaining 344 stars are given in the 
last three rows. No attempt is made to identify individual stars as members 
of clusters, but the above restriction should remove, for instance, all likely ‘Taurus 
members and about half of the Ursa Major group present in the original sample. 





Tas_e III 
a (1950) 8 (1950) l b 
x 279°8 ~ OF 353°6 - o4 
y 358°3 +59°5 $3°8 — 22 
z 193°6 +29°6 30°0 + 86°38 
x 286°5 + 33 5°7 — 36 
y 22°4 +61°0 95°7 —-o”7 
2’ 194°6 + 28°7 17°! + 86°3 


The galactic longitudes of the major axis of the velocity ellipsoid according to 
other authors are as follows: 


Nordstrém 357°°91 7°33 
Smart 14° 

MacRae and Nevin (First Method) 342°°3+10°7 
MacRae and Nevin (Second Method) 953° +145 
Delhaye 35t°2t 2°°3 
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6. Magnitude dispersion.—In Section 5, it is assumed that the absolute 
magnitude of an Ao star is +0™-14. Errors in the directions of the axes of the 
ellipsoid found by this analysis result from the fact that the true average absolute 
magnitude may differ considerably from the above figure and that the individual 
stellar luminosities have a spread about their mean. Interstellar absorption in 
the visual wave-lengths will cause further discrepancies, although the error in 
the derived value of M, due to this tends to reduce the effect on a given star. 
It will be assumed that the frequency function f(M) of the true absolute magnitudes 
is independent of the position and motion of a star. Firstly, consider motion in 
the U-component. Suppose we have a group of stars so situated that U is derived 
from proper motions alone (i.e. the stars are in a belt go° away from the U-axis). 
Let the true distance of a star assumed to be at distance d’ be d, where d’ = dt. 
‘Then ¢ has a certain frequency function g(t) which is independent of the position 
and motion of a star, since g(t) is dependent solely on the form of f(M). For 
a star with velocity U, the motion relative to the Sun is U— U,, where U, is the 
component of the solar motion. We interpret this motion to be ¢(U — U)) relative 
tothe Sun. If we use U,, +6U, for the value of the component of the solar motion, 
we write down t(U — U,) + U, + 6U, for the star’s component of peculiar velocity. 

Now the directions of the axes of the ellipsoid found in Section 5 are a function 
of the quantities [U*], [V7], [W?], [VW], [UW] and [UV]. In fact, only the 
ratios of these six numbers are significant for our purpose. Consider the term 
[U?]. We form [U?]’ given by the equation 

[U?]’ = [{t(U — Up) + Uy + 8U,}?] = [{tU + (9 —t)U, + 8 U,}?]. 
Irom the definition of the solar motion, we have 
[U]=o. 
‘There is good evidence from comparison with other results that 5U, is small, 
and this term is neglected. We thus have the relation 


vt ina tall il Ns deed 
i.e. [Uy = an g(t)? dt[ U2] + [. “ g(t)(1-2)2 atU,2. (3) 


We are, in practice, concerned with stars distributed over the whole sky. We 
shall suppose that equation (3) still holds if g(t) is replaced by g’(t), a function 
which allows for the fact that the part of the U-component derived from radial 
velocities is unaffected by magnitude dispersion. The function g’(t) depends 
on the distribution of the sample of stars relative to the U-axis. We shall assume 
that equations similar to (3) apply to the V and W components with the same 
function g’(t) inserted in place of g(t). Define k and / by the two relations 


. af atid b haaim 


0 
ay g' (t)t? dt. 
We then have the following equations for a group of N stars: 
[U?]=/[U?]' —kNU2 
[V?]=/[V?]' —RNV,? 
[W?]=/[W?]’ —kNW,?. 
Similarly, we obtain 
[VW]=l[VW]' —kNV,W, 
[UW]=/[UW]’-—kRNU,W, 
[UV]=I[UV]'—RNUV,. | 
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Hence, for a given set of [U?]’, [V?]’, [W?]’, [VW]’, [UW] and [UV]’, we can 
correct for the effects of magnitude dispersion on the directions of the axes if 
we know the value of k//, where 


x |. g’(t)(1 —t)* dt. 


If the absolute magnitudes follow a distribution of the form exp {—A|M—M,|} 
about some value 7,, where M, — M,=5 log,, 8, then we have 
g(t)=3Not*"B-* — (t< B) 
g(t)= Nat» (t>B) 


(where z= 5hlog,,¢=2°1715A). ‘The values of 


I "Hh 
— g —t)dt 
Ni], g(t) ) 


for various % are given in ‘Table IV. ‘he last two columns correspond to the 
two values of M, for which |M,—M,|=o"™-25. 


ki/l= 


"TABLE IV 


v p=1 p=1-122 p=o'891 
10°O 0°02 0°04 0°03 

7°O "05 08 05 

5'0 a 16 els) 
4°0 *20 29 “16 

3°0 55 0°74 "42 

es 0'90 I*I9 0°69 

2°5 1°40 1°83 1:08 

2°3 2°63 3°40 2°06 


The corresponding values .of 

tl og (t)\(1—ajedt 

= —bya 

N |, }& 

will be somewhat less, the exact values depending on the distribution of the 
stars concerned. If the value of k// is 1, we find from equations (2) and (4) that 
the galactic coordinates of the axes of the ellipsoid for the 475 stars are 


l b 
348 11 ~ 33 
79 °°1 10 +2 
50°°7 +785 
If k/l is 1, the corresponding Jongitude of the major axis is 343. ‘The minor 


axis is tilted as much as 22}° away from the Galactic Pole, but it is unlikely that 
this is really the case. If we stipulate that the galactic latitude of the true minor 
axis is greater than 75° and that |1/7, — M,| is less than o™-25, we must have that « 
is greater than 3. This value of x corresponds to a standard deviation in the 
absolute magnitudes about M, of about 1™-o0. 

The above analysis is rough, but it is expected that it simulates the effects 
of magnitude dispersion to a large extent. 
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7. Velocity distributions.—We shall take axes Oxyz defined by the principal 
directions found in Section 5, viz. 


x (1950) 8 (1950) 
x 279°°8 — 6%7 
y 358°3 +59°°5 
2 193°°6 +29°°6 


‘The peculiar components of velocity in these directions are computed for each 
star. The velocity distributions taken without regard to sign are given in 
‘Tables V, VI and VII. In each case, the Gaussian distribution with the same 
median velocity is quoted for comparison. 


"TABLE V 
X velocity Number of Gaussian 
(km/sec) Ao stars distribution 

o-5 1143 97 

5-10 87 gl 
10-15 615 794 
15-20 57 65 
20-25 40 50 
25-30 504 354 
30-35 244 24 
35-40 163 15 
40-45 11} 9 
45-50 5 5 
50-55 3 24 
55-60 I I 
60-65 2 4 
65+ I 4 

The median velocity is 13°0 km/sec. 
TABLE VI 
Y velocity Number of Gaussian 
(km/sec) Ao stars distribution 

o-5 160 157 

5-10 128 131 
10-15 87 gi 
15-20 41 534 
20-25 334 26 
25-30 144 10} 
30-35 4 34 
35-40 2 I 
40-45 5 ° 
45+ 44 ° 


‘The median velocity is 7-9 km/sec. 
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yal Tase VII 
Z velocity Number of Gaussian 
(km/sec) Ao stars distribution 
0-5 1994 1964 
5-10 1494 147 
10-15 68 83 
15-20 314 344 
r 20-25 16 II 
; 25-30 3) 24 
- 30-35 4 
ne 35-40 2% \ 
40-45 4 ; 
45+ ° 


‘The median velocity is 6-2 km/sec. 


In order to derive the constants of the velocity distributions, the median 
is to be preferred to other mean values since it is less affected by the presence 
of an excess of high-velocity stars. The medians and actual root-mean-squares 
of the velocities in the directions of the principal axes are given below. In the 
third row, the root-mean-square of a Gaussian distribution with the same median 
velocity as the actual distribution is in each case quoted. 


X Y Z 
Median (km/sec) 13°0 79 6:2 
Actual r.m.s. (km/sec) 20°3 13°4 10°I 
Gaussian r.m.s. (km/sec) 19°3 11°7 g:2 


Allowance can be made for accidental errors in the radial velocity and proper- 
motion data. The probable error of a radial velocity determination is taken to 
be 2:5 km/sec (8). The probable errors of the transverse velocities due to 
inaccuracies in the proper motions can be evaluated from the corresponding 
angular errors given in Boss. An average value of 2-5 km/sec is obtained for the 
sample of Ao stars taken. Correcting for these observational errors, we obtain 
the following median velocities : 

X Y Z 
Median (km/sec) 12°8 75 5°7 


The corresponding root-mean-squares 6,, 5, and 6, are derived on the 

assumption that the distributions are strictly Gaussian. In ‘Table VIII these 

values are compared with previous determinations. Smart assumes that the ; 
two minor axes are equal. 


TABLE VIII 


Author 3, 3, é, 6/8, 6,/6, 
km/sec km/sec km/sec 
Present 189 II‘l 8-4 0°59 0°44 
MacRae and Nevin 14°8 gl 9°4 ‘61 64 
Smart 19°9 11°8 - ‘61 - 
Nordstrém 7s 12°4 3°5 73 “50 
Delhaye - - - 0°63 O'sI 
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The values of 5,, 5, and 5, found by the present method depend on the 
choice of M,. Since the solar motion and the components of velocity of the 
stars relative to the Sun are derived from both proper motions and radial 
velocities, the effect of an error in M, on the distribution of the peculiar velocities 
is not straightforward. ‘The individual variations in absolute magnitude will 
also alter the derived components relative to the Sun, and the errors in the 
peculiar velocities will tend to be greater for directions in which the component 
of the solar motion is large. 

The root-mean-squares of the components in the directions é and » (defined 
in Section 4) are found to be 17-7 km/sec and 12-9 km/sec respectively. 

8. In Fig. 1, the projections of the peculiar velocities on the galactic plane 
are shown. ‘The horizontal axis corresponds to the direction of greatest mobility 
(in galactic longitude 354°) found in Section 5. ‘Ten of the Ao stars have 
velocities which fall outside the limits of the figure. 

+40 T 7 T T T T T T "7 





“hE 


o Y(KM/SEC.) & 











1 1 1 4 1 i 1 
— 40 —20 0 X (KM/SEC.) + 20 +40 


Fic. 1.—The projections of the peculiar velocities on the galactic plane. 





g. Conclusion.—The results for the solar motion are in reasonable agreement 
with previous work,, and this is additional confirmation that the value of M, 
used is near the true figure. The vertex deviation of the velocity ellipsoid found 
by previous authors is confirmed, and there is evidence to show that this is 
increased rather than decreased by the removal of the high-velocity members. 
The ratio of 5. to 5, is low compared with other determinations. This may be 
due to the fact that magnitude dispersion has less effect in the z-direction, for 
which the component of the solar motion is small, than in the x-direction. 

This work was carried out under the supervision of Dr R. v. d. R. Woolley, 
to whom the writer is greatly indebted for continual advice and discussion. 
Thanks are also due to Mrs V. A. J. Mitchell, who performed some of the 
calculations. 

Royal Greenwich Observatory, 


Herstmonceux Castle, Sussex: 
1957 October. 
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THE CLUSTER SYSTEM OF THE SMALL MAGELLANIC CLOUD 
E. M. Lindsay 


(Received 1958 January 3) 


Summary 


A cluster survey of the Small Magellanic Cloud has been completed 
using the ADH telescope and 30™ exposure 103a.O plates. Altogether 116 
objects are regarded as Cloud clusters. "These extend to the wing of the 
Cloud. Apparent magnitudes have been estimated on 10-inch Metcalf 
plates. ‘The cluster system gives dimensions of 8000 parsecs for the longest 
diagonal (including the wing) and 5400 parsecs for the smallest diagonal. 
‘There seems no evidence of spiral structure in the cluster distribution. 





Material for the Catalogue.—A previous paper (1), to be referred to as Paper I, 
contains a list of ninety-four objects which were considered to be clusters in the 
Small Magellanic Cloud. A further thirty-two objects, contained in ‘Table II 
of this paper, were regarded as possibly clusters. ‘The material was obtained 
from ADH plates of long exposure centred on the Cloud at R.A. o"35™, 
Dec. —72°-8 and one short exposure plate of 10™ centred at o!51™, —73°-0. 
Further material is now available for the outlying regions of the Cloud. This 
consists of ADH plates, types 103a.0 and 103a.E, of 30™ exposures. ‘The plates 
overlap so as to provide checks on doubtful objects. They cover an area around 
the Cloud such that clusters outside of this could not with certainty be regarded 
as Cloud members. On the eastern side the programme ties up with a simila: 
one, now nearly completed, on the Large Cloud. 

It is considered that Nos. 1, 3, 14, 28, 30, 31 and 32 of Table II, Paper I, 
should, in the light of the new material, be accepted as certainly clusters. All 
of these had been listed by the writer as probably clusters. ‘The last three were 
on the 10™ exposure plate referred to above. The last one is No. 5 of ‘Table I] 
of Jenka Mohr’s list (2) of outlying clusters of the Small Cloud found on 
long-exposure Bruce plates. Seven objects are contained in this Table. Six of 
these have been verified as clusters on the new ADH plates. No. 4 has not been 
verified ; in its position is an object which is almost certainly a galaxy. Miss Mohr 
had described this as ‘‘ hazy; somewhat doubtful as cluster’’ and it has been 
omitted from the present list. 

In addition to this total of twelve verified objects, ten new clusters have been 
found, most of which are distant from the main body of the Cloud. ‘Table Ia 
gives their number, approximate positions (1900), diameters to the nearest five 
parsecs and remarks. ‘Table Id gives the same data for Jenka Mohr’s outlying 
clusters, her No. 4 being excluded. The diameters were derived from estimated 
angular diameters on ADH plates, using de Vaucouleurs’ distance (3) of 
46 kiloparsecs. 

A note by de Vaucouleurs (4) on No. 7 (NGC 643), suggesting it was probably 
a Cloud member, appeared after the writer had included it in the above list but 
had not associated it with NGC 643. On a 30™ exposure ADH 103a.0 plate 
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it is a conspicuous object, clearly recognizable as a cluster with 20-30 faint 
resolved stars. The texture and distribution of the outlying stars suggest it is 
D globular; otherwise it would have been classified as certainly open. 
TABLE Ia 
New clusters in the Small Magellanic Cloud 
No. R.A. Dec. (1900) d Remarks 
(1900) (parsecs) 
h m E , 
I 00 184 —75 37°4 5 Open. 
2 00 43°4 69 283 5 Open; faint, poor; the most northerly cluster 
found. 
3 00 45°0 70 26:1 10 Open; faint, poor. 
4 OI 23:2 73 53°9 75 Open; nebulosity; stellar association with Nos. 1o1 
and 103. 
re or 28:8 76 34°7 5 Open; faint, poor. 
6 OI 31°5 74 40°4 5 Open. 
7 OI 33°5 76 04°4 25 Possibly globular, bright; stars resolved; NGC 643 
8 OI 34°5 75 585 5 Open; very faint, poor. 
9 o1 48°90 74 50°9 5 Open; very compact. 
I, 10 OF 55°5 78 08-6 5 Open; the most southerly and with No.7, Table I 4, 
the the most easterly cluster found. 
II 
ed TaBLe Ib 
;m, Miss Mohr’s outlying clusters in the Small Cloud (No. 4 omitted) 
*0. No. R.A. Dec. (1900) d Remarks 
his (1900) (parsecs) 
tes h m :# } ; 

d I 23 58:9 —74 O19 45 Circular, possibly globular. Many stars resolved. 
in The most westerly known cluster of the Cloud. 
led 2 OI OI1'! 70 52°9 5 Open; poor. 
ila: 3. -O1:-29°5 72 28:0, 10 Open. 

5 OF 32:2 73) 22°77 10 ©. Open. 
6 OI 47°9 74 13°3 30 ~=3=- Circular, possibly globular. Stars resolved. 
All 7  O1 §5°5 74 42°1 5 Open; very compact. 
— Table II of Paper I contained thirty-two objects which were considered to be 
I possible clusters; as mentioned above, seven of these have now been included 
= as certain, leaving twenty-five. On the new material there are an additional 
c of fourteen objects which resemble clusters, most of which are near the limit of the 
wen plates. Their positions are listed in Table II of the present paper. None of these 
ohr possible clusters has been included in the Catalogue of Table IIT. 
een 

TABLE II 

eon Additional possible clusters in the Small Magellanic Cloud 
: I a a R.A. (1900) Dec. (1900) No. R.A. (1900) Dec. (1900) 
five h m v ‘ h m . ‘i 
ying I 00 09'2 —70 36°0 8 OI 09°5 —75 42°9 
ated 2 oo 14°8 75 07°! 9 OI 09°9 71 29°6 
of 3 00 29°7 74 11°3 10 OI 20°5 75 31°6 

4 00 33°5 71 17°! II O1 3273 72 53°6 
bly 5 00 43°6 76 21°7 12 OI 39°3 71 39°9 
IDLY 6 OI 04°3 72 23°8 13 OI 40°2 71 46°9 
but 7 Ol 04°4 71 081 14 OI 40°9 73 51°1 
late 
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The Catalogue.—With the clusters of ‘Table I of Paper I, those of Table IT 
of the same paper regarded as certain, and the above objects of ‘Table La and I 4, 
116 objects are now considered to be clusters of the Small Cloud. To avoid 
confusion these have been collected in ‘Table II] and renumbered in order of 
right ascension (1900). ‘lhe right ascensions and declinations are approximate, 
being obtained with reference to stars of known positions plotted on squared 
paper. The fourth column gives the apparent photographic magnitudes as 
determined on MF (10-inch Metcalf refractor) plates. In some cases it was 
impossible to make any estimate of magnitude. The fifth column gives the 
diameters in minutes of arc. The last column gives a reference to published 
lists: I and II refer to Tables I and II of Paper 1; Ia and 1b to Tables Ia and Ib 
above; K to Kron’s list (§); the Shapley-Miss Wilson numbers (6) are those 
of HC 275 and 276; NGC or IC numbers are in brackets. 


TasBLe III 
Catalogue of clusters in the Small Magellanic Cloud 
No. R.A. (1900) Dec. (1900) My d’ Reference 
h m ne, 

I 23 58:9 —74 O1°9 II°7 4'1 Ibi 

2 00 08:2 74 02°7 16°0 4°5 Ilr 

3 oo 13°8 74 52°4 15°5 57 Ir 

4 oo 18-2 74. 183 13°9 13°6 Iz; K1; 276/80 

5 oo 18:4 75 374 15°7 4°5 lar 

6 oo 18°5 74 13°4 153 6:8 14; K4 

7 00 20°! 74 18-7 13'2 113 16; K5; 276/82 

8 00 20°3 74 28° 11°6 20°4 Is; K3; 276/81 
9 oo 21°! 74 378 14°7 4°5 17; K6; 276/84 
10 ©O 22°4 2 05°5 II‘l 6:8 13; K2; (121) 

II 00 23°3 73 20°7 14°0 10°2 18; K7; 276/85 
12 00 23°7 73 51°8 150 6:8 I9; K8; 276/86 
13 00 25°8 73 56°6 15°3 79 I10; Kg 

14 co 28°4 73 09°5 15°3 oy I12 

15 00 28°8 73 41°8 II'9 13°6 I11; Kro; (152) 
16 00 31°8 73 43°6 II'9 I1°3 114; K12; (176) 
17 co 31°8 74 09°6 — 13°6 I15; K13; 276/95 
18 00 33°! ay Ages — 26°1 116; K14 

19 00 33°9 74 27°7 _ 11°3 113 

20 00 34°2 73 02°1 14°8 5°7 113; K11; 276/94 
21 00 35°31 73 1671 13°5 5°7 117; K15; 275/2 
22 00 36°6 73 56°3 11°7 11°3 I20; K18; (220) 
23 00 36°7 73 184 15°4 3°4 118; K16; 276/98 
24 00 36°7 73 56°6 11°7 5°7 Iz1; K19; (222) 
2 00 37°! 73 54°7 11°8 g'I 123; K2o; (231) 
26 ©O 37°4 73 08-2 12°7 6°8 I19; K17; 275/3 
2 00 38°3 Pe me 11°9 13°6 I22; K21; 276/100 
28 00 39°5 73 08'5 13°6 4°5 124; 275/4 

29 00 39°9 73 59°! 12"! grI 125; K22; (242) 
30 00 42°3 74. 02°7 12'0 4°5 [26; K23; (256) 
31 00 43°2 73 15°8 13°9 3°4 127; 275/11 

32 ©0 43°4 69 28:3 —_ 3°4 la2 

33 00 44°0 73 21°6 14°3 3°4 128; 275/14 

34 00 44°0 74 O10 11°9 79 129; K24; (265) 
35 00 44°1 74 013 13°5 3°4 130; K25; 276/117 
36 00 44°6 74. 26°6 13°3 4°5 132; K27; 276/119 
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Taste III (continued) 
No. R.A. (i goo) Dee. (1900) Mpg d’ Reference 
h m ; 

ao 00 44°9 —74 04°0 12°2 4°5 131; K26 (269) 

38 ©O 45°0 70 261 16:2 gil a3 

39 00 45'5 73 541 11°8 17°0 133; 275/1 

40 00 47°! 73 29°! 52°93 12° 134; 276/133 

41 00 47°6 73 15°0 13°1 4°5 135; 275/27 

2 00 47°7 73 41°0 11°9 57 136; (290) 

43 00 483 12 225 II°g 10°2 137; K28; 275/30 
44 00 48°4 73 28°9 12°8 9'1 138; K29; 275/32 
45 00 49°! 72 43°2 12°7 11°3 139; K30; 276/139 
40 00 49°5 73 24°9 — 11°3 140; K31; 276/141 
47 00 49°6 73 53°3 12°7 g'1 141 

48 co 49°8 71 56°7 11°9 79 142; 275/42 

49 00 50°4 72 42°4 11°7 5°7 143; K32; (299) 
50 00 51°3 72 44°9 12°5 4°5 144; K33; (306) 
51 ©O 51°5 2 39°2 11°8 13°6 145 

52 00 52°0 74 O1°7 13°7 3°4 146; (156) 

53 00 §2°1 73 232 12°7 79 147; K34; 275/51 
34 00 52°8 73 O10 74 20°4 148; K35; (330) 
55 00 53°4 74 26°7 14°9 79 1114; 275/53 

56 00 54°! 72 49°0 II'5 4°5 149; 275/55 

57 00 54°4 73 57°72 — 23°8 152; K38; 276/165 
58 00 54°6 74 §2°2 14°] 57 151; K37; 276/164 
59 00 54°6 75 O10 12°0 17°0 I50; K36; (339) 
60 00 55°7 2 43°5 aaa 28°3 153; K39; (346) 
61 co 5674 72 §2°6 I2°1 9°! 154; K4o; (1611) 
62 co 56°6 72 54°9 I2°1 10°2 I55; K41; (1612) 
63 00 57:2 2 543 12°3 4°5 156; K42; 275/63 
64 ©0 57°7 73 52°7 13°5 4°5 157; K43; 276/175 
65 c0 57°38 73 17°0 13°4 g'I 158; 276/180 

66 00 58°6 73 05°7 I1'9 5°77 I60; 275/65 

67 00 58:9 72 08-9 12'0 17°0 162; K46 (361) 
68 00 59°0 74 27°0 11°9 17°0 159; K44; 276/184 
69 00 59°7 74 15°9 14°9 3°4 161; K45; 275/67 
70 00 59°9 72 48-6 12°3 34 163; K47; 275/69 
71 OI oor! 72 36°7 -- 79°3 164; K48; (371) 
72 o1 006 73 29s 11°8 79 I65; K49 

73 OI OF! 7O 52°9 rs 4°5 Ib2 

74 OIL O12 72. £39 12° 5°7 166; Ks50; 275/72 
75 OI oO1'9 12 37 — 17°0 167; Ks1; (395) 
76 OI 02°! 72 3eS 12°5 6:8 168; K52; (1624) 
77 OL. 63°13 73 49°5 13'0 ae | aa 
78 Of 03°7 72 13-4 13°0 rr-3 173; 276/214 

79 Ol 03°7 72 483 12°0 5°7 170; K54; 275/75 k 
So OI 04°3 73, 18:0 13°9 g'l 172; 276/212 

81 OI 04°6 72 38:8 14°1 CaF I71; K55; 276/211 
82 OI 04°7 72 17°8 11'O 12°5 177; K60; (411) 
83 O1 O5'1 72 52:8 12'0 113 176; Ks59; (416) 
54 OI 053 72. 3972 — 11°3 178; 275/79 

85 Ol 05°3 73 24°7 102 22°7 175; K58; (419) 
86 OI 05°3 73 4771 14°5 5°7 174; K57; 276/218 
87 o1 0673 72 174 12-7 4°5 179; K62; (422) 
88 o1 07°8 73 19°0 Pe7 3°4 I80; K63; 276/225 
59 Or 09°5 972 16°7 12°7 orl 185; K64; (1660) 
go or 09'S ae Sr2 12°5 10°2 183; (1655) 
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TABLe III (continued) 


No. R.A. (1900) Dec. (1900) Mpg a’ Reference 
h m - , 

g! o1 09°8 —73 37°7 14°8 5°7 182; 276/230 

92 o1 098 73 583 13°97 4°5 181; (1662) 

93 OI I0°0 73 59°5 —_— 5°7 1128; 276/228 

94 Or 112 73 481 _ 22°7 184; K65; (456) 

95 o1 11°8 71 §52°0 14°3 2°3 II30 

96 oI II°9 72 04°3 10°5 12°5 189; K69; (458) 

97 OI 12°0 72 49°8 —_ Q'I 186; K66; (460) 

98 OI 12°7 73 08°7 14°9 4°5 188; K68; 276/235 

99 O1 13°! 73 51°0 oe 90°7 187; K67; (465) 

100 OI 15°5 72 31°2 14°5 4°5 Igo 

101 OI 21°9 73 40°5 — 68:0 Ig1 

102 OI 22°1 71 42°3 14°3 4°5 II31 

103 Ol 22°7 73 45°7 — 34°0 Ig2 

104 O1 23°2 73 53°9 on 56°7 la4 

105 O1 27°3 74 04°3 —_ 34°0 193 (602) 

106 or 28:8 76 34°7 16°4 a3 Tas 

107 o1 28-9 73 55°6 oe 45°3 194 

108 OI 29°5 72 28:0 14°8 5°77 Ib3 

109 Or 31°5 74 40°4 I5‘I 3°4 la6 

110 Ol 32°2 7% 229 14°3 113 Ibs 

111 OI 33°5 76 04°4 13°4 20°4 Ia7 (643) 

112 OI 34°5 75 585 16-2 3°4 Ia8 

113 OI 47°9 Tm %35°3 II'9 22°7 Ib6 

114 or 48-9 74 50°9 12°1 34 lag 

115 Or 55°5 74 4271 — 3°4 Ib7 (796) 

116 OI 55°5 78 08:6 _— a3 lato 


Apparent photographic magnitudes.—Magnitudes were estimated on ADH 
plates (scale 68”/mm), one of 60™ exposure and three each of 30™, 10™ and 5™ 
exposures, and on four MF plates (scale 167"/mm) of 60™ exposure. Comparison 
stars were those of Sequence I. For the brighter stars c-h the magnitudes of 
Cox and Hallam (7) have been used; for stars k and fainter, King’s magnitudes 
(8). No high degree of accuracy is claimed for the magnitudes and the magnitudes 
of the sequence stars were taken to the nearest tenth. The ADH plates were 
first measured by making direct eye estimates. The measurements were repeated 
using a low-power eyepiece. This latter method was used for the MF plates. 
The mean differences between the MF results and others are given in Table IV. 


TABLE IV 


Mean differences between apparent magnitudes of clusters estimated on MF plates and other 
sources (Am= MEF magnitudes minus magnitudes from sources in column 1 of table) 


Source Am Remarks 
ADH 60™ and 30™ +2°9 Direct eye estimates: 73 clusters. 
ADH 60" and 30™ +1°2 Eye estimates with eyepiece: 45 clusters. 
ADH 10™ +0°5 Eye estimates with eyepiece: 48 clusters. 
ADH 5™ +02 Eye estimates with eyepiece. 
Shapley-Wilson +0°6 H.C. 275: 39 clusters: see below, 
Kron -0o'2 19 clusters; see below, 
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The quality of the 60™ plate was not quite as good as that of the 30™ plates 
and its limiting magnitude would be fainter by perhaps o™-5 or less. On all the 
long-exposure plates star y of Sequence I (m= 17-3) is about 1™-5—2"™ above the 
plate limit so that the latter is around m,,=19°0. ‘There is a slight systematic 
difference, less than o™-1, between the 60™ and 30™ cluster magnitudes, the 
former being the brighter, but all have been combined in ‘Table IV. 

‘The difficulties of estimating magnitudes of clusters by comparison with 
stars are accentuated with the pin-point stellar images of the ADH Baker—Schmidt. 
Looking directly at a plate, many of the clusters appeared quite stellar and it 
was rather surprising that, at least for these types, the differences from the MF 
plates were so large. ‘The differences might be largely subjective. However, 
Dr A. R. Hogg, during a brief visit to Armagh, kindly estimated the magnitudes 
of twenty clusters. ‘The mean difference between his estimates and the writer’s 
was less than o™-2, Hogg’s magnitudes being the brighter. In making direct 
eye estimates probably too much emphasis is placed on the size of a cluster. 
An eyepiece shows more of the structure; with it there is a tendency to 
concentrate on the nucleus, and attention is paid to density rather than size. 
By reducing the exposure time many clusters become more stellar-like. ‘This 
is not so for globular or circular objects; they remain disks with a different 
texture of the stars. It is much easier on MF plates to compare clusters with 
stars due to the chromatic aberration of the refractor which gives stellar images 
like the ADH out of focus. 

There are eighty-four objects in ‘Table II] common to the two papers of 
Shapley and Miss Wilson (ibid) and thirty-nine of these have magnitudes 
determined by these authors and the writer. ‘(he Harvard magnitudes are 
given as absolute magnitudes and allowance has been made for the different 
modulus used. ‘They are based on comparison with stars ‘‘in the numerous 
standard sequences’”’ set up in the Small Cloud. From the situation of the 
clusters with respect to Séquence I a guess was made as to which were based 
on this sequence. Due to the revision of the Harvard magnitudes it is likely 
that an average correction of —o™-1 should be added to the Harvard clusters 
here compared thus reducing the mean difference in the absolute magnitudes 
from +0°6 to +0°5. 

Considering the changes in appearance of clusters with exposure time and 
telescope it is not surprising to find differences between Armagh and Harvard. 
The latter were obtained from estimates on 8-inch Bache doublet plates (scale 
179"/mm) and 13-inch Cooke plates (scale 600”/mm). 

Magnitudes of twenty clusters in the Small Cloud have been published by 
Kron (9). He classified ten of these as globular and ten as non-globular. ‘The 
magnitudes of the twenty clusters were estimated on ADH plate 640. A sequence 
of four globular clusters established photoelectrically by Gascoigne and Kron 
(10) was used to find the magnitudes of the remaining six of the ten. It will be 
seen from Table IV that the mean difference between the writer’s MF magnitudes 
using the stars of Sequence I as standard and Kron’s ADH magnitudes using 
a sequence of globular clusters is —o-2. ‘That is, Kron’s absolute magnitudes 
are on the average fainter by this than the writer’s, whose values are roughly 
half-way between Kron’s and those of Shapley and Miss Wilson. In comparing 
with Kron, No. 54 (NGC 330) was omitted. The difference in this case 
amounted to +3™-1 but it was considered by the writer (and by Shapley and 
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Miss Wilson) to be too bright and scattered to measure. It was found extremely 
difficult to compare the globular clusters with stars, not only on account of the 
spread of the clusters but the different texture of their images. It was therefore 
surprising that the mean difference from Kron was the same for globular and 
non-globular objects. It is presumed that Kron compared the open clusters 
with stars of Sequence I and not with his globular cluster sequence. With one 
exception they lie beyond the range of the globular cluster sequence. Presumably, 
also, he used the magnitude system of Gascoigne and Kron (11) for Sequence I. 
If so, that would account for approximately o™-1 in the differences between his 
open clusters magnitudes and the present ones. Both Hogg and the writer 
estimated globular cluster magnitudes on a 30™ exposure ADH plate using 
Kron’s globular cluster sequence. ‘The same order of agreement was found 
as in Table IV. 

Considering the uncertainty in estimating integrated magnitudes of clusters, 
and that Kron’s magnitudes were determined from a single ADH plate compared 
to four MF plates of the writer, the agreement between the two is remarkably 
good. ‘This may be purely fortuitous. As will be seen from Table IV, different 
instruments and exposure times yield markedly different results and it is not 
possible to say what is the critical exposure time that gives correct magnitudes. 
Indeed what holds for one type of cluster may not hold for another. The 
behaviour of acluster image with exposure time is very different from that of astellar 
image. For example, on a 30™ exposure ADH plate, it would appear that cluster 
No. 85 (NGC 419) is very much brighter than the nearby star b of Sequence I 
(m=7:0); ona 5™ exposure ADH plate it is estimated to be over a magnitude 
fainter. As a matter of interest the behaviour of cluster No. 72 was followed on 
plates of different exposure until it was too faint to register. The results of the 
magnitude estimates are given in Table V. 


TaBLeE V 


Estimates of my, of Cluster No. 72 on different exposures 


Plate Exp. Mpg 
ADH 60™ 8-0 
ADH 30 7°9 
ADH ro™m 8-4 
ADH Cg g°0 
ADH* 208 12°0 

'‘ ADH** 5§ 130 
ers 60™ 11°8 


* On this plate the cluster was reduced to three images whose magnitudes measured 
on the Iris photometer were 13°08, 13°14 and 13°32 giving a total magnitude equivalent 
to m=12°0. 

** The cluster was at the limiting magnitude of the plate. So was star g of Sequence I 
(m= 13°0). 

*** Kron’s value is mpg = 12°0. 


Diameters. For the diameters of the clusters there is on the average, and 
individually with one or two exceptions, agreement with Kron. It is uncertain, 
however, how near to standard diameters they are; if a standard diameter is 
defined by the relative surface brightness they must certainly, at least for the 
globular clusters, be lower limits, For a given telescope the diameter depends 
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strongly on the exposure time. The behaviour of NGC 104 (47 Tucanae) was 
followed on ADH plates of exposures from 60™ to 58. ‘The measured diameters 
ranged from 60’ to 205” which, with a modulus of 14:17 (12), correspond to 118 
and 7 parsecs. On the 5% exposure plate the image of NGC 104 is similar to 
cluster images at long exposure in the Small Cloud. Due to the absence of 
background stars it is easy with this exposure to recognize some of the faint stars 
of the cluster outside the core and thus to classify it as globular. If this same 
image were, however, on a long exposure plate necessary to show many of the 
Cloud clusters the dense background of stars would mask the non-core stars. 
Only the core would be measured. From appearance it would be classified as 
a ‘‘cluster, probably globular ’’. 

On the 5% exposure plate the core was measured as 68”. ‘The exposure was 
too short to show any Cloud clusters. On a 60% exposure the globular cluster 
NGC 121 (No. 10) showed a measurable disk. Should this be similar to 
NGC 104, being seven times more distant, since the apparent diameters are 
inversely proportional to the distances, its core should be one-seventh. ‘The 
measured cores were 270” and 35”, or NGC 121 was one-eighth of NGC 104. The 
crude method indicates little difference between them. On a 30™ exposure the 
core of NGC 104 was measured as 680”. The mean of Nos. 8 and 15 was 163” 
on the same plate ; if we include No. 1, which is circular or globular, 177”. In this 
case the cores are only one-fourth that of NGC 104. Nos. 1 and 8 are the largest 
circular or globular clusters measured. ‘They are also away from the main body 
of the Cloud so that possibly more of their cores are recognizable. 

On a 30™ exposure plate the core of NGC 362 was measured as 340”. The 
mean of Nos. 27, 34, 43, 59, 67, 85 and 83 measured on the same plate was 88". 
Taking the modulus of NGC 362 to be 15-55 (11), if its dimensions are equal 
to this group its core should be between three and four times as large, which is 
what is found. The two globular clusters NGC 104 and 362 should be subject 
to the same amount of space absorption as the clusters of the Small Cloud. 

With a realization of the difficulties encountered in comparing galactic with 
Small Cloud clusters, Kron (9) has suggested that the latter may be smaller and 
less luminous than the former. His alternative is an increase in the modulus of 
the Cloud. There are so many uncertainties involved, including the estimation 
of magnitudes, the different methods of arriving at these for galactic and Cloud 
clusters, the different instruments involved, the difficulties of measuring diameters, 
and the effects of absorption in our galaxy, that it would seem at present 
impossible to come to any conclusion. From our present knowledge of the 
origin of clusters it would also not be possible to say whether one could expect 
a smaller galaxy would produce larger or smaller clusters than a larger galaxy. 

Shapley (13) considers that only one group of objects in the Small Cloud 
can be called a constellation in the sense as used in the Large Cloud. ‘This is 
the group NGC 456, 460 and 465 (Nos. 94, 97 and 99). ‘The group of Nos. tor, 
103 and 104 containing bright stars and nebulosity and with overall dimensions 
230 x go parsecs might be in the constellation class. Similar in appearance but 
less pronounced, also with bright stars and nebulosity, are Nos. tos and 107 
with dimensions 200 x go parsecs. On a 30™ exposure the nebulosity is well 
defined, measuring 45 x 30 parsecs. ‘These three groups form the outstanding 
connecting links with the wing of the Cloud. In Paper I, NGC 371 (No. 71), 
NGC 395 (No. 75) and NGC 346 (No. 60) were added to these. The name 
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stellar associations would probably be more appropriate, at least to these smaller 
groups, than constellations. Kron (iid) lists Nos. 18, 46, 60, 68, 71, 75, 94, 97 
and gg in this category. 

In Kron’s first paper (5) on the Small Cloud clusters he listed twelve clusters 
as globular, four questioned as globular and two as globular in appearance, 
In his second paper (9), applying a more rigorous test, namely that the brightest 
stars had to be of photographic magnitude of about 17, eight of the twelve and 
two of the four were retained as globular. ‘Taking into account Kron’s first list 
and the appearance of clusters as seen by Haffner and the writer, fourteen clusters 
were listed in Paper I as possibly globular. On the basis of integrated spectra, 
colours, magnitudes, apparent smooth structure and dimensions, Shapley (13) 
considers that only NGC 121, 416 and 419 can at present be listed as genuine 
globular clusters. NGC 339, 458 and H.C. 275/3 are classified by him as circular; 
he points out that M11 and NGC 2477, non-globular clusters in our own system, 
would appear circular to a remote external observer. 

Red magnitudes.—An attempt was made to get some idea of the redness of 
some of the clusters by comparison with unpublished red magnitudes of King’s 
of some stars of Sequence I. The plates were 103a.E with Wratten 25A filter; 
two of go™, one of 60™ and two of 30™ exposure. ‘The same effect was found 
as for the blue plates, namely, decreasing magnitudes with lower limiting magnitude 
of the plate. The differences between the go™ and 60™ plates were small and also 
between one of the 30™ plates which was of exceptionally good quality with a 
limiting magnitude about the same as the 60™. The clusters on the other 30™ 
plate averaged o™-6 fainter than on the first and this plate represents a fairly 
normal 30™ exposure red plate. The red magnitudes are not given in Table I] 
because it is not known what is a “‘standard’’ exposure time nor what are 
corresponding blue and red exposure times for comparison. Some qualitative 
ideas may, however, be obtained from the red plates. Colour indices were found 
by comparing the ‘‘normal’’ 30™ exposure red plate with a 5™ exposure blue 
plate. Due to the plates not being on the same centre and to difficulties of 
estimating magnitudes on one or the other or both, only twenty-seven clusters 
could be compared. All gave negative (blue minus red photographic) red 
indices. ‘They included the ten brightest open clusters listed by Kron but only 
one, No. 29, of his globular clusters. 

Kron’s globular cluster sequence consisted of the following objects; his P 
and V magnitudes are in brackets: NGC 419 (10°6; 10-1), 121 (1174; 10°7), 
416 (11-9; 11°4) and 339 (12:3; 11°8). It is clear from the red plates that the 
sequence is also in order of fainter red magnitudes. His other globular clusters 
were graded on the 103a.E plates with red filter. The results are given in 
Table VI which also includes Nos. 54, 82 and 96. 


TasBLe VI 


Estimates of red magnitudes (‘‘m,’’) according to Pickering’s method by comparison with Kron’s 
d r ° 
globular cluster sequence 


No. Miyy a No. Mpg ee No. Myo oe 

8 11°5§ 4167-339 34 12°0 ~ 339 67 12°0 ~ 339 
15 12'0 <50 43 12'0 <50 82 11-8 = 416°4°339 
31 12°3 <50 54 74 419°4°121 96 10°54 16°4°339 
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It will be seen from the Table that Nos. 54 and 96 are blue objects, the former 
being very blue. Nos. 8, 34, 67 and 82 are all about equal to NGC 339 in the 
red, Nos. 8 and 82 being somewhat brighter and Nos. 34 and 67 probably somewhat 
fainter. Since the first two are slightly brighter than NGC 339 in the blue and 
the last two about equal to it, it is likely that these four are red objects. The 
remaining clusters, Nos. 15, 31 and 43 are all about equal to NGC 339 in the blue 
but fainter in the red so it is not possible to say if there is or is not a few tenths 
of a magnitude of reddening. Little can be added to what has already been 
written. ‘There are other objects in the Small Cloud which could be considered 
as possibly globular; e.g. No. 58, a faint object in red, the texture of which 
resembles that of a globular cluster. For the present, Kron’s ten clusters together 
with Nos. 1, 82 and 113 are considered to be globular or circular. 

Distribution.—The clusters in the main part of the Cloud have all been plotted 
on Fig. 1. As measured from the cluster system the longest diagonal of the Small 
Cloud is 8040 parsecs; the shortest diagonal 5385 parsecs. The E—W extent is 
6445 parsecs; the N-S extent 7130 parsecs. These dimensions are larger 
than any so far published and probably represent the extreme limits except that 
star-counts would show a further slight extension in the direction of the wing. 
Shapley and Virginia Nail (14) have expressed some doubt that the wing may 
not be a part of the Cloud but a structure distinct from it. ‘The clusters extend, 
however, from the main body of the Cloud to the wing and end in it, leaving 
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Fic. 1.—Clusters in the Small Magellanic Cloud. 
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a gap between the two Clouds. A preliminary survey of ADH prism plates 
shows a bridge of emission objects, which follows closely the cluster bridge, 
from the main body of the Cloud to the wing and in it. There seems little doubt 
but that the wing is part of the structure of the Small Cloud. 

Little can be said about the distribution of the globular or circular clusters; 
their numbers are too small and their classification too uncertain. While it is 
easy in the majority of cases to recognize a cluster as open, there are a number 
of objects on the western side of the Cloud which are questionable. Neglecting 
the clusters towards the wing and the obviously outlying ones, an ellipse has 
been drawn to contain the main body of clusters. Within this are 63 per cent 
of all clusters; not counting twenty-two outlying clusters it contains 78 per cent 
of the remainder. The number of clusters in the four quadrants of the ellipse 
are 24, 17, 16 and 16 (see Fig. 1). ‘The excess of clusters on the wing side would, 
if distributed over the three sparse quadrants, give just an even number in all 
quadrants. 

Kron (15) considered that the non-globular clusters lie in a roughly spiral 
distribution and contained within an elliptic area. Within this are 57 per cent 
of his sixty-nine clusters. His field did not contain any of the outlying clusters. 
The ratio of his minor to major axis was 1 : 1-4 which, if interpreted as due to the 
tilt of a flat circular object, gave him a tilt of the Cloud of 45° and position angle 
of the major axis of the ellipse 60°. The position angle here of the major axis 
is also 60° but the ellipse drawn is much more elongated with ratio of minor to 
major axis 1:2°4, almost twice Kron’s. On the assumption of a flat circular 
object this would give a tilt of 65°. There is no sign of any spiral distribution 
of the non-globular clusters. From the distribution of clusters the Small Cloud 
could well be an irregular galaxy distorted by the tidal action of the Large Cloud. 
The tidal arm, culminating in the wing of the Cloud, while in the general direction 
of the Large Cloud, is not directed towards it but drops southward. This trailing 
arm might be a result of a past near approach of the two Clouds and the rotation 
of the Small Cloud. 


Armagh Observatory, 
Armagh, N. Ireland: 
1958 January 1. 
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Summary 


Some results from a survey of radio sources at 159 Mc/s are available for 
limited regions of the sky. These are discussed in relation to the spatial 
distribution of sources. 

A comparison of the individual sources shows that about one quarter of 
those in the previous survey at 81-5 Mc/s are confirmed, but that for the 
weaker ones the reliability is worse than was indicated in the list. A com- 
parison between the new positions and those found by Mills and Slee at 
85 Mc/s indicates fair agreement for sources with flux densities more than 
20 X 10°86 w.m.~* (c/s)~! at 85 Mc/s but poor agreement among the weaker 
sources. 

A simple plot of log N against log S, where N is the number of sources 
with flux density greater than S, gives a slope of —2-2 for one region and 
—2-7 for another, but these slopes may be affected to an unknown extent by 
confusion between nearby sources and by sources having finite angular 
diameters. 

A method of analysis depending on the statistical distribution of ampli- 
tudes on the records and independent of confusion errors shows that for the 
fainter sources the observations are consistent with a uniform distribution of 
sources. Between these and the intense sources a discrepancy in numbers 
remains. An argument may be devised to show that the 159 Mc/s obser- 
vations, when taken together with previous results at 81:5 Mc/s, are incon- 
sistent with a uniform distribution. This result is, however, based on data 
of which the statistical significance is not sufficiently great to be conclusive 
and further observations are needed. 





1. [ntroduction.—The most detailed surveys of radio sources which have 
been published are those of Shakeshaft et al. (1955) at 81-5 Mc/s (which will be 
referred to as the 2C survey), and of Mills and Slee (1957) at 85 Mc/s (referred 
to as the MS survey). More recently the aerial system used for the 2C survey 
has been adapted to a frequency of 159 Mc/s and a survev of radio sources (the 
3C survey) is approaching completion at this frequency. Some preliminary 
results from this survey, in particular those relating to the spatial distribution 
of radio sources, will be presented in this paper. 

The 159 Mc/s observations at Cambridge have been carried out, like those 
at 81-5 Mc/s, with the double interferometer described by Ryle and Hewish 
(1955). This has four aerials situated at the corners of a rectangle, 308 wave- 
lengths East-West by 27 wave-lengths North-South (at 159 Mc/s), so that the 
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overall reception pattern is that due to a single aerial, i.e. 8° by 1° between 
half-power points, but containing interference fringes in two planes at right 
angles. A phase-switching receiver (Ryle, 1952) is used so that any sources 
with angular extent comparable with the fringe separation are observed at reduced 
intensity. By phase-switching between the East pair and the West pair of aerials, 
sources with angular diameters greater than about 7’ are eliminated. A survey 
using such a technique will be called a ‘‘small diameter’’ survey. By phase- 
switching between the North pair and the South pair, it is possible to record 
sources with diameters up to about 1° and measure any angular diameters between 
the limits 2’ to 10’. This procedure gives a ‘‘ large-diameter ’”’ survey. 

The 85 Mc/s observations at Sydney were made with the Mills Cross aerial 
(Mills and Little, 1953), an instrument having a pencil beam reception pattern 
50’ by 50’ between half-power points when directed at the zenith. Sources of 
all angular diameters are detectable if their surface brightness is sufficiently 
high. Angular diameters can be measured by the increased width of the profile 
on the record, but the broadening is difficult to detect for sources much smaller 
than 50’. ‘Thus it is likely that some of the MS ‘“‘ point’’ sources will appear 
as extended sources in the 2C and 3C surveys. 

The principal purpose of the present paper is to set out the further information 
on the distribution of radio sources revealed by the 3C survey. Such information 
may be obtained by investigating the numbers of sources in different intensity 
ranges. If radio sources, for example, are distributed uniformly in space then 
N(S), the number with flux densities greater than S, is proportional to S~%?, 
so that a graph of log N against log S has a slope of — 1-5. Counts of 2C sources 
indicated a slope of nearly —3, implying a marked excess of faint sources, or a 
deficit of strong ones (Ryle and Scheuer, 1955). Mills and Slee (1957) on the 
other hand obtained a slope of — 1-8 and suggested that the discrepancy between 
this and —1-5 of a uniform distribution could be explained by instrumental 
effects. ‘hey also found poor agreement between the individual sources of the 
two surveys and suggested that both this and the steep slope of — 3 might be due 
to the effects of confusion between adjacent sources in the 2C survey. 

The various factors which may affect the counts of radio sources will first be 
summarized, and a brief comparison between the different surveys follows. 
This will indicate the extent to which agreement may be reached at present 
regarding the existence of individual sources. Finally the law N =constant x S~?? 
will be tested both by counts of individual sources and by a method (Scheuer, 
1957) which avoids the errors which may arise from confusion. 

2. Instrumental limitations and selection effects.—It is known that there are 
sources situated inside our Galaxy with a disk-like distribution near the galactic 
plane, among them being Cassiopeia A, the remnants of the supernovae of 
AD 1054, AD 1572 and AD 1604, and a number of extended sources of relatively 
low surface brightness (Cygnus X, Puppis A, IC 443, M 42, etc.). The 
inclusion of such sources would confuse the present issue since the log N — log S 
plot for a distribution of this type has a slope less steep than —1-5. The analysis 
is therefore restricted to regions away from the galactic equator. 

The factors affecting the detection of sources are discussed below. 

(i) Sensitivity.—The random noise on the records is not a serious limitation 
in either the 2C or 3C surveys; the noise level with a recording time constant 
of six seconds corresponds to a flux density of about 10-6 w.m.~? (c/s) 1 in each 
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en case*. Averaged over the time taken for a source to pass through the aerial 

ht beam, it becomes negligible. 

eS In the survey of Mills and Slee the r.m.s. noise, smoothed over the time of 

ed observation of a source (about three to five minutes) corresponds to a flux of 
Is, 3°5 x 10-*%6 w.m.~? (c/s) 1 and is the principal limitation of the survey. 
ey (ii) Confusion.—If the reception pattern of an aerial system is wide enough | 


to contain several sources of perceptible intensity, the records from the receiver 


rd no longer show the transits of a succession of individual sources but a ‘‘ confused ”’ 
- record is obtained on which several sources contribute to the deflection at any 
: instant. Only sources whose flux densities are some three or four times greater 
ial than that equivalent to the average level of confusion can be distinguished 
se unambiguously, and the measured flux densities of all sources are uncertain. 

of Confusion is the principal limitation of the 2C and 3C surveyst. The 
“j confusion level may conveniently be indicated by a flux density such that 
ile half the record shows deflections more, and half shows deflections less, 
er than the deflection corresponding to that flux density. In these terms the 
ar confusion level is about 15 x 10-76 w.m.~*(c/s)~? for the 2C survey, and about 

2x 10-*6w.m.~* (c/s)~4 for the 3C survey. 

- In the 2C catalogue the majority of sources having flux densities less than 
- 25 x 10-6 w.m.~*(c/s)~}, i.e. about one-third the total number, were bracketed 
Xl to indicate their unreliability. ‘The more recent analysis suggests that in fact 
cy the intensity at which individual sources can be reliably identified should be 

, increased to about 40 x 10°-*° w.m.~*(c/s)~!. This conclusion is confirmed by 
a a comparison of the individual sources in the 2C and 3C surveys. Of 143 
fe sources in the 2C catalogue between declinations 40° and 50°, 35 coincidences 


were found (allowing single lobe-shifts in right ascension and declination for 
the 2C positions) ; this result corresponds to a flux level of 40 x 10-** w.m.~*(c/s)~ 
a at 81-5 Mc/s. 


he Since the solid angle of the reception pattern of an aerial of given physical 
“ size is proportional to the square of the wave-length, nearly four times as many 
sources can be distinguished to the same level of confidence in the 3C survey 
¥ as in the 2C survey. The present analysis suggests that about 150 sources per 
~ steradian should be reliable in the 3C survey, a figure which is comparable with 
. the total number of sources in the MS survey, for which the confusion level 
(as defined above) and noise level are about the same. 
ws (ui) The effects of angular diameter.—'The equipment used for the MS survey 
™ detects all sources which produce a sufficient flux density within the reception 
rat pattern of 50’x 50’. Thus the limiting flux density is constant for sources 


of smaller than 50’ x 50’ and increases for sources with larger angular diameters. 
ly As explained above, the interferometer used in the 2C and 3C surveys 
“ normally observes sources of small angular diameter preferentially, i.e. in the 
Ss small diameter surveys sources with diameters greater than 6’ and 3’ respectively 
(in R.A.) are recorded with appreciably reduced intensity. ‘The corresponding 


1 ae ‘ . . ‘ 
limits for the large diameter surveys are 2° and 1° (in Dec.) respectively at the 
zenith and larger values at increased zenith angles. Intense sources can, of 
mn * Flux densities refer to the frequency of the relevant survey. 
nt + This limitation by confusion is not peculiar to interferometer systems, but is related in any 
h system to the solid angle over which sources may contribute to the output (cf. the blending of 


closely adjacent star images on an ordinary astronomical photograph). 
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course, be observed irrespective of their angular diameters by using one of the 
four arrays as a pencil beam instrument. 

Within any specified class of source, the nearest members will be the most 
intense and will at the same time have the largest angular diameters; in the 2C 
and 3C small diameter surveys there must therefore be to some extent a systematic 
bias against intense sources, and any apparent deficit of such sources must be 
corrected for this selection effect. 

There may be an added complication to statistical analysis, not confined to 
interferometric surveys, if radio sources occur predominantly in clusters. Such 
clusters would be resolved into individual sources when close by, but at great 
distances would appear as intense point sources. At intermediate positions the 
clusters would behave as sources of large angular diameter. An argument 
concerning the distribution which avoids these difficulties is presented in 
Section 4. 
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Fic. 1.—The ranges of sources visible on the 2C, 3C and MS surveys. The curves show the 
lower limits of the surveys. The approximate positions of several sources are shown on the diagram. 
[NGC 300 and NGC 6744 were observed by Mills during a special search (Mills, 1955).] 


(iv) Side-lobes of the reception pattern.—Since the most intense sources in 
the sky have flux densities over 1000 times greater than those of the faintest 
sources listed in the surveys under review, care must be taken to ensure that the 
passage of intense sources through subsidiary maxima of the reception pattern 
does not produce responses which are listed as true sources. ‘The danger of 
spurious responses is particularly high for aerials with ‘‘incompletely filled’’ 
apertures such as the Mills Cross, because, if the ratio of a subsidiary maximum 
to the principal maximum is R? (say) for a conventional aerial (or an interferometer 
using two similar arrays) it is only R for an “‘unfilled’’ aerial having a similar 
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resolving power. In order to reduce the spurious responses to an acceptable 
level it is therefore necessary to maintain a very much more accurate current 
distribution across the aperture in the latter case. 

In connection with the 2C survey a special investigation of the spurious 
responses was carried out and R? was measured accurately for subsidiary maxima 
down to one-tenthousandth of the main lobe. Mills and Slee were not able to 
make a similar check, and although they have taken considerable care to eliminate 
the side-lobe effects of intense sources it may be that some of them have escaped 
detection. Examples are given in Section 3 of two fairly intense MS sources 
which have not been found in the 3C survey: possible explanations of this are 
that the flux density of the sources falls exceptionally rapidly with frequency 
or that they are spurious responses. 

A further check on the reception of intense sources in subsidiary maxima of 
the 2C and 3C interferometer was made by comparing the apparent declination 
of the source as revealed by the periodicity of the interference fringes with that 
derived from the North-South reception pattern. 

The instrumental limits of the three surveys are illustrated in Fig. 1, in which 
the surface brightnesses of sources are plotted against their angular diameters. 
The regions of the plane accessible to the different instruments are marked. ‘The 
approximate characteristics of several sources are shown on the diagram. It 
will be seen that IC 443 is visible both to the 2C large diameter survey and the 
MS survey, but the Orion nebula (M 42), though observable by Mills and Slee, 
is outside the 2C limits. Mills and Slee also observe NGC 300 and NGC 6744 
which are undetectable on the 2C or 3C surveys. 

3. A comparison of the 3C and MS surveys.—The reduction of the 3C small 
diameter survey records has been completed for two regions: + 37° < Dec. < +52 , 
oo! <R.A.<24" and —10° < Dec. < +10°, 00" <R.A. <08", the latter region 
being common to the 2C, 3C and MS surveys. A comparison of the positions 
of the 2C and 3C sources shows little better agreement than that already found 
by Mills and Slee between their survey and the 2C survey*. A comparison of 
the 3C and MS surveys is more promising, however, and it may be noted that 
the limit set by confusion in the former is about the same as that set by noise 
in the latter (see Section 2 (ii) and (iii)). In cases where an MS source does 
not correspond to a source in the 3C list, the original 3C records have been 
consulted to check for possible errors. 

There appears to be a systematic difference between declinations quoted in 
the two surveys. If a correction, varying linearly between — 7’ at Dec. = + 10° 
and —19’ at Dec.= — 10°, is added to the 3C declinations, better agreement is 
found between the bright sources of the two surveys and in addition agreement 
to within the limits of error is reached between the 3C positions of Hydra A 
and Hercules A and those of the peculiar optical objects nearby. Because of 
the positional errors and the limited evidence this correction is not yet certain 
and it does not follow that all of the discrepancy lies in the 3C survey. 

A ‘‘coincidence’’ has been defined as agreement of the positions within the 
sum of the quoted errors in right ascension and within twice the sum of the 
quoted errors in declination. Owing to the possibility of picking the wrong 
interference lobe as the position of a 3C source, MS sources lying within similar 

* The reliability of all the surveys decreases towards the horizon, so that the agreement 


between the 2C and 3C surveys in the region —10°<Dec.<-+-10 is poorer than that for 
40° < Dec. < 50. 
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distances of the alternative positions on each side of the quoted one have been 
counted as coincidences. 
The results are as follows :— 


TABLE I (a) 


MS No. of No. of Mean expected 

Flux density MS observed no. of random 

(107-*86w.m.~*(c/s)~!) sources coincidences coincidences 
S'> 30 28 16 1°4 
30>S>20 49 13 3°6 
20>S> 15 59 4 4°3 
15>S 7 6 5°4 
Total — 207 39 14°7 


, 


When the ‘‘coincidence’’ range was extended to include sources within 
twice the sum of the right ascension errors and three times the declination errors, 
and a total of two lobe-shifts were allowed (but only one in declination) on each 
3C position, the increases in numbers of ‘‘coincidences’”’ were as follows: 


TABLE I (6) 


MS Observed Mean expected 
Flux density increase increase if random 
S> 30 4 I°4 
30>S>20 8 ee 
20>S>15 10 5°7 
I5>S 6 7°0 
Total 28 17°8 


In these tables the expected numbers have been calculated from the areas 
covered by the quoted errors for the sources. 

Eight of the MS sources in the highest intensity range are not accounted for 
in this comparison. One is almost certainly related to a 3C source of equivalent 
flux density about 1° away; another is near the southern limit of the 3C survey, 
where there are insufficient data to complete the analysis; and a third cannot 
be separated in the 3C survey from the side-lobes of Taurus A. Three are quoted 
in the MS list as being ‘‘ possibly extended ’’ and one of these has been found 
in the 3C large diameter survey. There is no 3C evidence for the remaining two. 

Assuming that the 67 observed ‘‘coincidences’’ out of 207 in Table I (a) 
and (b) are genuine, there are still 140 MS sources to be accounted for. In 
the case of 63 of these, the 3C evidence is inconclusive, but 77 sources definitely 
do not appear on the 3C small diameter survey records. This may mean that 
they have (i) angular diameters greater than 5’ or (ii) peculiar spectra, or (iii) that 
they are spurious side-lobe responses. In Section 4 evidence is presented to 
show that such a high proportion (77 out of 207) of extended sources is unlikely. 

A search has been conducted in the MS list around the positions of each 3C 
source to ascertain how many of the 3C sources can be explained in terms of 
MS sources. Of the 24 with flux densities >12 x 10-°* w.m.~?(c/s)~}, 16 can 
be directly related to MS “‘ point ’’ sources, 2 are presumably the sources referred 
to in the MS list as possible components of the ‘‘extended’’ source 00—o70 
(2C 50), 1 occurs close to the extended source o1 + 02, 3 are close to the source 
‘Taurus A and listed as unreliable, 1 is near 3 MS sources and is likely to be 
related to them, and 1 is not explicable in terms of MS sources. 
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Of the 138 fainter sources, 50 are confirmed directly by MS sources while 
a further 30 may originate in the same physical objects as nearby sources in the 
MS list but the evidence is not conclusive. The remaining 58 3C sources in 
the area have no corresponding sources in the MS list, but most of them have 
flux densities near the lower limit of the MS survey. 

Search has been made on the 3C large diameter records for MS sources 
listed as extended. Out of 18 sources, 3 were directly confirmed, in 7 cases the 
record was confused or there was a nearby source of about the right flux density 
but a small angular diameter, and in 8 cases no trace of the sources could be 
found. ‘Two notable sources absent on the 3C records are 06 +03 and 06+06. 
To investigate the discrepancy further, one of the four Cambridge aerials was 
set at declination +7° and a record made of the total power received. Relevant 
parts of the resulting trace are shown in Fig. 2. The size of the comparison 
sources reveals that the sensitivity should be quite adequate for the detection 
of 06 +03 and 06+06, so that their absence (if they are assumed to be genuine 
sources) can only be explained by very special peculiarities in their spectra. 











06+06 
he aan Hanne 
200,50) |(270, 
(120.31) +09°00' |+05°or’ 
+07°30' +077 14’ +0808 
2C12 |e ize! 
15-30 06-30 07:00 SID. TIME 


Fic. 2.—Parts of the record resulting from a total power run using one element of the Cambridge 
aerial as a pencil beam instrument, centred on declination -+-7° (as described in the text). On the 
right, the record in the vicinity of three MS sources is shown, together with their expected times of 
transit. On the left, for comparison is the record obtained from two 2C sources (2C1270 and 
2C1281). The source flux densities are quoted in brackets: for MS sources, tzo figures are quoted, 
representing the integrated (underlined) and peak flux densities: these flux densities are at a frequency 
of 85 Mc/s. The declinations of the sources are also shown. 


Whatever the final conclusions may be concerning the weaker sources, 
there is fair agreement on the whole about those with flux densities more than 
20 x 10-6 w.m.~? (c/s)~! at 85 Mc/s (10 x 107° w.m.~* (c/s)-! at 159 Mc/s). 

4. The relationship between the number of sources and the flux density.—The 
analysis of the 2C survey indicated that the numbers N(S) of radio sources 
above various levels S of flux density did not follow the law N(S)=constant x S~?2 
which would be expected if the sources were scattered uniformly through 
(Euclidean) space (Ryle and Scheuer, 1955). There appeared to be an excess 
of faint sources (or a deficit of intense ones). The data were of two kinds :— 


(i) Counts of sources—Owing to the effects of confusion, measured flux 
densities are subject to an uncertainty which becomes relatively greater for the 
faint sources (Bolton, 1956), and many of the weak ‘‘sources’’ recorded probably 
represent the joint effects of a number of sources with no dominant member, 
For these reasons the counts may be unreliable, 
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(ii) Counts of deflections.—The statistical properties of the deflections on - 
a chart which records the output of the receiving system can be predicted 
for any assumed number-flux density law (Scheuer, 1957). For the case 
N=constant x S~*? and a receiving system of the type used in the 2C and 
3C surveys the fraction of the record in which the-deflections have amplitudes 
between D and D+dD is P*(D)dD, given by 


P*(D) = x exp ( _ aw ?)\wDJ (wD) dw (1) 


where « is a scale factor determined by the constant in the number-flux density 
law and the reception pattern of the aerial, and J9(wD) is the Bessel function 
of zero order. 

Apart from the value of «, the relation (1) is independent of the form of the 
reception pattern; it is, however, affected if sources have finite angular diameters. 
By expanding the D axis and compressing the P*(D) axis by a factor 2a? in 
the graph of P(D), the scale factor « may be removed and the resulting graph 
is a unique one representing the function 


F(c) = , ot exp [—(4t)®2] Jq(ot) dt. (2) 


The parameter o= D/2a23 is a measure of the deflection as a multiple of the 
confusion level (which, when defined as in Section 2 (ii) by the median deflection, 
corresponds to o=0°83). 

For large deflections, each of which is generally dominated by the presence of 
a single intense source, the graph of P(D) against D is similar to a graph of 
—(dN/dS) against S, but for small deflections, where the overlapping of 
‘‘images’’ of sources on the record is important it approaches a Rayleigh 
distribution. 

By comparing the probability distribution (1) with the frequency histogram 
of observed deflections, the relationship NocS~*? may be tested in a manner 
independent of confusion and uncertainties in flux densities. 

Further data of both kinds are discussed below. 

(a) Counts of sources.—The 3C records for the region + 37° < Dec. < + 52° 
provide results free from interference and the spurious responses due to intense 
sources for an area of 0:25 steradians remote from the galactic equator. The 
graph of log N against log S for the 85 sources in this area has an average slope 
of — 2:2. Ne 

The second area for which the analysis is complete, that common to all the 
surveys, covers about 0-7 steradians, of which 0-5 steradians are further than 15° 10 
from the galactic equator. Fig. 3 shows graphs of log N against log S for this 
region (|b|>15°), both for the 2C and 3C sources. Each of the graphs is linear 
and the slopes are almost identical, viz. —2-5 for 2C and —2-7 for 3C. ‘The 
effects of confusion should be much larger for the 2C results, so that the agreement 
between the two slopes would seem to indicate that the systematic effects of 
confusion have not seriously affected them. There is the possibility, however, 
that the decreased confusion at 159 Mc/s may be compensated by an increased 
effect of the angular diameter of sources. 

(b) Counts of deflections.—Counts of deflections expressed in terms of o 
from both the 3C and 2C records are shown in Tables II and III. The areas 
of sky used do not contain any point within 15° of the galactic equator, and those 


2 mB OoOnpD et A 








XUM 





No. 2,1958 Evidence on the spatial distribution of radio sources 1g! 


parts of the records which show man-made interference, solar radiation, or the 
passage of very intense sources through side-lobes of the reception pattern* 
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Fic. 3.—Log N-log S plots for the two Cambridge surveys in the area common to the MS survey, 
excluding regions where |>| <15°. The solid angle involved is approximately 0-5 steradians, and 
N here refers to the actual number of sources observed. N= 10 corresponds to 1 source per 20 beam- 


zweidths. 
TABLE II (159 Mc/s) 
o © 0°869 2°197 4:355 6°505 8°65 10°64 12°67 14°61 > > 8-65 
to to to to to to to to to to 
0°869 2°197 4°355 6°505 8:65 10°64 12°67 14°61 20°45 20°45 8-65 20°45 
No. observed 1056 750 8153 29 8 2 I I ° 6 10 4 


No. computed 1056 749°5 137°8 29:2 12°3 Co 298 2:3 3°95 *4 21:2 15°85 


®=4°355 corresponds approximately to the deflection produced by a source of flux density 


10 10 *® w.m.~* (c/s)~! at the centre of the reception pattern. 


TABLE IIT (81-5 Mc/s) 


a oO 9°79 1°57 2°36 3°15 3°93 4°72 5°51 6:29 7°08 > > 39 

to to to to to to to to to to to 

O79 157 2°36 3°15 3°93 4°72 5°51 6:29 7-08 7°87 7°87 3°15 7°87 

No. observed 2443 2040 419 131 46 35 14 7; ¢ wg 176 119 
No. computed 2443 1816 497 175 85 48 30. 21 16 12 64 277 «213 


* The justification for the elimination of the side-lobe responses of the intense sources in 
Cassiopeia, Cygnus, Taurus and Virgo rests on the fact that their flux densities are far greater than 
those expected from any reasonable extrapolation of the log N-log S curve of other sources. In 
addition, only the source Virgo A lies more than 15° from the galactic plane. Whatever the de- 
clination setting of the aerial, most of the large deflections recorded are due to these four sources, 
and the statistics of large deflection amplitudes would be overwhelmed by their presence if all their 
side-lobe responses were included, 
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were omitted from the analysis. ‘he 159 Mc/s data comprise 2006 readings 
of deflection from an area of about 1-7 steradians in the region — 10 < Dec. < +56, 
08! <R.A.<16". The 81-5 Mc/s readings refer to a similar region but the 
amplitudes were noted at nearly ten points per beamwidth and the statistical 
errors are therefore about three times larger than the numbers would appear to 
indicate. 

The scale factor in (1) was adjusted by equalizing the observed and computed 
fractions of deflections in the lowest range of deflections in the table (about 
half the total number of deflections in each case)*. ‘Table II shows good 
agreement between the observations and the predictions for a uniform distribution 
up to o about 8-65 (the precision of the agreement is fortuitous). Counts of the 
small deflections at smaller intervals of D also show agreement with the predicted 
values. For the larger deflections there are fewer data but the observed numbers 
are consistently less than the predicted numbers except in the highest range 
(a > 20°45), which should perhaps be considered separately since the deflections 
are nearly all due to the intense source Virgo A (observed in the primary lobe of 
the reception pattern). It will be seen that there are only 10 deflection amplitudes 
with o> 8-65 instead of the predicted 21 and in the range o=8-65 to 20°45 
only 4 instead of 16. A similar effect can be seen in Table III with the 2C results 
but it sets in much sooner, at about o=3:15. ‘The following argument shows 
that a difference of this kind is to be expected if the discrepancy is due to a true 
deficit of intense sources. If the aerial, when operating at 81-5 Mc/s, has a 
reception pattern covering 1/Nth of the sky the confusion level (a = 0-83) corre- 
sponds approximately to the flux density of the Nth source in order of decreasing 
intensity. ‘The same aerial structure operating at twice this frequency has a 
reception pattern covering } Nth of the sky, so that o=0-83 corresponds to the 
flux density of the 4Nth source, and the Nth source corresponds to a higher 
value of «. Assuming for the sake of the argument that No S~*?, this value 
would be o=0°83 x 473=0°83 x 2°5. ‘hus the sources contributing most to 
the deflections in the neighbourhood of o = 3-2 at 81-5 Mc/s are the same sources 
which are important near o = 3:2 x 2°5 = 8-0 at 159 Mc/s. 

The method of analysis just described was developed especially to avoid the 
disturbing effects of confusion between adjacent sources. It does not allow 
for the effects of angular diameters. Ryle and Scheuer (1955) showed that the 
source counts on 81-5 Mc/s could not be explained by the effects of angular 
diameters alone, but since the counts were also very severely affected by 
confusion, the whole question must be re-examined. 

(c) The effects due to extended sources.—TVhe existence of extended sources 
at low galactic latitudes is well known but observations away from the plane 
have shown few extended sources and most of those investigated have had angular 
diameters less than the limits of measurement. ‘The published data are 
summarized below together with some new results provided by the 3C large 
diameter survey. Since the most intense sources in any given class will have 
the largest angular sizes these data may be used to derive the errors in the flux 
densities measured with an interferometer of given resolving power. 

(i) Morris, Palmer and Thompson (1957) have recently investigated the 
angular diameters of some sources in the catalogue of Hanbury Brown and 


*In an earlier paper by Ryle and Scheuer (1955) the curves were scaled by fitting at large 


values of D, 
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Ags Hazard (HBH, 1953). ‘The catalogue covers a region between declinations + 38° 

56, and +68°. ‘They find that ‘‘ with the exception of the Andromeda nebula all 

the the intense sources in the catalogue of Hanbury Brown and Hazard with galactic 

cal latitudes greater than 10° have angular diameters <5’’’. ‘They quote results 

to for five sources : | 
HBH 6 (NGC 1275) 2’4+0'5, 

ted HBH 10 <12", | 

a HBH 11 <12”, 

om HBH 12 4/20" + 2’, 

the HBH 18 <ev". 

ted The Andromeda nebula is the only normal galaxy of which the angular 

ers diameter is known. Baldwin (1954) finds it to be larger than the optical and 

nge equal to 2°-3 between half-power points. 

Ons (ii) ‘The angular diameters of other sources known not to belong to classes 

- of concentrated to the galactic plane have been measured: 

es Cygnus A: ‘Two disks each of about 40” diameter separated by go” (Jennison 

“45 and Das Gupta 1956). 

- Virgo A: Central concentration 5’ x 2’:5 (Mills, 1953) with 80 per cent 

mi of the flux surrounded by a disk about 50’ in diameter 

ne (Baldwin and Smith, 1956). 


Centaurus A: Central concentration 3’ in diameter (Mills, 1953), with extended 


“a source about 1°-5 in diameter. 
sa Hydra A: 1-6 (Carter, 1955). 
the (iii) Baldwin and Archer (private communication) have determined upper 
her limits to the angular diameters of the 27 brightest sources in the region 
lue —10° < Dec. < +47°, |b|>20°, and found that 25 have angular diameters less 
to than 5’, 21 less than 3’:5 and 10 less than 2:5. 
ces ‘The observations show immediately that unless faint sources differ radically 
from intense ones by having much larger angular diameters (the opposite being 
the expected) then the fraction of sources with angular diameters more than 5’ 
low is considerably less than one-third (see Section 3). ‘hey also suggest that for 
the the great majority of the relevant sources the reduction in apparent intensity 
lar due to finite angular diameters is small, even with a spacing of 308 wave-lengths. 
by For the 2C observations, in which the interferometer spacing was only 
_157 wave-lengths, the reduction should be smaller still. 
ces It is clear from an inspection of Fig. 1 that normal galaxies having a surface 
ane brightness similar to that of M 31 leave no trace on the 2C or 3C small diameter 
ilar survey records and therefore cannot affect the counts of individual sources. 
are A calculation shows that the effect of normal galaxies on the P*(D) histogram 
rge is also negligible. 
ave The results of Baldwin and Archer quoted in (iii) above refer to a repre- 
Jux sentative sample of sources and may therefore be used to estimate the effect 
of angular diameters on the data in Table I]. ‘The number of sources per 
the steradian is about the same as the expected number of independent deflections 
and with o > 8-65, per steradian, at 159 Mc/s, and it is reasonable to assume that 
-_ the sources which should be responsible for deflections with o > 8-65 have the 


same distribution of angular diameters as the sources observed by Baldwin and 


a3 











194 D. O. Edge, P. A. G. Scheuer and F. R. Shakeshaft Vol. 118 


Archer. With this assumption, it is found that up to 8 of the 21-2 expected 
deflections with o > 8-65 might be reduced below o=8-65 owing to the resolution 
of sources by the interferometer. The remaining discrepancy between 10 
observed and 13:2 expected large deflections in Table II would not be significant. 
The estimate just given is an upper limit to the effect of angular diameters, since 
the upper limits as quoted by Baldwin and Archer have been taken as actual 
angular diameters. On 81-5 Mc/s the interference lobes in the reception pattern 
are more widely spaced and a smaller reduction in the apparent intensity would 
occur, so it would be very difficult to explain the results in Table III as an effect 
of angular diameters. 

Quite apart from considerations of individual source diameters it is possible 
to use a general argument, based on the geometrical similarity between the 
reception patterns of the 2C and 3C surveys, to show that the data in Tables II 
and III taken as a whole cannot be explained on the hypothesis that the 
distribution of sources is uniform on a large scale, whatever the angular extent 
or clustering tendencies of these sources may be. For assume a uniform 
distribution and consider : 

(A) the sources in the spherical shell of thickness 5R at a distance R, 
observed at a frequency of 81-5 Mc/s; and 

(B) the sources in the spherical shell of the same thickness 5R at a distance 
(159/81-5)R, observed at 159 Mc/s. 

The two situations are geometrically similar, for the angular dimensions 
of the reception pattern and the interference lobes, of the sources and the clusters 
of sources (if any) are all reduced by the same factor (81-5/159) in case (B). 
The mean number of sources in the reception pattern, however, remains the 
same and the statistical properties of the receiver response are therefore the 
same in both cases. 

The identity holds if we consider the net effect of all the shells filling space 
in case (A) and the corresponding shells in case (B). These latter are spread 
throughout space but only fill a fraction (81-5/159) of it, since the thickness 
of each shell is the same as in case (A). We see therefore that the 159 Mc/s 
observations for a uniform distribution should be those which would have been 
obtained at 81-5 Mc/s from a distribution of sources with a space density 
increased by a factor (159/81°5). Now, if the density of sources is increased, 
the flux density which corresponds to any given o in Tables II and III is also 
increased, so that the discrepancy between the observed and computed numbers 
(if due to angular diameters or clustering) should become apparent at a lower 
value of o at 159 Mc/s than at 81-5 Mc/s; Tables II and III indicate that the 
reverse is true. 

This reasoning shows that the data in Tables II and III are incompatible 
with a distribution of sources in space which is uniform even if the sources have 
finite angular diameters or exist in clusters. ‘The data are compatible with a 
model in which there are fewer intense sources than expected on the basis of 
the number of faint ones. An estimate will now be made of the intensity level 
at which this deficit appears. In Table II the deviation occurs at about o = 8-65, 
corresponding to a flux density of 20x 10~%w.m.~*(c/s)-1. The sources 
responsible for these deflections must, in general, have somewhat greater fluxes, 
since the point at which a deflection is read on the record rarely coincides exactly 
with the transit of a source through the peak of the reception pattern, Thus the 
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d deficit of intense sources chiefly concerns those with flux densities greater than 
n about 30 x 10-*6 w.m.~?(c/s)—! at 159 Mc/s. 
6 If the deficit is isotropic there are in the whole sky some 100 sources brighter 
t. than 30 x 10-8 w.m.~* (c/s)! whereas, on the basis of a uniform distribution 
‘e and the observed numbers of fainter sources, about 200 would be expected. 
al 5. Summary.—An attempt has been made in this paper to present the 
n evidence concerning the distribution of radio sources in space which has been 
ld obtained from a new survey at a frequency of 159 Mc/s. When the reduction 
ct of the records has been completed, more flux densities will be available for 
plotting number counts, but it is doubtful whether this additional information 
le will alter the situation materially. 
1e The evidence may be summarized as follows :— 
i (i) Counts of sources in the small diameter survey give a log N-log S graph 
ne of slope —2-7 for a region of area 0-7 steradians, in good agreement with the 
nt slope of — 2-5 obtained from the 81-5 Mc/s observations of the same area. For 
m another region (0°25 steradians) the slope is —2:2. These counts are, however, 
subject to two important instrumental effects, namely confusion, and the 
R, resolution of sources with appreciable angular diameters, confusion being more 


severe on the lower frequency, and the effect of angular diameters on the higher. 


ce (ii) By a direct comparison of observed and calculated distributions of the 
deflection amplitudes on interferometer records, the errors due to confusion 
as can be avoided. A deficit of large deflections is found. 
ag (iii) The deficit of large deflections found on 159 Mc/s could possibly be 
). explained as the effect of angular diameters, but it would be very difficult to 
: account for the 81-5 Mc/s observations in this way. 
” (iv) A general argument shows that the 159 and 81-5 Mc/s data, taken 
- together, cannot be explained in terms of a uniform distribution in space, whatever 
- the angular diameters or clustering tendencies of the sources may be. 
pss (v) For the smaller deflections, the distribution of deflection amplitudes in 
c/s Table II agrees well with that expected from a uniform distribution of sources 
en in space. 
ity If complete reliance may be placed on the data, it follows that there are fewer 
ed, sources with flux densities > 30x 10-76 w.m.~*(c/s)“! at 159 Mc/s than one 
Iso would expect if the distribution of sources in space were uniform. This 
ers conclusion, however, depends on data of which the statistical significance is 
ver not very great, and which are subject to errors in reading and calibration. In 
the addition, doubts may be entertained about the procedure of excluding side-lobe 
responses of the intense sources in Cassiopeia, Cygnus, Taurus and Virgo, when 
ble counting deflections. ° 
ave It has been noted in Section 4, that the departure from a uniform distribution, 
ha if real and isotropic, shows up mainly among the brightest hundred or so sources 
: of in the sky. Only about half of these sources are situated in regions of the sky 
vel remote from the galactic equator and accessible to any one radio telescope. 
65, The use of elaborate statistical arguments in connection with a fairly small number 
ces such as fifty seems undesirable and should soon become unnecessary. ‘There 
xes, are now several very powerful radio telescopes in, or about to begin, operation 
ctly and it should therefore be possible to reach complete agreement concerning the 
the existence and flux densities of the brightest few hundred sources. 
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PESIOME IIOJJIMHHbIX DOKJIAIOB, B TMEPEBOQE HA PYCCKHH A3bIK 


MHTEPPEPOMETPHYECKHE HS3MEPEHHA JVIMH BOJIH. 
IV. TOUHOCTbh H3MEPEHHBIX JVIMH BOJIH 


M.T. Adam u C. Huxo.vc 


XapakTepucTHKa COJIHCUHBIX JIMHHK MOPIOUlCHHA HAMHOTO OTJIMYaeTCH OT COOTBCTCTBYIOILAX 
slaSopaTOpHbIix JIMHHHA HW TOUHOCTB ompetesICHHA JWIHHbI BOJIHbI NCPBbiX SHAYUHTCIIBHO HHPKC. 
IIpu 9TOM BO3SMO7KHO, UTO HCKOTOPAaH YaCTh KpaCHOTO CMCILIICHHA H H3MCHCHHE KpacHOoro CMCIIICHHA 
C MHTCHCHBHOCTbIO JIMHHH, BOSHHKaroT TIpH MHKPOMCTPH4CCKOM H3MeCpeHHH LIHpOKHX JIMHUL 
MOrsIOuUleHHA HW Y3KHAX JIMHHM IMHCCHH. HecnejjoBaHbl TOUHOCTA MH MOCIIC/OBaTCJIBHOCTL PeC3YyJIb- 
TaTOB TMIOJIYYUCHHbIX TIPH H3SMCpeHHH JVIHHbI JIYPOBbIX H COJIHCUHbIX BOJIH MCTOJIOM KOJIbIICBOH 
/[HarpaMMBIl. Ilo 9THM HCCHeOBAaHHAM OKa3bIBaeTCH, UYTO MOXKHO NOJYUYATh JOCTOBCpHbIc 
KpaCHbie CMeCLucHHH C TOUHOCTbIO nmpeBocxojAeH 0,001 A. 


MHTEPREPOMETPHUECKHE HS3MEPEHHA JUIMH BOJIH. 
V. MCTOJKOBAHHE KPACHBIX COJIHEUHbIX CMENIEHHH 
MPM MWOMOUM PAITMAJIBHBIX ILOTOKOB 


M. I. Adam 


PaHbule CUHTasIOCb, UTO TEOPHA BIIMAHHA OTHOCHTeIbHOLO 9¢pd*eKTa pajlHa/IbHbIX MOTOKOB Ha 
KpaCHbie COJIHCUHbIe CMCLICHHA IPHBOJAMT K TOTOKaM HCBO3MO7KHO Gobo CKOpocTH. OjtHako, 
OuaroqapxH HOBO MOjeNM aTMOCdepbI, GasHpoBaHHOii Ha rpaHyllax j{MaMeTpOM B 0,4 CeK H 
OoubUIOrO KOHTpacta APKOCTeii MO OTHOIWICHHIO K MHTeprpaHyJIOBbIM tmpoctpaHcTBaM, LIperep 
Halles, UTO CMCILICHHA MO?KHO OOLACHHTh MPH MOTOUHBIX CKOPOCTAX BCero B 2-3 KM,ceK. Ero 
feopeTHUeCKHe MpeCKa3zaHHA JLOCTaTOUHO COriacyloTCA C HaOJOAaeMbIMH H3MCHCHHAMH CMCICHHs 
OT 1eHTpa K MepHdepHu. OTy TeOpHIoO MO7KHO TaKxKe IIpOBePHTb MPH MOMOLIM OKCHopACKHX 
H3MepeHHit {VIMHbI BOJIH B WeHTpe COHeUHOTO JHCKa MW B BakyyMHOii lyre. TlocneqHue alot 
BO3MO7KHOCTh OMPeWCIHTh PaBHOEHCTBYIONLYIO MOTOUHOM CKOPOCTH B 3aBHCHMOCTH OT HHTeH- 
CHBHOCTH JIMHHM; HO NOJyYeHHbIe TaKHM OOpa30M CKOPOCTH He CoriacyloTcH C MpejckKasaHHAMii 
IIIperepa. GaTem paccMoTpeHa 3aBHCHMOCTb CMeLI[CHHIi OT JVIMH BOJIH, HOJIb3yACh H3MepeCHHAMH 
aJliereHCKoit O6cepBaTopHuH vu Gropo ctaHAapToB CIIIA; gtoT onbiT Tak7Ke He NOs WeprKMBacT 
TCOPHH pajMaJIbHbIX MOTOKOB. 


MNEPHODbI MW KPMBBIE CBETHMOCTH HECKOJIBKUX IIEPEMEHHbIX 3BE3JI 
B CKOINJIEHMW NGC 121 


A. J. Taxkepat 


Ilomb3yxicb 82 CHHMKaMH WiapoBoro CkomieHHun NGC 121 Bo BHemHel uactH Manoro Maren- 
JaHOBa OOaKa 3a JICBATHJICTHHA NepHon, XOTA GONbUIMHCTBO HX ObIIO CHATO B TeueHHe 3% 
NOCIE€LOBATeIbHbIX FOOB, BbIBECHbI MePpHOAbI HM KPHBbIec CBCTHMOCTH NATH MepeMeHHbIX 3Be3/1, 
HaXO{AIMXCH B pesesax yr B 1’ OT WeHTpa ckonIeHusA. M3 Hux TpH ABIIAIOTCA NepeMeCHHbIMH 
KOPOTKOrO NepHojla (HX MepHOsbI—OT 0,50 {10 0,64 jtHH), THNa a HIM / no bese, c OObIuHbIMH 
COOTHOWICHHAMH = M@KJIY MepHOAOM, aMMIMTyoi Hu dopmoii KpHBO cBeTHMocTH. I[lepnosbI 
OCT@JIBHBIX JIByX, IIPHHajWIe7KalllMX K CaMbIM KpaCHbIM HM CBETJIbIM 3BC34aM CKOMJICHHA, OKa3bl- 
BaloTCA 140,2 H 112,4 cyToK. M3mMeHeHuA CBCTHMOCTH XapaKTepH3yIOTCA CTYMCHUAaTbIMH OLCHKa- 
MH 110 OTHOUICHHiO K ( OCHOBHOMY )») TIOPAKY MOCJIeOBaTeJIBHOCTH. 

Ha cHUMKax OOHApyxKeHbI elle 4YeTLIPe MePeMeEHHbIX 3Be3{bI KOPOTKOTO MepHojla, pacnos1o- 
XKEHHbIC CIMUIKOM {1a/1CKO OT CKOMJICHHA IIA TOrO, YTOOLI MpHHasiexkaTb K HEMY. OTH 3Be3/IbI 
AOCTaTOUHO CaGbI, YUTOOLI MOXKHO OBUIO HX PaCCMATPHBaTL Kak NepBbIi MpHMep MepeMeHHbIX 
RR Lyrae B o6ujem none Manoro MaresinanoBa o6nakKa. 
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AHAJIN3 PAIMAJIBHbIX CKOPOCTEM ITAJIBHHX 3BE3], THIIA B 
M. B. ®ucm u A. A. T3xxepeii 


PajiMasIbHbie CKOpocTH 314 3Be3j, THNa B, paHHero noApasyCIeHHs, C BbIUMCJICHHbIMH pac- 
CTOAHHAMH (i110 (boOTOJIEKTPHYeCKHM BeTaM WM KJlaccudbuKaunMn Mopraxa-Kunana) mpeppiuaromn- 
MH I KIC, AHAIMSMPYIOTCA C IICJIbIO ONIpeyereHHA WUddepenuMasbHOrO rasaKTHUeCcKOrO Bpalllenns. 
Bkarovenve cKOpocteif 1O>KHbIX 3Be3jI, COrmacHO OMmpeyeseHHAM payikuIMdpcKoi oOcepBatopun, 
jlaeT ropa3yo Goulee OHOpoAHOe pacnpesenenMe 10 WONroTe, 4YeM MPH MpeAbIAyUMx aHassax. 
Ha paccToAHHAX MIpeBbIIMalOuIMX 2 KIC jlaxike 3aMeTHO He XBaTaeT 3Be3/, CeBepHOroO TosTyMapnA, 

TiulaTenbHO MpoaHasIM3MpOBaHbI 3Be3HbIe HM Me*KLY3BC3HbIe CKODOCTH py pa3lIMuHDix 
paccTosHnit (oTmes §, TaOsHMua ID). Orpyyarenbypii tepmun K (orgen 6) nmomyuaetca yA 
OOaCTH OT I {0 2 KIC, HO COrmacHo GOo.iee JjJasleKHM rpylilaM HeT yka3aHHA, UTO STO ABJIACTCA pesysE- 
vaTOM OOujero cKaTHA Tanaktuku. Bkpatue oOcyxkJ[aeTCH BO3MO>KHOCTb CBA3H MexKity Hads10- 
aCMbIM ABJICHHeEM H OTKJIOHeHHeM BepumMHbI. Tepmun K sp oOnactu oT I M0 2 Kic 3HauMTesIbHO 
VMeHbIUaeTCA KOrla IIpHHUMaecTcH B COOOparxkKeHHe CHcTemaTHUeCKOe HcMpaBJIeHHe, UTO BHMMO 
IIPHMeHHMO K cKOpocTsaM JInKKa ca6pix 3Be37, THNa B (oTyem 3). [log BamMaAHMeM 9TOrO UcMpaB- 
JICHUA TAKOKeE YBCTHUMBAIOTCA BeJINUHHbI, BbIBCJICHHbIe B Pe3yJIbTaTe IIPHMCHCHHA MOCTOAHHOI 
A Oopta. Tpynna 3Be3y yqasieHHbix Gosee 4eM Ha 3,0 KIC, 1aeT HCKJIOUNTeEJIbHO HH3K0e 3HAaUeHHC 
A, HO 9T0 OObACHACTCA, TIABHbIM OOPasoM, Kak pesyJIbTaT OTOOPAa, MIPHBOAUIero K HecopasmMepHO 
OoubIUIO!T MpONOpuHH B ITO pyle mpeyBesIM4eHHbIX paccToAHHi. 

PaccesHue CKopoctelt BO BHYTpeHHelt uactu TanakTHKH He ABJIAeTCH HECOpasMePpHbIM (OTe 7). 

IIpumensierca MpHOIM3HTeIbHbI MHelHEIM pacuér Kamna (oTy[es 8) It BbIBOJIa BeJIMUHHbI 
uoctosHHOH A Oopta. IIpuustra okoHUaTeIbHaA BesIMUMHa A~ 17,5--1,5 KM/CeK.KIIC. 

Jia Bprpoya BenmmunHb! Ro, paccronHund oO weHtTpa TasakTHKH, IPHMeHAIOTCA J{Ba YaCTHUHO 
HCSaBMCHMbix Meroua. Handoliee BepoxATHaH BesIMUMHa—-8,9 KIC (OTe 10, Tad. IX). Orno- 
IeHve Ocelt 9JUIMMCOUa CKOpoctelt JIA asIbHUXx 3Be3;, THNa B BHAHMO CJIeMyeT CUMTATh PaBHbIM 
0,5. 

ITposepvetcH KammMOpaitHwA urKalibr cBeTumMoctnh Moprana-Kuuana (ores 11) JIA 3Re3;, B 
30He oT Bs 0 Oa,s. 


Bprpoyutcs KpuBan WuddepenyHasibHoro BpaleHHA w(R) (oTAes 12). Merkjty 9THM penieH- 
MeCM HM peSVJIbTAaTaMH paj{MO-HaGINO/CHHM HET CYI[CCTRCHHOI! paSHHIbI H MmpeywouaraetcH, 4To 
w” NOUO*KNTeIbHa B OOACTH UccuIestyeMbIx R, 


SBESJIHbIE PPYIIbL LL. PPYIMLbI SBESIL BOJIBIIONW CKOPOCTH ¢ PEPKYJIECA, 
« HHTEEUA UH 61 JIEBETA 


Oauns A. Deeen 


I[pocmotpeH KaTasIor CeMHCOT 3Be3, OObLUOM CKOPOCTH, C I@JIbIO HaMTH 3Be3j[bl C ABHO 
ONpeesIeHHbIM JIBYOKEHHeEM B MpoctpaHcTRe, mosoOnEImM £ Tepkymeca (23 3Be3yb1), « Mnyecua 
(14 3Be3q) u 61 JIe6ens (15 sBe3n). JiBwKeHuA B MpocTpaHcTBe 9THX rpynn MoryT GpITb 0600- 
IN[€HbI CulexyIOUuAM OOpasom (rpynupbt Puag wu Cupuyca BKINOUeHbI JIA CpaBHeHHA): 


Ipynna A D Vt U V W 
, KM/CeK 
Cupuyca 204 44M — 42,7 18,4 f—14 fe) —12) 
Muay 64 18M t+ 795 44,0 | +40 —18 — 2 
¢ Her 54 48m — 16,5 7455 +54 —45 —26 
« Ind 64 44M + 6,7° 88,5 +79 —39 + 6 
61 Cyg 64 28M + 0,3° 108,0 +92 —§3 -— 6], 


rye 3Hauenna (A, D) nperctaBsa1oT coGoii TOUKM CXOAMMOCTH KaxKyIIMxcA ABMKeHHi, a V;— 
OOLLYIO CKOPOCTh B MpocTpaHcTBe KaxkKOi rpynbl no oTHOWeHHIO K Conny (ipHuem HalipaB- 
JICHMA BEKTOPHasIbHBIX CKOpocteii B mpocrpanicree: /=148°, b=o0°; 1=58°, b=o° ub=-+90 
cooTBetcTBeHHo yin U,V u W). Tlouckn B Tom xke KaTasiore 3Be3, MOTyII{Mx MpHHajVIenkaTh K 
CHHTeETHUCCKHM rpynam c TOUKOI CxOAHMOCTH B b= —90’, a V;=74,53; 88,5 M 108 KM/CeK, j[aloT 
ykKa3aHua Ha TO, UTO BCerO 2 JIM 3 ClIYUaMHBIX 3Be3/[bI MOryT ObITb BKJIOUCHBI B KarK{yI0 Ppylilly. 

Jiuarpamma WBeta-CBeTHMOCTH 3Be3, BCeX Tpex rpynn OonbUO CKOpocTH MoA0GHa jHarpamMe 
rasakTaueckoro cKommeHHa M 67. 
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QJLIMNCOMT, CKOPOCTEH JUIA 3BE3], THIIA Ao 
Aloe. B. Aaexcandep 


PajWasIbHbIe CKOPOCTH 3Be3/, CeKTpoBoro THNa AO, BOJIM3H BepuInHbI TpaeKTOopuHH Comnua, 
COMOCTABJIAIOTCA C TaHTeCHI[HaJIBHbIMH CKOPOCTAMH 3Be3)I TOFO 2Ke THMa B 30HeE HeOa, OTKJIOHEHHOII 
Ha 90° OT BepLUIHHbI, UTO MIPHBOAMT K KHHeMaTHYeCKON cpejHei aOcOMOTHOM BeJIMUHHe B+0",14. 
[IpeqnomaraA, UTO Kaka 3Be3jja THA Ao oOsayaeT OHO H TOM xe AOCOJIOTHOM BeJINUMHO!, 
MOXKHO HalHiTH MpHONMSHTeJIbHbIe pacCTOAHHA H3 HaOJWOMAeCMbIX KaKYIWIMXCA BeJIMYHH. OTOT 
MeTO/, IIPHMeHAeTCA K 475 3Be3aM THMa AO, ACHOCTb KOTOPbIX MpeBocxogzUT 6"™,5, uToObI 
OJIYUMTL IKBATOPHAJIbHbIe COCTAaBJIAIOUIMe JIMHeMHOM CKOpOCTH M3 COOCTBeEHHBIX JIBH>KeHHit 
H pajiMasIbHbIx CKOpocteii. 

SuaueHuaA ckopocTtu CosHya M lapamMeTpoB 9JIMMCONya CKOpocTei BbIUHCJIAIOTCA MH COMO- 
CTaBJIAIOTCA C Pe3yJIbTAaTAaMM Jpyrux HCcIe7OBaHHiK. 

TanakTuyeckad jourota SOoubwio OCH MosyuaeTcA B 354°, a CpeyHMe KBafpaTHUHBIe coocT- 
BeHHbIX CKOpOCTeii B HallpaBJICHHH rlaBHbIxX Oceii COOTBETCTBCHHO 18,9 KM CeK., I1,I KM/CeK. H 
8,4 KM/CeK,. 


CHCTEMA CKOINIEHHM B MAJIOM MATEJIJIAHOBOM OBJIAKE 
FE. M. JIundcer 


BpiosHeHa cbemKa cKkomeHHi B Masom MarennaHopom oOnake MpH MOMOIM TesIecKoma 
ADH (Apmax-J[ancunk-XapsBapf) 4 mactuHoK 103a.O c BbIepxKKoi B 3C MMH. B oOnmlem 116 
TeJI PaCCMATPHBAaIOTCH Kak CKOMJICHHA B OOIaKe. CHCTeMa STHX CKOIMIeHHi pactiIpoctpaHseTcA 
yO Kpbuia OOsIaka. OmpesesieHbI KaxkylliMecH pa3smMepbI Ha 10 JIOMMOBHIxX (250 MM) MJIaCTHHKaXx 
Merxanpd. OOsnauHad cucTeMa Npoctupaetcs Ha 8000 Mc No HanOoubUIeH AMaroHasM (BKITOUaAA 
KPbluIO) MH Ha 5400 lic MO HaMMeHbuIeit WHaroHasm. Huuto He CBHeTCIECTBYeT O CIMpasIbHoil 
CTpyKType pacnipejlesIeHHA CKOMIeHHII. 


HAJIHGHE PACHIPEDTEJIEHHA B IIPOCTPAHCTBE PAJIMO-HCTOUHHKOB 
B PESYJIBTATE CbhEMOK HA UWACTOTE B 159 MITII 


Al. O. dnc, IT. A. T. Metiep u Joc. A. Metixwagm. 


PaccMOTpeHbI H€KOTOPbI€ pe3yJIbTaTbl CbhéEMKH PaHO-HCTOYHHKOB Ha “acTOTe B 159 MITIL, 
{WIM OPPaHHYeHHbIX YUaCTKOB HeOa, OTHOCHTEJIBHO paclipeeseHHA ITHX HCTOUHHKOB B MIpoct- 
PaHCTBe. 

ConoctaBiieHHe OTJ{CJIBHbIX HCTOUHHKOB MOTBeEprKMaeT PeaJIbHOCTh OKOJIO YeTBEPTH HCTO- 
UHHKOB, Halij[eHHbIX IPH MperxkKHei CbéMKe Ha “aCTOTe B 81,5 MITIL, HO HajlexKHOCTh Olpes[eseHHA 
ciaObIx MCTOUHHKOB Xy?Ke, 4YCM OHO ABJIATIOCb MO pexkHeMy cnucky. CpaBHHBad HOBbIe 
NOJOHKCHHA C NOMOHKCHHAMH, ONpeesIeHHbIMH Musiipcom 4 Com Ha uactote B 85 MITU, HaiijqeHo, 
YTO MIOJI0xKCHHE HCTOUHHKOB, IWIOTHOCTh MOTOKa 9HCPrHH KOTOPbIX Ha uacTOTe B 85 MITIL MpeBo- 
CXOJMT 20.10°*° BT.m *. riy!, MOCTAaTOUHO COPacyIOTCA, HO COrslacOBaHHe HeCyJOBJIETBOPHTeIbHO 
Juin Oonee cnabbix HCTOUHHKOB. 

Ilocrpoenne rpacpuka log N B 3aBucumoctTn oT log S (rae N—uncs0 HCTOUHMKOR, TMJIOTHOCTh 
NOTOKa 9HEPrHH KOTOPbIx MpeBocxoguT S) AaeT HAKIIOH B —2,2 JIA OHOrO yuacTKa HW —2,7 JIA 
“Apyroro. Ho 9TH HaKJIOHbI MOryT ObITb HCKaxKeHbI JO HeH3BeCTHOH cTeneHH, OmarogapxA 
HeACHOCTH M@?KJIy COCeJIHHMH HCTOUHHKaMH A OylarofapA TOMY, UTO YrOBbIe MaMeTPbI MCTO- 
4YHHKOB KOHCYUHBI,. 

AHasIM3 TPH MOMOUIM MeTO/a, 3aBHCALIerO OT CTaTHCTMYeCKOrO paciipeeseHuA aMIVIMTYyJ 
Ha CHHMKax, HO He 3aBHCHUlero OT OUIMOOK M3-38 HeACHOCTH, yKa3bIBaeT, UTO HaOOeHHA 
dosee CnaGbIX HCTOUHHKOB He MpPOTHBOPeYMT FHMOTEse OAHOPOAHOTO paciipe;esIeHHA HCTOUHHKOB. 
Ocraetca mpoTHBOpeune M@KJIy 9THM OHOPOHbIM pacipeseseHueM MH UHCJIOM CHJIBHBIX 
HCTOUHHKOB. 

MookHo WOoKa3aTb, UTO HaGJOe€HHA Ha YACTOTe B 159 MITL, COMOCTaBJICHHbIe C HaOJOeCHHAMH 
Ha yacToTe B 81,5 MITI, IpoTHBOpeuaT rHMOTese OMHOPOAHOTO paciipeyesweHHA. OTOT pesysbTaT, 
OHaKO, OCHOBaH Ha JI@HHbIX, CTATHCTHUeCKOe€ 3HAYCHHE KOTOPbIX HEOCTATOUHO, UTOObI aTb 
HCUePNbIBaIOUIHe BLIBOJIbI M HEOOXOAMMbI JaibHeuimMe ACCIesOBaHHA, 





